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This paper reviews progress on the Fractional Quantum Hall Effect (FQHE) based on what we 
term hamiltonian theories, i.e., theories that proceed from the microscopic electronic hamiltonian 
to the final solution via a sequence of transformations and approximations, either in the hamil- 
tonian or path integral approach, as compared to theories based on exact diagonalization or trial 
wavefunctions. We focus on the Chern-Simons (CS) approach in which electrons are converted to 
CS fermions or bosons that carry along flux tubes and our Extended Hamiltonian Theory (EHT) in 
which electrons are paired with pseudo- vortices to form composite fermions (CF) whose properties 
are a lot closer to the ultimate low-energy quasiparticles. We address a variety of qualitative and 
quantitative questions: In what sense do electrons really bind to vortices? What is the internal 
structure of the Composite Fermion and what does it mean? What exactly is the dipole picture? 
What degree of freedom carries the Hall current when the quasiparticles are localized or neutral or 
both? How exactly is the kinetic energy quenched in the lowest Landau level and resurrected by 
interactions? How well does the CF picture work at and near v = 1/2? Is the system compressible 
at v = 1/2? If so, how can Composite Fermions be dipolar at v = 1/2 and still be compressible? 
How is compressibility demonstrated experimentally? How does the charge of the excitation get 
renormalized from that of the electron to that of the CF in an operator treatment? Why do 
Composite Fermions sometimes appear to be free when they are not? How are (approximate) 
transport gaps, zero-temperature magnetic transitions, the temperature dependent polarizations 
of gapped and gapless states, the NMR relaxation rate 1/Ti in gapless states, and gaps in inho- 
mogeneous states computed? It is seen that though the CS and EHT approaches agree whenever 
a comparison is possible, results that are transparent in one approach are typically opaque in the 
other, making them truly complementary. 
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I. INTRODUCTION 

Twenty years ago, von Klitzing, Dorda, and Pepper (1980) made the first discovery in what has proved to be a 
vital and exciting subfield of condensed matter physics to the present day, that of the quantum Hall effects (QHE). 
Their discovery that the Hall conductance a xy of a two-dimensional electron gas is quantized at integer multiples of 
e 2 /2irh, the quantum unit of conductance, is known as the integer quantum Hall effect (IQHE). Soon after, Tsui, 
Stormer, and Gossard (1982) discovered the even more puzzling fractional quantum Hall effect (FQHE). Since the 
first observation of the fraction |, many more have been seen (see, for example, Stormer, Tsui, and Gossard 1999). 

This review focuses on what we term the hamiltonian theories of the FQHE, by which we mean theoretical ap- 
proaches that begin with the microscopic hamiltonian for interacting electrons and try to obtain a satisfactory de- 
scription of the underlying physics through a sequence of transformations and approximations in the operator or path 
integral formalism. We bother to give a special name to what appears to be business as usual, because an alternate 
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approach, pioneered by Laughlin (1983a, b) and based on writing trial wavefunctions, has proven so extraordinarily 
successful. While we will of course discuss the wavefunction approach here because the hamiltonian approach is 
inspired by it, the discussion (of this and any other topic) will be aimed at serving our primary goal: Providing a 
cogent and critical description of the hamiltonian approach in one place, including all the hindsight and insight the 
intervening years have provided. Of necessity, many topics will have to be omitted or treated summarily. Fortunately 
many excellent reviews exist 1 and the reader is directed to them. 

The hamiltonian theories described here address both qualitative and quantitative issues. They give a concrete 
operator realization of many heuristic pictures that have been espoused and make precise under what conditions and 
in what sense these pictures are valid. They allow one to compute (within reasonable approximations) a large number 
of quantities at zero and nonzero temperatures at equal and unequal times. They open the door for a treatment of 
disorder. 

The hamiltonian theories themselves fall into two categories. The first, which we call the Chern-Simons (CS) 
approach, consists of making a singular gauge transformation on the electronic wavefunction that leads to a composite 
particle which is the union of an electron and some number of point flux tubes. The composite particles have 
a nondegenerate ground state at mean-field level. They are coupled to a gauge field that is fully determined by 
the particle coordinates. This formalism is best suited for computation of a variety of response functions. It is 
particularly effective at v = 1/2, where adherence to gauge invariance is of paramount importance if certain low- 
energy phenomena pertaining to the overdamped mode, coupling to surface acoustic waves, or the compressibility 
are to be properly described. On the other hand, the composite particles do not exhibit in any transparent way the 
quasiparticle properties (such as charge e* or effective mass to*) deduced from trial wavefunctions, though in principle 
they would surface after considerable work. It is also very hard to obtain a smooth limit as the electron mass to — ► 
in this approach. The Extended Hamiltonian Theory (EHT), which we developed over the years addresses some 
of these issues. 2 By the adjective extended we signify both that our work is an extension of older CS work and 
that our hamiltonian is defined in an extended or enlarged Hilbert space with additional degrees of freedom. The 
enlarged space allows us to introduce a quasiparticle that is a much better approximation to what we expect, to 
easily disentangle low-energy (Lowest Landau Level) and high-energy physics, to compute a variety of gaps, finite 
temperature properties, to describe inhomogeneous states and so on. On the other hand, the proper treatment of 
the additional degrees of freedom is very difficult to ensure in the computation of certain very low-energy quantities. 
Both the nonzero compressibility and an overdamped mode at v = 1/2 (and its experimental consequences) are all 
but invisible in this approach, though they can be extracted with some effort. No contradictions exist between the 
predictions of the CS and EHT approaches, which by and large tend to make predictions in complementary regions. 
Their predictions can be shown to agree in overlapping regions, but only with some effort. 



A. The experiment 

Figure 1 shows in schematic form the experiment in a rectangular geometry, with the current j x , magnetic field B z 
and electric field E y in mutually perpendicular directions. By definition 

jx &xyEy. (1) 

If we multiply both sides by the width of the sample we obtain 

Ix = O'xyVy (2) 

as the relation between current and voltage, reminding us that in this d = 2 problem the Hall conductance and 
conductivity are the same. 3 

Let us ask what one expects for a xy based on the simplest ideas, so as to place the experimental discoveries in 
perspective. Let us ignore interactions and disorder. Then we can assert that in the steady state, the electric and 
magnetic forces balance: 

eE = evB (3) 



1 The earliest comprehensive introduction to the full range of the quantum Hall effects is Prange and Girvin (1987). Huckestein (1995) 
concentrates on the integer quantum Hall effect. A more recent comprehensive review is Das Sarma and Pinczuk (1997), while a more 
focused treatment of the enormous body of research on Composite Fcrmions is summarized in Hcinonen (1998). 

2 We used to refer to this as THE hamiltonian formalism, but were persuaded to consider a name change that better reflected the state 
of affairs. Our compliance could not have been more total- the new acronym is the exact reverse of the old one. 

3 The longitudinal conductance and conductivity are related by the (aspect) ratio of length to width. 
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FIG. 1 The Hall experiment in a rectangular geometry. 



where e is the charge of the electrons, and v is their velocity. (Vector indices are omitted when obvious.) Therefore 

neE I A \ 

j = nev = -— (4) 

where n is the number density. Thus 

ne 
~B 



°*v = -5- ( 5 ) 



This result is unaffected by interactions and relies only on Galilean invariance (relativity). To see this, let us perform 
a boost to a frame in which the electric field (to leading order in v/c) 

E' = E - v x B (6) 

vanishes, that is, to a frame with v = E/B. In this frame j = 0. Boosting back to the lab frame we obtain j — neE/B 
and regain Eqn.JHJ). 

Figure 2 is a caricature of what is measured at T = in the thermodynamic limit. The key feature is that a xy , 
instead of varying linearly with ne/ B, is quantized at steps or plateaus. This article focuses on steps given by 

= V > (?) 

where 

V 



2ps + 1 



p= 1,2, s = 0,1,2,.. (8) 



The IQHE corresponds to s = 0, v — p, while FQHE refers to s > 0. 4 We explicitly display the h which will soon 
be set equal to unity. 

On each step the longitudinal conductances vanish. Thus the conductivity tensor assumes the form 

v 



Consequently the diagonal part of the resistivity tensor is also zero and the transport is dissipationless. In the real 
world, due to nonzero temperature T and finite sample size, the transitions between steps acquire a finite width, with 
d xx > therein. 



4 Fractions with denominator 2ps — 1 require only minor modifications and are left as an exercise. 



ne/B 



FIG. 2 Schematic form of measured Hall conductance at T = 0. The straight line is the Galilean invariant result, while the 
steps describe the data. 



Why is the quantization of conductance so surprising, given that we have seen so many instances of quantization 
of observables before? The answer is that the observable in this case does not refer to an atom or molecule, but to 
a macroscopic sample, with sample specific disorder. The role of disorder is truly paradoxical in the QHE: on the 
one hand, without disorder, we cannot escape the Galilean invariant relation a xy = ne/B; on the other hand, despite 
disorder (which surely varies from sample to sample), the conductance in any step is constant to better than a few 
parts in 10 10 ! A fairly detailed explanation of this now exists and will be described soon. An integral part of this 
resolution is also why a xx — 0, for only then can we understand why o~ xy is so well quantized even though the leads 
used to measure V y may not be perfectly aligned- there will be no unwanted contributions from the longitudinal 
voltage drop. 



B. What is special about the steps? 

In searching for the physics underlying any one step or plateau, it makes sense to begin with the points where the 
straight line (the Galilean invariant result) intercepts the steps. Not only are these points singled out by experiment, 
but they also have a Hall conductance that we could hope to understand without including disorder. 5 

At these points we have 

ne e 2 p 

~B = 2^h2ps + l ( ' 

which we rewrite as 

^m- 2s+l p 

Thus the ratio of B, the flux density, and n, the particle density, has some special rational values at these points. For 
example in the case of the fraction 1/3 there are three quanta of flux per electron. 



In the absence of disorder any state would have the desired Hall conductance. We are looking for a nontrivial correlated state that is 
robust under the inclusion of disorder and capable of dominating the plateau it belongs to. 
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To understand what is so special about these values, recall the following textbook results (Shankar 1994) about a 
single particle of mass m and charge e moving in two dimensions in a perpendicular magnetic field. 

• The energy is quantized into Landau Levels (LL's) located at E = (ehB /m)(n + i) = fiuj^in + i), where ujq is 
called the cyclotron frequency. 

• Each LL has a degeneracy equal to &/&o, the flux in units of the flux quantum $o = 2irh/e, or : 

Degeneracy per unit area of each LL = ^ . (12) 

• The wavefunctions of the Lowest Landau Level (LLL) are, in the symmetric gauge, 

V> LLL = z m e - |z|2/4 ' 2 z = x + iy m = 0,l,... (13) 



where I — ^JnjeB is the magnetic length. We will often drop the gaussian factor in LLL wavefunctions. Higher 
LL wavefunctions depend on both z and z. 

Thus Eqns. (|ll|l and Eqn. 1121 )imply that at the special points 

Number of states in a LL Flux quanta B 1 _ 1 

Number of particles particles n(2nh/e) p 

Note that v : which stood for the dimcnsionless conductance, is thus the number of occupied LL's. It is called the 
filling factor. In the FQHE v is a fraction, restricted in this article to be < ^ unless otherwise mentioned. 

If v < 1, there are more LLL states than particles, and in the noninteracting limit they can all be fit into the 
LLL with room to spare. For example at v = ^, there are 3 LLL states per electron. This macroscopic degeneracy 
in the noninteracting ground state means we cannot even get started with a perturbative treatment of disorder and 
interactions. This is the central problem. While this problem exists for any v that is not an integer, experiments 
suggest that v~ x = 2s + ~ is somehow preferred by nature. Indeed at such points there is a way out, but it is not 
simple perturbation theory. This approach is what this article is all about. 

An alternative to perturbation theory, the Hartree-Fock (HF) approximation, also does not work, if applied directly 
to the electronic hamiltonian. Let us see why. In HF, one takes the interaction (quartic in fermion operators) and 
obtains an expression quadratic in the fermion operators by replacing various bilinears by their ground-state averages, 
initially taken as free parameters, and dropping some higher-order fluctuations. The quadratic hamiltonian is then 
solved and the ground state energy is then minimized as a function of the assumed averages. At the minimum, the 
actual averages will self-consistently come out equal to the assumed averages. Here is a toy model that illustrates the 
main ideas. The model has just two Fermion operators c and d: 

H = e c c c + Ed d) d + no c'cd'd (15) 

Let us separate out the bilinears into averages (c'c) = A c and (cftd) = X d and fluctuating parts : c^c : and : d'd : 

Sc = : c'c: +(c f c) =: c f c : +A C (16) 

dU = :dU: +(dU) =: Sd : +X d . (17) 

We now rewrite H as 

H = e c Sc + e d d^d + u [Ajc+c + \ c d) d - X c \ d + : c f c :: d) d :] (18) 

and neglect the last term, quadratic in the fluctuations, to obtain the HF hamiltonian. 

Let Eq(A c , X d ) be the energy of the ground state |0) of the HF hamiltonian. From the Feynman- Hellman theorem 

which you can prove using the fact that (0|0) = 1 has no A derivative, it follows that at the minimum, A c = {c'c) and 
\ d = (dU) . 

Why does this fail in the FQHE problem? For integer filling this procedure gives the expected state of an integer 
number of LL's filled (at least for not too strong interactions). However, it turns out that for generic fractional filling 
translational symmetry is spontaneously broken in HF, and the solution corresponds to a crystalline (usually a Wigner 
Crystal) state. Due to the breaking of translational symmetry, the LL and angular momentum indices (n, m) are 
no longer good quantum numbers. The single-particle states are now bands in the crystal, parametrized by a Bloch 
quasimomentum k and a band index rifc. The HF hamiltonian for FQHE will of the form H = X^riik 6 ™^^) c n i ,(^ ; ) c nb(l <; ) 
up to c- number terms (the analogs of A c A a ). The self-consistent solutions were worked out by Yoshioka and Fukuyama 
(1979), Fukuyama and Platzman (1982), and Yoshioka and Lee(1982). The HF solution is compressible, translationally 
non-invariant, and has no preference for any particular density. Thus, it does not describe the FQHE phenomenology. 



7 



C. The IQHE - A warm-up 

We begin with a brief look at the IQHE case s — 0, or v = p, which paves the way for Jain's (1989) view of FQHE 
as the IQHE of entities called Composite Fermions. 

In the IQHE, exactly p LL's arc filled in the noninteracting limit. There is exactly one totally antisymmetric 
ground state which we denote by \ p . The simplest example is p — 1, with just the LLL filled up. The corresponding 
wavefunction is 

Xi = Det 



Gaussian — \\{zi — Zj) exp 



i<j 



-E 



4/ 2 



(20) 



This nondegenerate ground state is separated from excited states by a gap equal to the cyclotron energy uj . Though 
now we have a starting point for perturbation theory, we still need to actually carry out the perturbative calculations 
and in particular understand why the Hall conductance is unaffected by these perturbations and stays at the Galilean 
invariant value. We also need to understand why the conductance is unchanged as we make small changes in density, 
i.e., why there are steps. 

The explanation of the IQHE, discovered over the years, can be most easily understood in the noninteracting limit. 
The single-electron problem with certain special types of disorder can be exactly solved (Prange, 1981, Aoki and 
Ando, 1981, Prange and Joynt, 1982), and the exact solution shows that the Hall current is independent of disorder. 
Moving beyond this, Trugman (1983) showed in a seminal paper that for generic disorder in very high magnetic fields 
the electronic states follow the equipotentials of the disorder potential, and each LL gets broadened into a band. The 
problem has many fruitful analogies with percolation. In particular, there is exactly one energy at the band center at 
which extended states can exist, while states at all other energies are localized. When the chemical potential n crosses 
this energy, the Hall conductance changes by an integer. Contrariwise, when n lies in energy range corresponding to 
localized states, changes of a (changes of filling) produce no change in a xy , which explains the steps. 

From a somewhat different point of view, Laughlin (1981) and Halperin (1982) showed by a gauge argument that 
to have a quantized Hall conductance in the noninteracting limit all one needs is a xx = 0, which is assured if the 
chemical potential lies in the region of localized states. This argument shows that the Hall conductance (in units 
of the quantum unit of conductance e 2 /2irh) is a topological integer. There is also a field-theoretic analysis of the 
noninteracting quantum Hall problem in the presence of disorder which yields similar results (Pruisken, 1984, 1985a, b, 
and Levine, Libby, and Pruisken, 1983, 1984a, b,c). 

The IQHE can also be understood in the other limit, where interactions dominate, as follows: In this case if one 
is not exactly at v — integer, the excess (or deficit of) particles form a Wigner Crystal due to interactions. First 
turn off the interactions and impose an external periodic potential of the same period as the Wigner Crystal. In this 
case Thouless, Kohmoto, Nightingale, and den Nijs (1982) showed that the dimcnsionlcss Hall conductance has to be 
an integer, and that this integer is topological. Later work by Thouless (1983), Niu and Thouless (1984), and Niu, 
Thouless, Wu (1985) showed that since this integer is topological, it is robust to adiabatically introducing disorder and 
interactions. As long as the charge gap does not close, and there are no ground state transitions, the Hall conductance 
will remain unchanged 6 . This allows one to turn off the external periodic potential as the interactions are being turned 
on. For further details on the IQHE the reader is referred to an excellent review by Huckestein (1995). 



D. The FQHE 

We turn now to fractions, such as |. There is nothing special about disordered single-particle states at such a 
filling: They are expected to be localized except right near the middle of the Landau band (Trugman, 1983). If 
we ignore disorder and interactions, we are left with the macroscopically degenerate manifold of many-body states. 
Presumably interactions will select the ground state. How is one to find it? 

Here the story branches into two trails. The one blazed by Laughlin (1983a, b) consists of writing down inspired 
trial wavefunctions; the other is the hamiltonian approach, which starts with an assault on the degeneracy problem. 



By definition, an integer cannot change continuously. 



8 



1. Laughlin's answer 

After some experimentation, Laughlin wrote down the following celebrated trial wavefunction for v = I /{2s + 1) : 

= - ^) 2s+i cx p(- e i^i 2 /4i 2 ) (2i) 

j<i i 

which was shown to have nearly unit overlap with the numerical solution for small systems for generic repulsive 
interactions (Laughlin, 1983a, b, 1987). We will refer to v — l/(2s + 1) as Laughlin fractions, as per convention. 

Of the many remarkable properties of this function, we list those that have the greatest bearing on what follows. 
First, it is from the LLL and obeys the Pauli principle under particle exchange, since 2s + 1 is odd and spin is assumed 
polarized. Halperin (1983,1984) zeroed in on one of its central features: it has no wasted zeros, by which he meant 
the following. Consider tp(zi), which is ^(21, . . . , Zn) as a function of any one variable, randomly chosen to be z\, 
with all others held fixed. Given that the sample is penetrated by N/v — N(2s + 1) quanta of flux, the phase of 
ip(zi) has to have an Aharonov-Bohm phase change of 2ttN(2s + 1) per particle, or N(2s + 1) zeros given the LLL 
condition of analyticity. Only N of these had to lie on other electrons by the Pauli principle. But they all do lie on 
other electrons, thereby keeping the electrons away from each other very effectively, producing a low potential energy. 
(The kinetic energy is of course the same for any function in the LLL.) 

By showing that quantum averages in these ground states are precisely statistical averages in a one-component 
plasma (Baus and Hansen, 1980, Caillol, Levesque, Weis, and Hansen, 1982), Laughlin (1983a, b) showed that the 
system is an incompressible fluid, which means a fluid that abhors density changes. Unlike a Fermi gas which increases 
(decreases) its density globally when compressed (decompressed), an incompressible fluid is wedded to a certain density 
and will first show no response to any applied pressure, and then suddenly nucleate a localized region of different 
density (just the way a Type II superconductor, in which a magnetic field is not welcome, will allow it to enter in 
quantized units in a region that turns normal). 

Laughlin (1983a, b) also provided the wavefunction for a state with such a localized charge deficit. If one imagines 
inserting a tiny solenoid at a point zq and slowly increasing the flux to one quantum, one must, by gauge invariance, 
return to an eigenstate of H, and each particle must undergo a 27r phase shift as it goes around z . 8 This condition 
and analyticity point to the ansatz 

= \{{z l - Z )*2.s+1 (22) 

i 

This is a quasihole. There is a more complicated state with a quasiparticle. 

The prefactor is a vortex at zq. Since vortices play a significant role in the FQHE let us digress to understand them 
better. Consider first a real flux tube inserted into the sample at a point z . Clearly every electron will "see" this 
tube at zq, i.e., ip{zi), the wavefunction, seen as a function of any one coordinate, chosen to be z\, will see a 2tt phase 
shift as z\ goes around zq . The flux tube clearly has a reality and location of its own, independent of the locations of 
the electrons. The vortex is an analytic LLL version of the flux tube: the 2ir phase change is accompanied by a zero, 
i.e., there is an analytic zero at z for every coordinate. The zeros associated with the vortex have the key feature 
that their location does not depend on the location of the particles, i.e., it too has an independent reality and location 
like the flux tube. 

There are also vortices in Laughlin's ground state wavefunction. The only difference is that instead of sitting at 
some point zq, they are anchored to the particles themselves. To see this, consider v = |, and focus on the part 
involving just Z\\ 

n^-^i) 3 (23) 

If we freeze z\ and view it as a parameter like zq, we see that there is a triple vortex on particle 1, and by symmetry, 
each particle. Apart from the vortex mandated by the Pauli principle, there are two vortices bound to each electron. 
(The Pauli zero comes with the turf, for being a fermion, and is not included in the count of vortices attached to 
enforce correlations.) The term vortex is again appropriate here since the location of the vortex is independent of any 
coordinate except, of course, the electron to which it is attached. 

Contrast the zeros that constitute a vortex to the zeros of a generic antisymmetric analytic polynomial of the same 
degree as the Laughlin wavefunction. Once again, in tp(z\), there is one zero at Z{ (i > 1) independent of all other 



7 The discussion is independent of this choice. 

8 The flux must be inserted slowly enough to prevent transitions to other excited states, but fast enough to prevent a lapse to the original 
ground state at the end, a possibility that exists in finite systems with disorder. 
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Zj by the Pauli principle. This zero is part of a vortex, since all particles will "see" (in the wavefunctions ip) the 
Pauli zeros anchored on the other particles. The non-Pauli zeros of tp(zi), by contrast, will be parametric functions 
of all other particle coordinates i — 2, ..TV (what else can they be?). If we now consider ip( z 2) its non-Pauli zeros will 
in general bear no relation to those of ?/>(zi). Since the location of the non-Pauli zeros move parametrically with all 
the z's, they do not belong to or form vortices. The reader must remember these distinctions in order to follow the 
discussion on the internal structure of CF's that follows shortly 9 . 

Returning now to the quasihole, another important property is that it represents a charge deficit of l/(2,s + 1) in 
electronic units. One way to show this is to employ the plasma analogy (Laughlin 1983a, b, 1987). Here is a more 
general way (Su and Schricffcr, 1981, Su, 1984) that only depends on the state being gapped, incompressible, and 
having a quantized Hall conductance. As the flux quantum $ is adiabatically inserted to create the quasihole, the 
charge driven out to infinity is given by integrating the radial current density j produced by the Hall response to the 
induced azimuthal E field 



For non-Laughlin fractions, the charge driven out by inserting one flux quantum is that of p quasiparticles, each of 
which has a charge equal to l/(2ps + 1) (Su, 1984). Note that the fractional a xy is the cause behind the fractional 
charge. The fractional charge of the quasiparticles has been confirmed experimentally (Goldman and Su, 1995, 
Goldman, 1996). The quasiparticles also have fractional statistics, as was pointed out by Halpcrin (1984), and 
explicitly confirmed in a wavefunction calculation by Arovas, Schrieffer, Wilczek, and Zee (1985) 10 . 

The insensitivity of the Hall conductance to disorder in the FQHE can be established by an extension of the 
Laughlin-Halperin gauge arguments for the IQHE. To show this in toroidal geometry, the old arguments have to be 
supplemented by the fact that at fractions of the form p/q, there are q ground states, and the assumption that upon 
adding q flux quanta one returns to the starting state. There seems to be a consensus that Hall conductance is a 
ground state property that is impervious to a reasonable amount of disorder, very much the way superfluidity is. 

Laughlin (1983a, b, 1987) explained the plateaus as we move off the magic fraction as follows. Suppose we move 
slightly off 3, to say v = 3 ± .01. The system has two choices. It can either go to a FQHE ground state at this 
fraction, or it can be at | plus some number of quasiparticles (holes) of charge ±|. It turns out the latter has a 
lower energy. The reason o~ xy is locked at 3 is that the quasiparticles/quasiholes, being sharply defined entities, get 
localized by the disorder potential and do not contribute to the Hall current. Thus the Hall conductance is expected 11 
to be unaffected by disorder. As we move further off |, the system switches to another ground state with its own 
quasiparticles, and o~ xy jumps to the neighboring plateau. It is clear that if the fate of the quasiparticles/quasiholes is 
to get localized, their fractional charge cannot explain the fractional Hall conductance. As seen earlier, the connection 
is the other way around. 

Thus a fairly complete description of the FQHE emerged from the Laughlin wavefunctions for the ground and 
quasiparticle/quasihole states. However numerous other issues surfaced and were tackled by subsequent work in the 
intervening two decades. We now steer our attention to those bearing directly on this paper. 

Sticking for a moment to the Laughlin fractions, one important question that was raised was whether the system 
was really gapped, i.e., could there be lighter excitations than Laughlin's quasiparticles and quasiholcs? Perhaps a 
quasiparticle and quasihole could bind to form a gapless excitation. Haldane and Rezayi (1985) verified by exact 
diagonalization of small systems that the u = 3 system was gapped. Girvin, MacDonald and Platzmann (1986) 
explored the question of neutral particle-hole excitations using an analogy with Feynman's work on superfluids and 
found a way to calculate the dispersion relations of the magnetoexciton within the LLL and showed that they were 
gapped in this sector for v = ^. 



If the Laughlin function is perturbed slightly, 1/1(21) will have one Pauli zero on every other electron and two others nearby. As long as 
the perturbation is small we can say, at length scales bigger than the excursion of the zeros, that electrons are bound to double vortices. 
For stronger perturbations that take us far from the Laughlin ansatz, there will be no simple description of correlations in terms of 
vortices. 

Fractional statistics was shown to be an emergent property of Composite Fcrmions by Goldhaber and Jain (1995). 

There is currently no way of rigorously calculating transport coefficients in the FQHE in the presence of disorder, hence the qualifier 
" expected" . 




(24) 



2irk e 2 
e 2irh 



(25) 
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2. Jain's Composite Fermions 

Let us now proceed to the non-Laughlin fractions. Here the central question is this: Given that non-Laughlin 
fractions like | are seen in experiment, what is their wavefunction? Simply replacing the factors (zi — z 3 ) 2s+l by 
(zi — Zj) 5 / 2 in Laughlin's ansatz is not acceptable since this does not produce the mandatory change of sign under 
particle exchange. In the hierarchy approach (Haldane, 1983, Halperin, 1983, 1984) the quasiparticles of the Laughlin 
states condense into their own FQHE states, whose quasiparticles in turn do the same thing and so on. While this 
approach gives a natural way to generate additional fractions, it does not give explicit wave functions in terms of 
electrons. 

Explicit trial states for fractions v = 2p s+i were provided by Jain (1989,1990,1994), who also explained why they 
were natural in terms of objects called Composite Fermions (CF's). 

Jain's scheme hinges upon the seminal idea of flux attachment, which plays such a central role in the FQHE that 
it merits a little digression. It was introduced first by Leinaas and Myrrheim (1977) in terms of wavefunctions, and 
explored in the language of Chern-Simons field theories by Wilczek (1982a, b). The crux of the idea is that in d = 2 one 
can have particles (dubbed anyons by Wilczek) that suffer a phase change e tS upon exchange, with 8 = , ir ( mod 2tt) 
corresponding to bosons and fermions respectively. To obtain one of these particles one takes a fermion and drives 
through its center a point flux tube, the amount of flux being decided by the desired 9. In particular, if this tube 
contains an even/odd number of flux quanta, the composite particle one gets is a fermion/boson. 

Jain (1989,1990) exploited flux attachment for v~ x = 2s + -, as follows. Suppose we trade our electrons for CF's 
carrying 2s point flux quanta pointing opposite to the external B. On average the CF's effectively see ^ flux quanta 
per particle and fill up exactly p LL's. At the mean field level, this approach gives the following trial wavefunction: 



p/2ps+l 



n 



~(Zi - Zj)~ 



Zi - Zi 



■X P (z,z) 



(26) 



where \y is the CF wavefunction with p- filled LL's and the prefactor takes the CF wavefunction back to the electronic 
wavefunction. (This will be made clearer shortly.) 

Jain improved this ansatz in two ways and proposed: 



i> P /2 P 8+i = V Y[(z t - z 3 ) 2s ■ x P {z, z) 



where 



is called the Jastrow factor, and 



i<3 



i<j 



\2s 



(27) 



(28) 



V-.z-.2V 



dz 



(29) 



projects the wavefunction to the LLL by eliminating all z's in x P by turning them to ^-derivatives on the rest of the 
wavefunction, except for the gaussian. (We shall see 12 that in the LLL [z,z] — —2l 2 .) 



In making the change \\ 



(fj-fl) 



— > Jli<j ( z i — z j) , Jain replaces flux tubes by vortices. Although in both 
cases each particle picks up an extra phase shift of 47rs on going around another, only the vortex factor has the desirable 
multiple zero that keeps the particles apart and lies in the LLL. Jain does not need to justify these modifications 
since in writing a trial wavefunction, he is free to resort to any changes that improve its energy. However, in the 
hamiltonian approach, the replacement of flux tubes by vortices and the projection to the LLL prove difficult but not 
insurmountable. 

At p = 1, since \( z , z) — Yl i<3 X z i — z i) ' Gaussian, we do not need V to get back Laughlin's answer. For p > 1 we 
have concrete expressions for ^/ in terms of electron coordinates, although the action of V can be quite involved, and 
has a big impact on the wavefunction as we will see soon. 

Thus while the degeneracy of the noninteracting problem is present for any v < 1, at the Jain fractions one can 
beat it by thinking in terms of Composite Fermions. As we move off the Jain fractions, the incremental CF's (particles 



For a nice review of states, operators, and matrix elements within the LLL see Girvin and Jach (1984) 
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or holes) get localized, giving rise to the plateaus. This is the sense in which the CF approach allows one to think of 
the FQHE in terms of the IQHE of the CF's. Since the IQHE can be understood without invoking interactions, it is 
sometimes suggested that CF's are free. Later we will see why this cannot be so. 

The reader will have noted that both Laughlin and Jain wavefunctions make no reference to the inter- electron 
potential. This feature, which permits them to work for a whole class of potentials also renders them insensitive to 
specific features, including even the range of the interaction. 



3. Vortices, zeros, and the dipole picture re-examined 



Let us return to the vortices that are attached to electrons in the Jastrow factor. We have already discussed 
Halperin's observation (Halperin 1983, 1984) that electrons are bound to vortices in the Laughlin wavefunction. This 
concept was taken up and vigorously pursued by Read (1989). In particular Read made a clear distinction between 
flux attachment and vortex attachment, which were often loosely interchanged. Attaching electrons to flux tubes is a 
mathematical trick; flux tubes are unphysical and neutral. Vortices on the other hand are physical excitations of the 
Laughlin ground states, and electrons would naturally bind to them since the two are oppositely charged. The reason 
the CF sees a weaker field is not due to any mysterious capture of flux tubes. Indeed, the external field is uniform 
and not quantized into point flux tubes. What really happens in the Laughlin case is that each electron pairs with 
2s vortices and when the composite object goes around on a loop, it sees a phase change of — 2tt(2s + 1) per enclosed 
particle due to the external B and 27r(2s) per enclosed CF due to the vortices attached to each one. 

Next, since the 2s-fold vortex has a charge — 2s/(2s + 1) as per the flux insertion argument (Su, 1984), Eqn. I|25|l . 
the vortices reduce e down to the quasiparticle or CF charge 



e* = 1 - = (30) 

2s + 1 2s + 1 V ' 

In a Fermi liquid an added electron is also screened by a correlation hole; here the difference is that the possible 
response of the Laughlin liquid to an extra electron is limited to an integer number of vortices of quantized charge. 

Read (1994, 1996) next applied the notion of electrons binding to vortices to the case v = ^ and derived what 
is called the dipole picture. We now turn to a critical review of this concept. Since the dipole picture comes from 
the projection involved in the v = ^ wavefunction, it is instructive to start with understanding the impact of V in a 
general Jain fraction where higher CF-LL's {p > 1) are filled. 

Consider N electrons at v = |. Prior to projection, each electron in the Jain wavefunction has a double vortex 
sitting on it due to the Jastrow factor. As we explained earlier, the term vortex is merited since every other particle 
sees it at the same place. At this stage it is correct to view the CF as an electron plus a double vortex 13 . All this 
changes after projection to the LLL. Now there can only be 5N/2 zeros of the wavefunction as a function of, say, 
z\] A of them must lie on other electrons by the Pauli principle (forming single vortices anchored to electrons), with 
3 N/2 left over. Clearly they cannot all be on electrons, or be associated with them uniquely (since we have 3/2 zeros 
per electron). If the zeros are not on electrons, their location must depend parametrically on the locations of all the 
electrons and they cannot organize themselves into vortices. 

Thus the wavefunction analysis leaves us with the following quandry. The CF cannot be viewed as an electron- 
vortex complex in the projected non-Laughlin states, but the quasiparticle charge, e* = 1 — 2 ps+i 1S r °bust under V 
since e* is tied to o~ xy which is presumed robust under projection. What, if any, is the entity that binds to the electron 
to bring e down to e* ? How does this entity enter the theory? What makes it bind to the electron? We will see that 
the EHT answers such questions. 

Let us now return to v = ^. In the unprojected wavefunction the double- vortex charge fully cancels the electron 
charge so that e* = 0. Read (1994) has argued that the neutral CF of momentum k has a dipole moment d* = kl 2 , 
based the wave function at v — i, also called the Rezayi-Read (Rezayi and Read, 1994) wavefunction, in which \v m 
Eqn. (|27(1 is replaced by the filled Fermi sea: 



*i ='P\\{z l ~ z 3 ) 2 Det\e l ^^\. (31) 



i<j 

He considers 



There is also another zero, generally non-analytic, in \v which is antisymmetric. This is not part of the vortex count. 
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e ikr = cxp - (kz + kz) k = k x +ik y . (32) 

Since under the action of V, z acts on the analytic part of the wavefunction as 2l 2 d/dz, e ll2k & causes the shift 
z —* z + ikl 2 in the Jastrow factor. This motion of the vortex off the electron produces the dipole moment d* = kl 2 . 
The energy needed to separate the vortex from the electron (the Coulomb attraction), must begin with a term 
quadratic in the separation or momentum. This then gives the LLL fermions a kinetic energy or an effective mass m* 
resulting from interactions. 

Attractive though this picture is, closer scrutiny reveals that the preceding line of reasoning is incorrect in the 
following ways: 

• Since every Zi gets translated, (zi — Zj) 2s — > (z^ — Zj + iki — ikj) 2s and particle i sees a multiple zero associated 
with Zj at Zj + i(kj — kj). We cannot call this a vortex since the location of the zero varies with the label i. 
Different particles see it in different places. This problem is actually moot because of the next, more serious 
one. 

• Even this multiple zero is there for one particular assignment of k 's, or one term in the determinant. Upon 
antisymmetrization, we cannot relate the zeros of the sum over the N\ terms to the zeros of individual terms. The 
situation would have been different if we had been talking about poles which can survive such a sum. In fact, one 
of the two zeros that moved off the particles (in each individual term) should return upon antisymmetrization 
to lie exactly on the particles by Pauli's principle, leaving another zero to depend parametrically on all other 
coordinates, and all the fcj. Thus there is no reason to expect any simple relation between the location of the 
electrons and the non-Pauli zeros in "fi or to conclude that these zeros form vortices. We commend to the 

2 

reader experimentation with the limited but nonetheless instructive case of just two particles, to see some of 
these ideas in concrete form. 

• If one looks at v — | in isolation, the fact there is one non-Pauli zero per electron (in the thermodynamic 
limit) may tempt one to suggest that perhaps in this case they organize themselves into vortices. Besides our 
analysis of the gapped fractions which shows that this is extremely unlikely, there is the general argument that 
zeros that do not lie on particles or external flux tubes (i.e., locations with an independent reality) must vary 
parametrically with all coordinates and therefore cannot be organized into vortices. 

Thus at v — i, there seem to be two choices: Either use the unprojected wavefunction in which the Jastrow factor 
explicitly has two vortices per electron, but the dipole moment is zero (since the vortices are on the electrons), or go 
the LLL where the vortices disintegrate into a smaller number of ordinary zeros, not correlated with the electrons in 
any simple fashion. So what happens to the dipole? 

Our EHT, which provides an operator realization of the CF, will show that the dipole picture is more robust 
than the wavefunction-based arguments pointing to it. However one must look for it not in the wavefunction, but in 
operators and correlation functions. Even then, the value of d* will be sensitive to the details of the wavefunction. 
In this respect, it is quite unlike a xy which is robust under changes in the ground state wavefunctions, (including 
projection to the LLL), and depends only on the filling factor. 

E. Hamiltonian approach 

Let us begin by reviewing the challenges facing any hamiltonian theory. In the FQHE there is more than enough 
room in the LLL for all the electrons, leading to a macroscopic number of degenerate ground states in the absence of 
interactions. This preempts many of the standard approximations. It is expected that when interactions are turned 
on, there will emerge from this degenerate manifold a unique ground state, separated by a gap from other low-lying 
excited states. While the true ground state and low-lying excitations will of course contain an admixture of higher 
LL's, one expects that in the limit of bare mass m — > 0, or wo ^ oo, there will emerge a low-energy sector spanned 
by LLL states. In other words, in the limit m — > a nonsingular low-energy theory must emerge. The m-independent 
Laughlin wavefunctions for the ground state, quasiparticles and quasiholes illustrate this point. 

If all memory of m is lost, the only energy scale is set by the interactions. How is one to isolate the LLL physics, 
staring with the full Hilbcrt space? How is one to battle the degeneracy of the noninteracting ground state? How is 
one to get rid of the 1/m dependence of the kinetic energy? These are problems for the hamiltonian approach. 

Next, there is one vestige of m-dependence even in the LLL, discovered by Simon, Stern and Halperin (1996): if 
the external field B varies slowly in space, the zero point energy eB/2m, can no longer be eliminated by redefinition 
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of the zero. In fact the particles behave as if they have a magnetic moment /i* = e/2m coupled to B. The theory 
must reproduce this moment. 14 

We also do not want to ban higher LL's too soon, since the Hall conductance necessarily involves them 15 . Finally, 
Kohn's theorem (Kohn, 1961) assures us that despite interactions, the cyclotron mode at q = must be at u>o = eB/m 
and it must saturate certain sum rules at small q. The theory must pass this test. 

The hamiltonian theory must succeed in keeping the m- dependence where it belongs and exorcise it elsewhere. 

The next set of issues concern the quasiparticles, the CF's. A heuristic picture that we saw arising from the 
unprojected wavefunctions is that the quasiparticles are composites of electrons and 2s vortices. How is this actually 
realized in the hamiltonian approach? After all, this is not analogous to the statement that mesons are made of 
quarks and antiquarks, for unlike antiquarks, which appear in the hamiltonian, vortices do not. Vortices are not 
independent of electrons and to speak of them and electrons at the same time surely paves the way for overcounting. 
In the non-Laughlin cases, the issue is further complicated by the fact that there are not enough zeros to form 2s 
vortices per electron, and in any case the zeros are not organized into vortices. Yet some object of the same charge 
as the 2s-fold vortex seems to bind with the electron since in the end e gets reduced to e*. Somehow this object has 
to enter the theory and give rise to the CF. 

Other issues surround the CF. How strong are interactions between CF's? Why do they sometimes appear to be 
free when we can give persuasive arguments for why they cannot be? How is disorder to be included? 

Besides addressing these questions of principle, we need answers for quantitative questions. For example, what is 
the gap at v = i? Now, very precise answers can be given for such a question in the trial wavefunction approach 16 
or by exact diagonalization . While these results are founded on the microscopic hamiltonian, the intermediate 
steps are computer-intensive. It would be nice to have a theory which displays transparently the features of the CF 
deduced from the study of wave functions, and furnishes quantitative results, to say 10% accuracy. The same goes 
for polarizations in various states, and transitions between them. 

A specially interesting set of questions appear at and near v = 1/2. Can the notion of the CF survive in this 
region without a robust gap? In particular, will they really behave like particles of charge e* in this region? Since the 
average effective magnetic field B* = at v — 1/2, we have fermions in zero magnetic field, except for fluctuations. 
Will this Fermi system form a Fermi liquid after including fluctuations? What are its response functions? Will it be 
compressible? If so, how do particles with e* = manage to be compressible or even have a nonzero Hall conductance? 
What are the collective modes and how can they be detected? How is disorder to be handled? Exact diagonalization 
and the wavefunction based approaches are not very helpful for unequal-time correlations, and finite-temperature 
physics. Hamiltonian theories should fill this void and answer all the questions raised above. 



F. Organization of this review 

Section II introduces the reader to some notation and the physics of the LLL, including Kohn's theorem which 
is often referred to in our work. Section III describes the Chern-Simons theory of flux attachment including both 
composite bosons and fermions. Section IV deals with the region at and near v = 1/2 and is mainly about the work 
of Halpcrin, Lee and Read (1993). Section V explains our extended hamiltonian approach, EHT, formulated in an 
enlarged Hilbert space. It is shown that there are two approaches to solving our final equations: the conserving 
approximation in which the constraints are respected, but the CF physics is hidden, and a shortcut in which the 
opposite is true. The compressibility paradox at v = 1/2 is discussed. Section VI illustrates the conserving approx- 
imation in the EHT via the calculation of the structure factor at small q in the Jain series and the magnetoexciton 
dispersion relations for 1/3 and 2/3. Section VII is devoted to the computation of gaps and comparison to numerical 
work and experiment. Sections VIII deals with magnetic phenomena at T = for gapped and gapless systems. The 
question of whether CF's are free or not is discussed and answered in the negative, along with an explanation of why 
sometimes they appear to be so. Section IX addresses physics at T > 0, describing the computation of polarization 
and relaxation at v = ^ and polarization at | and |. The results are then compared to experiment. Section X deals 
with inhomogeneous states. Section XI gives a critical evaluation of the EHT. A summary and a discussion of open 
problems follow in Section XII. 



This magnetic moment is actually a compact way to describe an effect that is orbital in origin. 

The current operator goes as 1/m while a xy does not. The 1/m in the current is cancelled by the energy denominator of order u)q that 
comes from virtual transitions to the higher LL's. For a more detailed explanation, sec Girvin, MacDonald, and Platzman (1986), and 
Sondhi and Kivelson (1992). 

For a review of this method with complete references, see Jain and Kamilla (1998). 

Due to limitations of space we present only the earliest references: Girvin and Jach (1983), Haldane and Rezayi (1985), Su (1985), 
Haldanc (1985), Yoshioka (1986), Morf and Halperin (1986, 1987), with the exception of Morf et al (2002). 
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II. PRELIMINARIES AND NOTATION 

This section focuses on a single electron in two dimensions. The hamiltonian is 

Here A is the vector potential that leads to the external magnetic field B = V x A = — Bz, p is the canonical 
momentum of the electron, — e is its charge, and m is its bare or band mass. Note that B points along the negative 
z-axis. 

Although one needs to pick a gauge for A in order to find the wavefunctions, we can obtain the spectrum without 
making that choice. Let us define a cyclotron coordinate 

77 = Z 2 z x n (34) 



where I = y/h/eB is the magnetic length. Despite the name, the two components of rj arc not commuting but 
canonically conjugate variables: 



ft follows that 



describes a harmonic oscillator with energies 



[Vx ,Vy] = H 2 - (35) 



E = (n + i)w (37) 



where n is the LL index. 

The LL's are highly degenerate because another conjugate pair that commutes with rj, and called the guiding center 
coordinate 

R = r- ij (38) 

is cyclic. The components of R obey 

[R x , R y ] = -il 2 . (39) 

Thus I 2 plays the role of h. Since in the LLL (r; 2 ) = I 2 , R roams all over the sample, whose area L 2 plays the role of 
phase space. The degeneracy of the LLL (or any LL) is, from the Bohr-Sommerfeld quantization rule, 

D-— — - - — 

where $o is the flux quantum. This leads to the following result worthy of committing to memory: 

LLL states D flux quanta density eB 



particles N particle density 2imh 
At points where steps in a xy cross the straight line dictated by Galilean invariance, 



(41) 



ne e 2 e 2 p 
xv = — = v= , (42) 



B 2irfr 2tt11 2ps + 1 ' 

and 



LLL states flux quanta density eB ^ 1 _ l 
particles particle density 2-Khn p 



(43) 



Hereafter we set h = 1. 
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A. Gauge choices 

There are two famous choices for gauge. In the Landau gauge A = —}Bx the Hamiltonian is cyclic in y, and hence 
has the eigenfunctions e lky . Making the ansatz 

il>(x,y) = e lky <j>(x) (44) 
we find that cj> obeys the equation for a displaced harmonic oscillator 

" 5^ + 5""*« < 45 » 

The LLL corresponds to putting the oscillator in its ground state. The degeneracy of any LL can be computed by 
considering a sample with sides L x and L y , with periodic boundary conditions in the y direction. This forces k = ^p.. 

Since the wavefunction is centered around x = kl 2 < L x , we must demand < j < ■ Thus the degeneracy D of 

each LL, is (going back to a square sample) D = in agreement with our prior gauge-invariant counting of states. 
We shall use this gauge in the computation of the magnetoexciton. 
In the rest of this article we employ the symmetric gauge in which 









A = 




(46) 


f) = 


1 ,2- 

-r + rz x p 


(47) 


R = 


1 ,2- 

-r-^z x p. 


(48) 



B. Getting to know the LLL 



As with any simple harmonic oscillator, we can construct ladder operators from the canonically conjugate pair n x 
and rj y (with l 2 playing the role of K) 



4 = kM 



at = A^iVx- iVy) 



The LLL condition a n \LLL) = gives 
where z = x + iy. A basis for ip is 



ip = e 



f(z) 



m = 0, 1, . . . 



(49) 
(50) 

(51) 
(52) 



The Gaussian is often suppressed. The state has angular momentum L z = m. 

If v = 1, (one electron per LLL state) there is a unique noninteracting ground state, which may then be perturbed 
by standard means, 



Xi = ~ z j) ' Gaussian = Det 



i<j 



r o _1 r 2 

z i z i z i 

r o _1 r 2 

^2 z 2 '-2 



■ Gaussian 



(53) 



For v < 1, we often want to focus on the limit w ^ oo and work entirely within the LLL. If in H = T + V we set 
T equal to a constant (eB/2m per particle), all the action is in V. Why is this a problem if V is a function of just 
coordinates? After all, 



and 



V = 



lJ d2r J d 2 r'p{v)v{v - r>(r') 



P(q) 



J d 2 rp{v)e-^ = £\ 



(54) 

(55) 
(56) 

(57) 
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The point is that if one sets T — eB/2m, the LLL value, one must project the operator r to the LLL. The 
coordinates x and y, which commute in the full Hilbcrt space, no longer commute in the LLL and V now contains 
noncommuting operators. 

C. Projection to the LLL 

Let V denote projection to the LLL 18 . Then 

V : r = R + r7^R. (58) 

The projected components do not commute: 

[Rx , Ry] = -H 2 or [z , z] = -2/ 2 in the LLL. (59) 

As for the densities 

V : e- lq r (e^ 4q ^) LLL e- ,q R = e -9 2(2 /4 e -«j-R ( 60 ) 
Thus the projected problem is defined by 

q 

p(q) = £ e -^ R , (62) 



The projection p is the Magnetic Translation Operator 19 , which differs from the projected density by a factor e~ q 1 / 4 . 
We shall often refer to it as the density, but take care to include the gaussian in Ean. l|61|) . The commutation rules of 
p define the magnetic translation algebra: 



[p(q),p(q')] = 2 * sin 



(q x q') I 2 



p(q + q') (63) 



which was thoroughly exploited in the work of Girvin, MacDonald, and Platzman (1986). 

There is no small parameter in H and the overall energy scale is set by v(q). This is why the FQHE problem is 
unique. As mentioned earlier the Hartree-Fock solution is also not an option at fractional filling. 

D. Kohn's Theorem 

By Kohn's theorem we shall mean the following results. Consider, in a translationally invariant system, the 
density-density response function K(q,uj), the Fourier transform of 

K(d ,t) = i9(t)(0\ [p(q,t),p(-q,0)] |0>. (64) 

In the frequency domain we define 

/CO 
K{^t)e iut dt. (65) 
-OO 

^From just the canonical commutation rules it can be shown that K obeys the sum rule 

[°°T K( \ dbJ 

/ Im K(q,uj) — = — . (66) 
Jo 71" 2m 



A very nice introduction to LLL physics appears in Girvin and Jach (1984) 

These operators and the projective group they form have a long history, and were first used to describe symmetries of the noninteracting 
electron Hamiltonian in a magnetic field. For references, see: Peterson (1960), Brown (1964), Zak (1964a, b). To the best of our 
knowledge, Girvin, MacDonald, and Platzman (1986) were the first to concentrate on the algebra of this operator 
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Kohn showed that K must have a pole, the magnetoplasmon, at the cyclotron frequency loq = eB/m with a residue 
that saturates the above sum rule as q — > 0. 

It follows that p(q), the restriction of p(q) to the LLL, cannot have (transition) matrix elements that are linear in 
q for small q. 

The structure factor S(q, ui) 

S(q,uj) = ^|(0|/;( ?) 0)|n)| 2 5(a; - E n ) (67) 

is related to K (q, u>) for lu > by 

-7m K(q,u>) = S(q,u). (68) 
20 Kohn theorem tells us that if we limit ourselves to the LLL, S(q) ~ q 4 for small q. 



III. HAMILTONIAN THEORY I - THE CS APPROACH 

The hamiltonian for electrons in a vector potential A is 

g = £ (p i+: A (rj .))> +v m 



2m 

3 

where V is the electron-electron interaction, say Coulombic. Disorder is not included. As stated earlier, the field A is 
such that there are 2s + - flux quanta or LLL states per electron and the attendant degeneracy frustrates perturbative 
analysis. 

The first step in the CS approach is to deal with the degeneracy by resorting to flux attachment. For the Laughlin 
fractions there are actually two options involving either composite bosons (CB) or composite fermions. 



A. Composite bosons 

Historically, the first treatment came from Zhang, Hansson and Kivelson (1989) who considered Laughlin fractions 
v = I /{2s + 1). They traded electrons for CB's carrying 2s + 1 flux quanta in opposition to the applied field so that 
at mean-field level the bosons saw zero field and had a unique ground state 21 . The trading is done by introducing a 
Chern-Simons wave function ^cs defined as follows: 

= II T7 Z Zj) \2s + +1 = GX P(( 2S + ^ E hi) ( 70 ) 

\Zi Zj\ . 

„ \ - (p t + eA(rQ+a cs (r a )) 2 

H °s = ^ + V (71) 

i 

where <j>ij is the phase of the coordinate difference Zj — Zj. Since ^ e changes sign under particle exchange and the 
prefactor produces 2s + 1 extra minus signs, ^cs describes bosons. 

The CS gauge field, a cs , comes from the action of p on the prefactor (which is just the phase of the Jastrow factor) 
that multiplies ^cs'- 



a cs {r, t ) = 2sV^2<f>ij (72) 

a cs (r l ) • dYi = 2s <j)^2 ^4>ij ' drt (73) 

= 2tt(2s + 1) (number of particles enclosed) (74) 

V x a cs = 2vr(2s + l)p (75) 



To see this, introduce a complete set of exact eigenstatcs |n) of H between the two factors of p in Eqn. 1641 . express the Heisenberg 
operator p(t) in terms of p(0) and H ,a,nd do the time integral. 

This mean-field idea was first applied to anyon superconductivity by Laughlin (1988). Many more works on this topic followed: Fetter, 
Hanna, and Laughlin (1989), Hanna, Laughlin, Fetter (1989), Chen, Halperin, Wilczek, and Witten (1989), Halperin, March-Russell, 
and Wilczek (1989), Lee and Fisher (1989), Hanna, Laughlin, Fetter (1991), Dai, Levy, Fetter, Hanna, and Laughlin (1992). 
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Eq. I|75(l shows explicitly that the flux quanta density is 2s + 1 times the particle density. This is what one means 
by flux attachment. Eqs. (|71|l and (|75|> define a Chern-Simons theory. 22 The possibility that the FQHE would be 
described by a CS theory was presaged by Girvin (1987). 

Since the idea of flux attachment is to cancel the applied field on average, Zhang, Hansson, and Kivelson separate 
a cs and p into average and fluctuating parts: 

V x (a cs ) + Vx : a cs := 2vr(2s + l)n + 2tt(2s + 1) : p : (76) 

This gives 

„ \- (p + eA+ (a ca )+ : a cs 

tics = > 7T- 1" V 

— ' Am 



E 



2to 

(p+ : a cs :) 2 
2m 



V (77) 



upon using the fact that the flux due to eA precisely cancels that due to (a cs ). Bosons in zero field have no degeneracy 
problem (assuming they have some repulsive interactions) and allow one to describe much of the FQHE physics in 
the familiar language of superfluids. For example, it has been shown by computing response functions, that the 
superfluidity of the bosons implies the FQHE for electrons (Zhang, Hansson, and Kivelson , 1989), that the vortex 
in the superfluid is Laughlin's quasi-hole (Lee and Zhang, 1991) The nature of the collective modes has also been 
explored (Kane, Kivelson, Lee, and Zhang, 1991). 

Neglecting : a cs : (the mean-field approximation) and the interaction, ^cs the wavefunction for bosons (in Eqn. 
{7UJ|) is just unity and that of electrons is 

/ _ \ 2s+l 

•-n(it^) - 1 ™ 

which is the phase of Laughlin's answer. Kane, Kivelson, Lee, and Zhang (1991) showed that if long-wavelength 
gaussian fluctuations are included, the full Laughlin wavefunction is obtained at long distances. But the same 
fluctuations also reduce the long-range order in the boson field down to the power-law order found by Girvin and 
MacDonald (1987) who analyzed a gauge-transformed version of the Laughlin wavefunction. 23 Read (1989) then 
showed that for Laughlin fractions, one could form an operator (which was a composite of an electron and 2s + 1 
vortices) which was neutral and had true long-range order. The corresponding Landau-Ginzburg theory for the order 
parameter was, however, very complicated. Constraints of time and space prevent us from describing CB's any further. 
We refer the reader to primary sources and excellent reviews (Zhang, 1992, Karlhede, Kivelson, and Sondhi, 1993). 

There are, however, some shortcomings in the CB approach. First, it is restricted to Laughlin fractions. Next, since 
bosons have to be interacting to be stable, a noninteracting starting point does not exist. Finally, there is singular 
dependence on to as to — > and it is hard to carry out quantitative computations. Nonetheless, this work has served 
as a paradigm for the hamiltonian approach. 

B. Composite fermions 

We turn now to composite fermions. It was seen in the review of Jain's work that by trading electrons for CF's 
which carry 2s flux quanta in opposition to the applied field, one could get a fermionic system that on the average 
sees a field A* that is just right to fill p Landau levels. Lopez and Fradkin (1991,1992,1993,1998) were the first to 
accomplish this in the hamiltonian approach. 24 They traded the electronic wavefunction 'J'e for vE'cs defined as 
follows: 

* e = II u'Z^.p ^cs = exp(2wj^) * CS . (79) 



22 Note that these manipulations could just as well be done in the path integral formulation. We have chosen to use the first quantized 
operator version for all our discussions, in the interest of uniformity. 

23 This is to be expected given that the gauge transformation of Girvin and MacDonald (1987) is the same CS transformation of Eq. 1701 . 

24 Note that we do not distinguish between the functional integral formalism which they employed, and the operator approach which we 
used in our work, and in the interest of uniformity, throughout this paper. 
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H CS = £ (P* + eA ^+ a ^)) 2 + V . (80) 

i 

Since \I/ e describes fermions, so does ^cs since the phase factor is even under particle interchange. The CS gauge 
field a cs now obeys 

V x a cs = Ansp. (81) 
Separating a cs and p into average and fluctuating parts 

H cs = V (P + eA+ (a C8 )+ : a cs Q| | y 



2m 

E (n ± _^0i + v m 



2m 

n = p + eA + (a cs ) = p + eA* (83) 

V x (eA + (ac S )) = -eB + Aims (84) 

pR A 

-eB* (A* = -) (85) 



2ps +1 v 2ps + 1 

I* = = /v/2p.s + l (86) 

yeS* 

The following important results emerged from the work of Lopez and Fradkin (1991,1992,1993): 

• If we ignore : a cs : and V, the composite fermions see ^ flux quanta each (since 2s + ^ — ► ^ under flux 
attachment) and have a unique ground state \v °f P filled LL's. Excitations are given by pushing fermions into 
higher CF-LL's. 

There is however a problem: If we excite a fermion from level p to p + 1, the energy cost (activation gap) of the 
particle- hole pair is A = eB* /m plus corrections due to neglected terms. This divergent dependence on m flies 
in the face of the nonsingular m — > limit we have argued must exist 25 . We want A ~ e 2 / 'el in the Coulomb 
case. 

• At the mean- field level, the CF wavefunction \p, transformed back to electrons is 



i<3 



Lopez and Fradkin showed that fluctuations at one loop give the square of the wavefunction (at long distances) 
for Laughlin fractions. (The factors \zi — zj\ 2s in the denominator of the mean-field wave function are eliminated 
by fluctuations.) 

• They calculated time-dependent density-density response functions in the Random Phase Approximation (RPA) 
(which will be explained in the next section). The cyclotron mode appears with the right position and residue. 
However, between the cyclotron mode (eB/m) and the LLL excitations, there are many spurious modes at- 
tributable to the ubiquitous presence of m which prevents a clear separation of LLL and non-LLL energy scales. 
This is a problem common to all CS theories and was well appreciated by the authors. 

The Lopez-Fradkin work paved the way for subsequent work to which we now turn. 



IV. PHYSICS AT AND NEAR v = \ 

In the early days the CF, there was a widespread belief that its utility was confined to fractions with a robust 
gap: the all-forgiving gap allowed one to neglect, in a first approximation, interactions, disorder and gauge field 
fluctuations. It was therefore quite a surprise to see that the CF survived even when p — > oo or v — > l/2s, and the 
CF cyclotron gap eB* /m = eB/(2ps + l)m approached zero. 



Jain does not have this problem since he does not use Hqs or Xp or its excitations directly. For him the CS picture is a step towards 
getting electronic wavefunctions for the ground and excited states by attaching the Jastrow factor and projecting. The energy gap is 
computed as the difference in (V) between the ground and excited electronic wavefunctions. 
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A. Physics at v = i 

Kalmeyer and Zhang (1992) were the first to discuss the case v = 1/2. Their work emphasized the following 
important point. One may expect that the effect of disorder will be rather small since the electron donors lie at a 
respectable distance from the electron gas itself. However any small charge inhomogcneity induced by disorder is 
accompanied by a corresponding flux density in a CS theory, and this can cause significant scattering. 

A landmark study of the region at and near v = 1/2 was made by Halperin, Lee and Read (1993), (HLR) who 
seriously pursued the remarkable possibility that a Fermi liquid could be hiding deep in the FQH regime. HLR started 
with the following hamiltonian at v — 1/2 : 

^ = E (p+ ;^ :)l +y ( 88 ) 

i 

where : acs ■ is related to charge fluctuations by the CS condition 

Vx : a cs := 4tt : p : . (89) 

Note that CS theory dictates that to must be the bare mass, since flux attachment by minimal coupling is a gauge 
transformation performed on the bare electron. If one wants to take the view that this is an effective low-energy 
theory with an effective mass to*, one must also be prepared for the possibility that the coupling of the fcrmion to 
the gauge field is more complicated. 

The cornerstone of HLR's work is the computation of the electromagnetic response functions. Let us recall some 
general results on this topic so we may better understand and appreciate their work. In writing this description we 
were greatly aided by the reviews of Halperin (in Das Sarma and Pinczuk, 1997) and Simon (in Heinonen, 1998). 

If an external four-potential eA ext (\, u) is applied to the system, it will generate a four-current j M (q, to) (whose 
components are the number current and density) as per 

j M (q,w) - eK^(q,u)Al xt (q,Lo). (90) 

Linear response theory tells us that 

/oo 
dte^H8(t)(0\ [j M (q,t),.7„(-q,0)] |0> (91) 
-oo 

where |0) is the vacuum state. 

For pedagogical purposes consider just K a0} the density-density correlator. In free- field theory, K ao is given by 
Kq , the particle- hole bubble, in which the particle and hole created by one p are absorbed by the other. The full 
theory of course requires us to include interactions. Once again, for pedagogical purposes, let us begin with the 
case with just the Coulomb interactions. The exact K 00 is given by an infinite sum of Feynman graphs in which the 
electrostatic propagator v(q) — 2ne 2 /q appears in all possible ways. It is possible to organize the sum as in Figure (3). 
Each bubble K^ , called the Coulomb-irreducible response, has the property that it cannot be cut into two disjoint 
pieces by snipping just one Coulomb propagator v(q), denoted by the wiggly line connecting the irreducible bubbles. 
Performing the geometric sum one finds 

Km = floo = i . (92) 

The response K^ has the following significance. Consider some conducting system with an applied potential e<j) ext . 
An electron inside the conductor feels in addition the potential generated by the charges themselves, i.e., it feels the 
total potential e<f) T = ecf) ext — v(q)p(uj 7 q), and Kq Q is the response to this total field. To verify this let us write 

p(u, q) = eKSo [e4> ext - v(q)p(u>, q)] (93) 

solving which we get 

Who cares about the total field? Answer: The voltmeter we use to measure the drop in a wire. The voltmeter 
responds to the total electric field and not to the externally applied one that would exist in the absence of the 
conductor. The longitudinal conductivity a xx can be related to Kq as per 



V 



21 



FIG. 3 The diagrammatic relation between the full response function Koo and the irreducible part Kqo for a problem with 
just Coulomb interactions, denoted by the wiggly line V = v(q). In RPA K^o is approximated by the particle-hole bubble K®o- 



where we have used the continuity equation qj = up and E — iqcj). 26 

In the RPA Kq is approximated by Kq , the free particle-hole bubble. Thus RPA takes into account Coulomb 
interaction via the internal field the particle themselves generate, but ignores all vertex and self-energy corrections 
coming from exchanging v(q) inside the particle-hole bubble. Figure (3) shows the connection between the irreducible 
and total responses. 

The HLR work differs from this illustrative example in that there are now two types of gauge fields, the Coulombic 
v(q) and the CS gauge field, so that the wiggly line in Fig. (3) is described by a 2 x 2 matrix propagator. The HLR 
version of RPA consists of summing repeated bubbles irreducible with respect to both propagators, and approximating 
the irreducible part by the free-field response. Thus the CS field and v(q) are included only to take into account the 
internal fields the induced charges and currents produce. In matrix notation it is still true that 

[K]- 1 =[K°]~ 1 + U (96) 

where U is a matrix propagator for the Coulombic v(q) and the CS gauge fields. Let us focus on just the 00 element 
-Koo, which assumes the form 

^oo = 1 . _, 2 (97) 

v{q) + [iCF 1 + *?! 

where is the free-field transverse-current response, and </> is the number of flux quanta attached. Though 4> = 2s 
in our notation, we still use HLR's notation here to help readers who may wish to consult that work for more details. 
Let us now examine this expression in various regimes. 



1. Static compressibility 

The behavior of -Koo at u> = 0, q ^ 0, gives the static compressibility of the system. For u> = and q << hp, they 
find 

K °° = v(q)+2TT(l + ft/6)/m (98) 

which shows that just as in a Fermi liquid, the compressibility is nonzero if v(q) is short ranged and vanishes as q if 
Coulombic. This vanishing only means that any applied external force is unable to change the density because the 
field inside the medium is strongly screened by the medium. 



2. Cyclotron mode 

The poles of -Koo define the natural modes of oscillation of the system, for at these poles we get a response with 
no applied potential (see Eqn. I|90|) ). If at any q there is a pole at some ui(q), it means that there is a mode of energy 
to(q). For example, in a 3D electron gas, a plasmon pole appears at the plasma frequency as q — > and then moves as 
a function of q. At larger q, when decay into particle-hole pairs is possible, the pole position acquires an imaginary 
part denoting a finite lifetime. 

In the present case, as q — > 0, at very high w, -Koo has a pole at the cyclotron frequency ujq = Ann/m (where m is 
the bare or band mass) with a residue in accordance with Kohn's theorem. 



The extra e in front comes in because p and j refer to the particle density and not charge density. 
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3. The overdamped mode. 

At very low to and q, one finds 

m 12 V 3 / 9*f 

The zero occurs at the overdamped mode with the dispersion relation 

w ~ ig 3 w(g) ~ <? 2 (100) 

for the Coulomb case. If w(g) is short ranged, the mode is sub-diffusive. The reason charge diffuses rather than 
moving ballistically (even in the clean system) is that the magnetic field makes the charge move perpendicular to the 
Coulomb force, which tries to even out the density gradient. The overdamped mode and the effects it produces lie at 
the heart of the HLR work. 

Note that if one first sends m — > 0, one would miss the over damped mode. Of course in a real system we can 
always find a q such that l/q dominates 1/m. It is however a theoretical problem that we cannot send m — ► here. 
A way out will be discussed later. 



4. Longitudinal conductivity 



Upon examining Eqn. I|99l) . we can see that the Coulomb- irreducible part (obtained by dropping v) is dominated 
by the last term in the limit of small q and u>. Thus 



n 

e q 



The above result holds only for q^> l/l m , where l m is the mean free path for the CF. For q <C 1/1,- 

^2 



It is interesting to see how one arrives at Eqn. i|102|) . 
Let us define a CF conductivity by 



ej = o-cfE 



(101) 



(102) 



(103) 



where E T is the total field, which is the sum of the applied field and the internal field e generated by the fermions 
themselves. To compute e imagine a particle current j and a unit length perpendicular to it. In one second j particles 
cross it carrying with them 2j flux quanta of CS flux. However the CF makes no distinction between real flux and 
CS flux since they enter Hqs only via their sum. It will therefore sense a CS electric field e = — z x j. Thus 



ej = o C f 



4-7T 

E H z x j 

e 



(104) 



The electron resistivity tensor p (with components p a p and not to be confused with the density operator) defined by 



is then 



where 



E a = p a p (ejp) 



P = PCS + PCF 



PCS 



47T 
72" 



1 

-1 



At mean-field level CF's do not see a field of any kind. The effect of disorder is to produce 

TIG T TIG Ijyi k p6 Ifn ^ k ]? 



0~CF, 



0~CF. 



-in 



(105) 
(106) 

(107) 
(108) 
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where r is the elastic scattering time and l m is the mean free path. If we now use Eqn. I|106[) and assume fcpi OT 1 
we obtain Eqn. I|1U2I) . 

Note that the resistivity of the electron is the sum of the resistivities of the CF and the CS term. If we move off to 
a general Jain fraction and ignore disorder, we will have 



Pcf 



2n 
72" 



i 
1 p 
-A 

p 



(109) 



since the CF's fill p LL's of their own. When the resistivity matrices are added, one obtains the correct Hall resistivity 
for electrons at v = p/(2p+ 1). 



5. Surface acoustic waves (SAW) 

When a surface acoustic wave is coupled to the electronic system, it is predicted to undergo a velocity shift and an 
attenuation described by 

5v„ in a 2 12 

— = (110) 

v s q l + io- xx (q)/cr m 

where a is a piezoelectric constant, v s is the sound velocity, k describes the attenuation, o~ xx (q) — a xx (q,uj — qv s ) 
where v s is the sound velocity, and a m — Theory fits the experiments of Willett et al (1990,1996) (for a review 
see Willett (1997)) with a a m that is about five times larger than expected. The reader is strongly urged to consider 
HLR and the reviews for more details. 



6. Mass divergences 

HLR (1993) predicted a divergence in the effective mass m* at the Fermi surface arising from the fermion self-energy 
diagram involving the emission and absorption of the overdamped mode. For the Coulomb case one has 

m*(w) ~ baw. (Ill) 

(Shorter range interaction lead to more violent divergences.) Assuming this mass can be used near v = 1/2 one 
expects that the gaps will be given by 

E v = ^ (112) 
m 

where m* is self-consistently defined by m* (ui — E v ). For v = pj (2p+ 1), this implies, as p — > 00, that E p ~ l/(plnp). 
The log divergent effective mass will also imply a specific heat C(T) ~ TlnT. It has not been possible to confirm 
these logarithms in numerical or experimental work. 

Since the nature of the mass divergence depends on the range of the potential, it cannot be reproduced by trial 
wavefunctions, which make no explicit reference to the potential, in particular its range. 

It has been shown that the mass divergences do no affect bosonic (e.g., density-density) correlations 27 . 



B. Physics near v = 1/2 

The HLR (1993) predictions transcend v = 1/2 and describe its immediate vicinity. The key idea is that in this 
region, it is useful to think of the CF as a particle seeing a weak effective magnetic field B* — B/(2p+ 1). This means 
that it will describe a cyclotron orbit of radius 

R* = ^f (113) 
eB* v ; 

with kF = V47TO (given that spin is fully polarized). Note that this result is independent of the fermion mass, whose 
treatment is quite tricky. 



Kim, Furusaki, Wen, and Lee (1994), Kim, Lee, Wen, and Stamp (1994), Kim, Lee, and Wen (1995), Stern and Halperin (1995). 
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The SAW results of Willett et al (1990,1996) found that away from v = 1/2, there was a resonance in the velocity 
shift when the wavelength of the surface acoustic wave coincided with 2R* . 

Confirmation of the R* concept was also found in the experiments of Kang et al (1993), Goldman, Su, and Jain 
(1994), and Smet et al (1996). To visualize the Goldman et al experiment, imagine a semi-infinite system in the upper 
- half plane. If a current is introduced at the origin up the y-axis, it should bend and return to the x-axis at x = 2R* , 
after completing one semi-circle. It would then bounce off and start the next semi-circle. It follows that if a return 
path is provided on the x-axis, the maximum current will flow if the drain is located an integral multiple of 2R* from 
the source. Or if a return location is held fixed, and B is varied, a maximum is expected and (found) whenever the 
spacing between source and drain is a multiple of 2R* . 

Kang et al built an array of antidots, where each antidot is a region where electrons are absent. It was found 
that the sheet conductance had peaks when the antidot lattice constant equaled 2R* . While a detailed formula for 
conductance is not known in this context, the correlation is very suggestive and is what was seen with ordinary 
electrons in a weak field. The fact that CF's, which entered the theory as a mathematical devise, manifest themselves 
so clearly in transport is a stunning affirmation of the theoretical framework. 

C. HLR - Room for improvement? 

Despite its many remarkable and experimentally confirmed predictions, the HLR theory leaves room for improve- 
ment mainly because of the CS approach on which it is based. 

One such area has to do with the dependence on bare mass to. We have argued that at the starting point Hcs 
must contain the bare mass since flux attachment by minimal coupling the CS field is an exact transformation done 
on electrons. This choice certainly helps to get the cyclotron frequency in accord with Kohn's theorem. However in 
our extraction of various low-energy long wavelength quantities like the over-damped mode from Eqn. 199fl . we had 
to assume that 1/q dominates over 1/m. While this is certainly valid for realistic and fixed values of to, one would 
like to be able to extract what is evidently the correct physics even if to = is imposed first. 

Simon and Halperin (1993) and Simon, Stern, and Halperin (1996), proposed a way out. They assume, as in 
Landau theory, that Hcs is an effective theory with an effective m*. Kohn's theorem and the over damped mode can 
both be salvaged if a suitable Landau parameter F\ is introduced. This is reasonable since Kohn's theorem relies on 
Galilean invariance and Landau's theory uses this principle to related m, m* and F\. However this leaves the origin of 
to* inside a black box. 28 It is also not clear why in the effective theory the fermion the CS field should be minimally 
coupled. 

In the case of Coulomb interactions, one can also argue that the bare mass is swamped by the renormalized 
to* ~ hxu) generated by the exchange of the over-damped mode. Again it would be nice to be able to follow in 
detail the separation of the low energy physics controlled by an m* generated by interactions and high-energy physics 
controlled by to. 

Another shortcoming is that there is no evidence of the neutral fermion one expects at v = 1/2 , and more generally 
a fermion of charge e* = e/(2ps + 1) at other Jain fractions. This was to be expected since the CF in Halperin et al 
and in earlier CS work of Lopez and Fradkin was an electron bound to two flux tubes which carry no charge. 

Also missing was the effective magnetic moment /i* = e/2m that Simon, Stern and Halperin (1996) would later 
argue must go with each CF. 

In short, the CF that appears in the CS theory does not have, in obvious form, the right charge or energy scale 
of the ultimate quasiparticle. While the CS transformation is exact and can yield all these in principle, they are not 
manifest in the RPA. 

Finally, Lee, Krotov, Gan, and Kivelson (1997, 1998) have raised the following issue. Suppose we push all LL's 
up to infinity except the LLL and assume particle-hole symmetric disorder at v = i Then it can be verified that 
a nonzero a xx for electrons implies for CS-fermion a a X y = — e 2 /(47r). However, both the mean-field approximation 
and the RPA corrections to it, give cr^y — 0- We are not aware of a resolution of this issue, arising from the fact that 
resistivities add. 



Why do we view this inability to calculate m* as a weakness, when in Landau theory it is simply accepted as a fact of life? As we shall 
see, in the FQHE we can go a long way towards computing m* from the interactions by using the EHT. 
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V. HAMILTONIAN THEORY II- EXTENDED HAMILTONIAN THEORY (EHT) 

We now turn to our extension of the CS formalism 29 . In this extension, a CF with all the known properties can 
be made manifest and a variety of LLL quantities computed with no 1/m singularities. But the computation of 
certain low energy long wavelength quantities (like the compressibility), straightforward in the CS approach, become 
extremely delicate. Hopefully we will impart to the reader a sense of which approach to use for what purpose. 

We want to eschew the historical route and furnish an axiomatic presentation of our work with the minimum of 
preamble, setting the stage for exact or approximate calculations as rapidly as possible. However, to establish the 
context, we will begin with a brief sketch of our earlier work so as to take some of the mystery out of the end product 
and provide some degree of motivation. 

Let us recall the work of Bohm and Pines (Bohm and Pines, 1953) on the electron gas in three dimensions. At 
small q the spectrum consists of particle-hole pairs at low energies and the plasmon at high energies. The particle 
and hole are part of the original hamiltonian while the plasmon comes from summing an infinite class of diagrams 
in the density-density response. It is not an independent entity, even though it is a sharply defined excitation that 
can be experimentally produced and detected. Bohm and Pines showed that by introducing extra canonical oscillator 
degrees of freedom at small q, one could describe plasmons as independent objects in an enlarged Hilbert space. In 
order to prevent double counting, they imposed constraints on state vectors of the form 

x(q)|physical state) = (114) 

with one x(q) for each q at which a plasma mode was introduced as an independent canonical oscillator. The plasmons, 
initially coupled to the fermions, were approximately decoupled, leaving behind fermions with a renormalized mass, 
and constraints that essentially froze out any putative plasmons with small q. Such a simple description of plasmons, 
or isolation of high-energy physics, would have been impossible within the confines of the electronic Hilbert space, 
wherein plasmons are complicated collective excitations of the electrons themselves. 

What we originally did for the FQHE (Shankar and Murthy, 1997, Murthy and Shankar, 1998a) was similar in 
many ways. We enlarged the Hilbert space to include at each q, a new set of independent canonically conjugate 
variables - a transverse vector field a(q) and a longitudinal vector field P(q). Using these it was possible, by a unitary 
transformation, to get rid of the dependent CS vector potential acs- While all the operators in our hamiltonian were 
now independent, the physical sector was defined by the (CS) condition 

(V x a - 47rsp) | physical state) = 0. (115) 

The conjugate variables a and P formed oscillators near, but not exactly at, the cyclotron frequency. They were 
coupled to fermions. We found a way to decouple the oscillators in the limit ql — > by a second unitary transformation. 
The decoupled oscillators now complied with Kohn's theorem for pole position and residue. The formula for p(q), the 
electronic charge density projected to the LLL, was derived to order ql, as was the constraint x(q) which now acted 
only within the particle sector. At this point a guess was made (Shankar 1999) to extend the small ql answer to all 
orders by exponentiation of leading-order terms. The final theory of the LLL sector took the following form: 

H = ^^(qM^e-^/^-q) (116) 
q 

H = (117) 

= x(q)|physical state ) (118) 

The hamiltonian H is just the potential energy V projected to the LLL, written in terms of the projected electron 
charge density p(q) 

p(q) = ]TcxpHq-R e ,) (119) 

3 

I 2 

R e = r z x n (120) 

1 + c v ; 

where we define a very important and frequently occurring variable: 

c 2 = (121) 
2ps + 1 V ' 



Murthy and Shankar (2002), Shankar and Murthy (1997), Murthy and Shankar (1998a,b, 1999, 2002), Murthy (1999, 2000a,c, 2001b,c), 
Shankar (1999, 2000, 2001), Murthy, Park, Shankar, and Jain (1998). 
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From the commutation relations 



(122) 



it is clear that R e is just the electron guiding center coordinate but now expressed in terms of CF variables r and 
n = p + eA*. 

The physical states are annihilated by the constraint 



x(q) 



^ exp(-iq • H vj ) 

3 

I 2 



c(l + c) 



z x n 



(123) 
(124) 



The pseudovortex coordinate R„ describes and object of charge 



-2ps/(2ps + 1): 



and commutes with R e : 

[Re } R«] = 0. 

Thus the constraints commute with H , and form an algebra: 

'i 2 (qxq') 



R(q),x(q')] 



-2i sin 



2c 2 



x(q + q')- 



(125) 



(126) 



(127) 



The problem is just like Yang-Mills theory. 

The expressions for R e and R„ in terms of r and II which were already encoded in the small-q theory, jumped out 
upon exponentiation. They, together with the constraints, lie at the heart of our approach. 

We will now show how one can get to the final result, Equations $116l - \TT8\) . if we simply make a certain enlargement 
of the electronic Hilbert space and follow it with a change of variables. Any reader who wants to know more about 
what makes us introduce these variables should consult our earlier work. Those who just want to use the results can 
go ahead. 



A. The axiomatic introduction to the EHT 



Let us begin afresh with the primordial hamiltonian in terms of electronic variables (which carry the subscript e 
to distinguish them from other coordinates to be introduced): 



H = 



^ 2ml 4 
3 



i,3,1 



v(q)e 



iq-(r 



= H + V. 



(128) 



This hamiltonian contains complete information about the problem. It can be used to study LL-mixing and the 
computation of Hall current which requires higher LL's in an essential way. These topics will be discussed in due 
course. We first focus on the main challenge of hamiltonian theories: extracting the m-independent physics of the 
LLL. As discussed in Section (fTTf) . projecting to the LLL as such is no problem: one drops the first term and makes 
the replacement r — * R in the density operator and v(q) — ► v(q)e~~ q 1 I 2 . The catch is that the projected hamiltonian 
lives in the highly degenerate LLL, frustrating both perturbation theory and the Hartree-Fock approximation. 

In the CS approach one resorts to flux attachment to beat the degeneracy of the kinetic term, but that can be done 
only in the full electronic Hilbert space. Consequently m gets into everything and the low (LLL) and high-energy 
sectors get hopelessly entangled. What we really want to do is work within the LLL and attach flux tubes, which in 
the LLL translates into vortices, by analyticity. However when we say vortices, we do not mean zeros of the wave 
function, for such a thing does not exist as a degree of freedom within the hamiltonian and as we have seen, there 
are not enough of them in the wavefunction to go around anyway. Instead we mean by a vortex some object which 
has the charge of the 2s-fold vortex and corresponds to an excitation that can be created by inserting 2s flux quanta 
into the Hall system. Since such an object does not exist in the original Hilbert space (as an independent entity) we 
enlarge it to make room for this entity, which we call the pseudovortex, the vortex to emphasize its similarity to the 
vortices in the wavefunctions, and the pseudo to emphasize its differences. 
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The enlargement of Hilbert space can be explained in terms of just one electron, with cyclotron and guiding center 
coordinates rj e and R e . Let us temporarily focus on just the LLL physics and ignore r] e which does not participate 
in the change of variables, and will be reinstated subsequently. 

First we introduce an extra guiding center coordinate R„ (the pseudovortex), defining it algebraically by its 
commutation relations which represent a charge — c 2 

if 

[R vx ,R vy ] = —. (129) 
c z 

Next we combine R e and R„ to form the CF space. (Note that it takes two canonical pairs to make one regular 
fermion in two dimensions.) This CF space can be defined in terms of either the position r and velocity II of the CF, 
or in terms of its cyclotron and guiding center coordinates, r\ and R. Note that henceforth variables carrying 
not indentifying subscript will refer to the CF. The CF coordinates r\ and R obey the commutation rules of 
the cyclotron and guiding center coordinates of an object of charge e* = 1 — c 2 = l/(2ps +1): 

[Th,Vy] = it* 2 — i~~c? (130) 

[R X1 R V ] = ~il* 2 . (131) 
The rule for forming the CF from R e and R t , is the following: 



R = Re 1 ° 2 f v (132) 
1 — c z 

V = — — j(Rt, — R e ) (133) 
1 — 

The first equation says that the CF guiding center is the weighted sum of its parts. The second can be found by 
demanding that r\ be linear in R e and R„, commute with R, and have an overall scale that produces the right 
commutator. 

The inverse transformation is 

R e = R + r/c (134) 
R„ = K + rj/c (135) 

In terms of r and II, the CF coordinate and velocity operators 

I 2 

R e = r-- -z x n, (136) 

(1 + c) y J 

I 2 

R„ = r + z x n (137) 

c(l + c) 

which are just the expressions encountered in the brief historical review. 
Ignoring the zero point energy, here is where we stand in the LLL sector: 

H = i^>(g) e-^ 2 /2 e iq-(R.i-R«) (138) 

= I ]T v(q) e-^l 2 exp (iq • [(R, -B d ) + c{ Vi - Vj )] ) (139) 

*,i,q 

While it is true that we have managed to get rid of m and isolate the LLL cleanly, the reader may ask what we 
have gained, since algebraically the problem is the same as in electronic coordinates. The answer is that now there is 
a natural nondegenerate HF ground state in the extended space. This is because the HF hamiltonian is now written 
in terms of CF operators R and r) and the particle density is just right to fill exactly p- filled CF-LL's, i.e., \p 1S the 
ground state. The proof of its HF nature is found in Appendix[n] This key step opens up all the usual approximation 
schemes. 

Depending on what we want to compute, there are two distinct schemes. Both rely on the nondegenerate HF ground 
state, and both acknowledge a huge symmetry group of H , which comes from the following fact: H , as embedded in 
the CF space, does not depend on R„ = R + r)/c, the pseudovortex coordinate. Equivalently, it depends on R and 
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rj only through the combination R e = R + cq. The other (commuting) combination R„ = R + rj/c can be used to 
define a pseudovortex density 



x(q) = E e ~ iq ' R " ( 14 °) 

3 



which commutes with H: 

H ,x(q) 



(141) 



Let us understand this symmetry. If we view H as a function of 77 and R, or r and II, its eigenfunctions will 
depend on two coordinates, which can be chosen to be one component each of rj and R or just x and y. If we view 
H as a function of R e and R„, it depends only on R e . The energy eigenfunctions will be of the form 

ip(z e , z v ) = ip e (z e )ip v (z v ) (142) 

where ip v (z v ) is arbitrary since nothing in H determines it. This degeneracy is of the same type as that of the nonin- 
teracting electron hamiltonian which depends on 77 e but not R e . However, since R v is unphysical, and all physical 
observables depend only on R e , this is a gauge symmetry. 

We deal with the gauge symmetry as usual, by selecting a representative from each orbit, which in this case means 
eliminating the degeneracy due to the arbitrary function tp(z v ) that tags along for the ride. Let us make a particular 
choice ip v o(z v ). All we require is that it be translationally invariant so that (x(q)) = in this state. The gauge-fixed 
Hilbert space now consists of functions of the form ip e (z e )ip v o(z v ). In this sector x(q) — m the weak sense: any 
Green's function involving a string of x(q)'s will vanish since (i) x(q, i) = x(q, 0) = X(q),(ii) (x(q)) = 0> 30 m the 
one-dimensional space spanned by ip v o(z v ). (Imagine inserting the projector to this state between any two x(q)'s m 
the Green';s function.) 

In the path integral language we can make x(q) vanish weakly, that is to say, vanish whenever it appears in a 
Green's function, by a method similar to what is done in gauge theories. Imagine writing a path integral in the full 
CF space and then inserting a delta function imposing x(q) = for all q at any one time. 31 Since x(q) does not 
change with time, this restriction will hold automatically for all times and specify the fate of x(q) the way a gauge 
fixing terms does for the longitudinal degrees of freedom. 

The theory is thus defined (in schematic form) by the equations 

H = e-f^mK-ti (143) 

q 

jj , X (q)] = (144) 
X(q) - (145) 

where ~ means vanish weakly. 

Now we turn to the two approximate ways of dealing with this problem: the conserving approximation and the 
shortcut. 

B. The conserving approximation 

Given a hamiltonian, a reasonable first approximation to try is Hartree-Fock. It is shown in Appendix [C] that 
CF-LL states are Hartree-Fock states of H. How good is this Hartree-Fock approximation likely to be? For example, 
how will we fare if we compute the activation or transport gap in a fully polarized sample as per: 

A = (p + PH\H\p + PH) - (p|if |p) (146) 

where |p) stands for the Hartree-Fock ground state with p-filled LL's and PH stands for a widely separated particle- 
hole pair? 

There are at least two good reasons to expect that this naive HF result will require fairly strong corrections. First, 
if we compute the matrix element of the projected electron density between any two HF states, the answer will be 



30 An exception occurs if the string contains only x(0) = N. This does not affect what we plan to do. 

31 Since x(*j) involves r and p, this will have to be done in the phase space path integral. 
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linear in q, whereas in the exact theory, and within the LLL, it must go as q 2 as per Kohn's theorem. To see this 
note that 

e- iq - R * = l-iq-(R + CT7)+0(g 2 ). (147) 

While R has no transition matrix elements between different CF Landau levels, tj does. The second problem is that 
as ql — > 0, the projected electronic density (which reduces to . exp(— iq ■ Tj)) has unit contribution from each CF 
while we would like it to be e*. Evidently the HF result will receive strong corrections that will renormalize these 
quantities till they are in line with these expectations. 

These shortcomings are to be expected since the HF solution does not obey the constraint, or equivalently, does 
not factorize into the form ip = ip e (z e )ip v (z v ). 

The conserving approximation (Anderson, 1958a, 1958b, Rickayzen, 1958, Nambu, 1960, Baym and Kadanoff, 1961) 
is a sophisticated procedure for improving the HF state with additional diagrammatic corrections so that x(q) — 
in Green's functions. For v — 1/2 Read (1998) showed that this procedure restores Kohn's theorem, and reveals a 
dipolar structure for density-density correlations. We shall say more about this in connection with the compressibility 
paradox. 

Now we discuss the other approximation, the shortcut. 



C. The shortcut: The preferred charge and hamiltonian 

Consider the exact solution to the gauge fixed problem. Suppose, in the hamiltonian and elsewhere, we replace 
p(q) by the preferred combination 

p p (q)=p-c 2 x- (148) 

This makes no difference (to the computation of anything physical) in an exact calculation, since x is essentially zero. 

However, in the HF approximation it makes a big difference if we start with the hamiltonian written in terms of 
p p (q). To see why, consider its expansion in powers of ql: 

P P = J2 e " qrj (^TT~ ^ 2q >< n , + o • (q x n,) 2 + ■) . (149) 

• The transition matrix elements are of order q 2 between HF states because coefficient of q is proportional to the 
CF guiding center coordinate r — l* 2 z x II with no admixture of the CF cyclotron coordinate. This is more 
transparent if we use R and r] to write 

^ p (q) = (l-iq-(R + cT7) + ..)-c 2 (l-iq-(R + 77/c+...) (150) 
= (l-c 2 )(l-q-R + (q 2 ). (151) 

• With no further fixing, we see that the electronic charge density associated with ff(q) is now 1 — c 2 = e*. 

• We see from Eqn. (|149fl that when v = 1/2, the preferred density couples to an external electric field like a 
dipole of size I z x p. 

Thus there is not apriori need for strong corrections, at least in the long- wavelength limit. 
The hamiltonian we work with 

h p = ^p p (i)"(9)^ ,!,2/2 p p (-q) (1 52 ) 

q 

subsumes a lot of the right low-energy physics and CF properties. Unlike H which commutes with x(q), H P is weakly 
gauge invariant, that is 

[H(p>) ,x(q)]^0 (153) 

where the ~ symbol means that it vanishes in the subspace obeying x(q) = 0. Thus neither H{p p ) nor p p will mix 
physical and unphysical states. 

The significance of H(p p ) is the following. If the constraint x = is imposed exactly, there are many equivalent 
hamiltonians depending on how \ is insinuated into it. However, in the HF approximation, these are not equivalent and 
H(p p ) best approximates, between HF states and at long wavelengths, the true hamiltonian between true eigenstates. 
In contrast to a variational calculation where one searches among states for an optimal one, here the HF states are 
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the same for a class of hamiltonians (where x is introduced into H in any rotationally invariant form), and we seek 
the best hamiltonian: H(p p ) encodes the fact that every electron is accompanied by a correlation hole of some sort 
which leads to a certain e* , d* and obeys Kohn's theorem (q 2 matrix element for the LLL projected charge density). 

Note that when we use the preferred charge and hamiltonian we will make no further reference to constraints, 
and simply carry out the Hartree-Fock approximation. This is based on the expectation that even if we found some 
way to include the effect of constraints, it will make no difference in the small ql region. This is because the leading 
renormalization of e to e* and suppression of q matrix elements down to q 2 that are achieved by the conserving 
approximation by summing ladder diagrams following Read (1998) or using the time-dependent HF (Murthy, 2001a, 
Section VI) are built in here. 

It is in this approximate, operator sense, where we use p p in place of p that the binding of electrons and pseu- 
dovortices to form CF's is realized in the Hamiltonian theory. Since the pseudovortices have coordinates that are 
independent of the electrons, there is no double-counting here. The problems that we encountered with vortices in 
the wavefunction approach for non-Laughlin fractions (with the limiting case being the dipole picture of v = i) are 
also absent since we are not talking about those vortices. Since the pseudovortices per electron are independent of 
electrons, their number does not change when v changes (though their charge, tied to c 2 = 2sv, does.) Whereas 
antisymmetrization fragmented the vortices in the Jastrow factor into ordinary parametric zeros (except for the Pauli 
zero), it does nothing to the pseudovortices. Antisymmetrization is accomplished by simply writing the operator p(q) 
in terms of second quantized (composite) fermion operators. 

The reader will recall that any simple picture of quasiparticles, whether it be in Landau's Fermi liquid theory, or 

in BCS theory, is best captured by approximate and not exact descriptions. The quasiparticles are all caricatures 

of some exact reality and therein lies their utility. Similarly the CF in our extended formalism appears only in the 
— v 

HF approximation to H . Recall that we brought in the coordinate R„ to become the electron's partner in forming 

the CF. However R„ was cyclic in the exact hamiltonian H. Thus the exact dynamics never demanded that R„ be 

bound to R e or even be anywhere near R e . However, in the HF approximation, since we wanted the right charge 

and transition matrix elements of the density operator (Kohn's theorem) to be manifest, we needed to replace p(q) 

by p p (q), and trade H for H , the preferred hamiltonian. In H , R v is coupled to R e . The HF approximation and 

this coupling go hand in hand. The exact eigenfunctions of the original H are factorized in the analytic coordinates 

z e and z v and presumably reproduce the electronic correlations of the FQHE states. On the other hand, in the HF 

approximation to H , the wavefunctions (e.g., p- filled LL's) mix up z e and z v , and H , the preferred hamiltonian, 

dynamically couples R e and H v . The net result is that, at least at long wavelengths, these two wrongs make it right 

and mimic what happens in the exact solution. 

— p 

Another advantage of H is that it gives an approximate formula for to* originating entirely from interactions. 
This is best seen at v — 1/2. When we square p p (Eqn. I|149(l . we get a double sum over particles whose diagonal 
part is the one particle (free-field) term: 



This is not a hamiltonian of the form k 2 /2m*. However if the potential is peaked at very small q, we can expand 
the sine and read off an approximate 1 /to* 



which has its origin in electron-electron interactions. However we can do more: we have the full Hq as well as the 
interactions. The point to emphasize is that H is not of the traditional form (p 2 /2m + V ) and that there is no reason 
it had to be. 

The reader should verify that if we use H (p) instead, the one-particle piece will be a constant with no momentum 
dependence. The entire energy will be due to the Fock term, as in Read's (1998) conserving calculation. 

For the benefit of readers who may be overwhelmed by seeing too many approaches, we present the key equations 
of the CS and EHT approaches in Table [I] 

To summarize, in any context where p p can be reliably employed, we can say that the CF appears to be the bound 
state of an electron and the pseudovortex. We shall see that this includes calculations of the gap and magnetization 
properties at zero and nonzero temperatures, even at v — 1/2. We shall turn to these later on. But first, we list some 
cases where p(q) p and H fail to capture the right physics. 




(154) 




(155) 
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TABLE I Basic equations 

CS Theory 

H cs= ^EiCP + eA + ac,)? + V 
V x a cs = 4nsp 
Extended hamiltonian Theory (LLL only) 

#=§£ q «(?) e -^ i2/2 ^(q)p(-q) 
[ff ,X(q)] =0 

x(q) - o 
p(q) = E, ex p(-*q- R ej) 

R e = r — ji-^z x II = R + r)c 
X(q) = E, exp(-iq- R„j) 
R v = r + ^rqrjjz x II = R + rj/c 

# P = #(p p ) 
^(q) = ?(q)-c 2 x(q) 



D. The conserving approximation and compressibility paradox 

HLR predicted that the v = 1/2 system has a static compressibility that is nonzero for short range forces and 
vanishes as q for the Coulomb interaction. This result appears paradoxical in the present approach in which the 
dipolar nature of the CF has been transparently exposed (using ]f ). Imagine coupling the system to an external 
potential The dipole will couple to the gradient of $ and the resulting response will be a dipolar density whose 
divergence will give the induced charge. These two italicized factors imply a q 2 in the response even for short range 
interactions. Indeed, this is what we first obtained (Shankar and Murthy, 1997, Murthy and Shankar, 1998a) upon 
doing a simple RPA calculation. How is this to be reconciled with the HLR (1993) result, assuming their (compressible) 
answer is right? 

Halperin and Stern (1998) first raised this question and provided the key to its resolution. They first established 
a matter of principle, namely, that dipolar objects could be compressible, by considering the following hamiltonian: 

2 n,n,Q 



This hamiltonian arose in our earlier work (Murthy and Shankar, 1998a) when we decoupled the magnetoplasmon 
oscillators from the fermions. There we had chosen the upper cut-off Q — kp so that the i — j term from the 
double sum contributes — p 2 /2m to each particle and cancels the first sum, rendering I/to* = 0. Halperin and Stern 
considered the limiting case Q — *■ 0. In this limit H takes the form 

H _ (Pi - Pj) 2 _ ^ (p* ~ i gj Pj) 2 (l57) 

^ 2m ^ 2m 

and is invariant under the simultaneous shift of all momenta. This symmetry called K-invariance had also been pointed 
out by Haldanc (1995) in unpublished work and arose as part of a gauge symmetry in our work. The symmetry implies 
that it costs no energy to move the Fermi surface as a whole. Consequently there are some very soft modes that could 
lead to a singular density response which can offset the q 2 from the dipolar factors, provided these soft modes are not 
merely gauge artifacts that couple to nothing physical. The detailed calculation of Stern, Halperin, von Oppen and 
Simon (1999) (hereafter Stern et at) which we now describe, demonstrated that gauge invariant soft modes do exist 
and lead to nonzero compressibility. 

The first order of business for them was to start with a hamiltonian that had i^-invariance for small Q and not just 
Q = 0, since one needed to consider the response functions at small but nonzero q. The hamiltonian in Eqn. I|157|) 
had to be augmented by more terms to assure this. These terms were derived by Stern et al (1999) as follows. They 
go a step back in our calculation and start with the following hamiltonian and constraints for the coupled oscillators 
and particles: 

H = Y, ij ^ )1 +y ass) 



j 



y Zl_y-g( q) ^^ a( _ q) 

^ 2m ^ Sl4i 2mii CS V HJ 
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d 2 r [n cs (r)|a(r)| 2 ] +V 



g(q) = XlV ""• 

3 

= (V x a — 4irncs) I physical state) 



(159) 
(160) 
(161) 



where ncs is the fermion number density operator. The only minor difference from our work is that a x and a y (instead 
of the longitudinal and transverse components a and P) are canonically conjugate. The cut-off on q at Q means the 
Chern-Simons fermions now carry flux tubes, ones smeared over a distance 1/Q (sometimes called "fat" flux tubes, 
see Halperin, 1992). 

Next Stern at al (1999) approximate ncs by n, the average density and obtain 



H _ v 3 

' 2m 2mn 



n 

2m 



a(q) 



5(q) 



(162) 



Now they invoke the unitary transformation we employed to decouple the a fields from the fermions in the small 
ql limit: 



U — exp 



i 

47m 



Q 



J^g(q) x a(-q) 



(163) 



Since a x and a y are conjugates, U is just the shift operator a — g/n — > a'. Here comes the big difference. While we 
kept just the H from Eqn. (|157fl they augment H with additional terms to ensure that the constraints 



p(q) = — ^-q x g(q) 

commute with H aug to the desired order. The reader should verify that this constraint is just 



(q) = E e " iq ' R " = E ex p(" if i • ( r i + x = 



(164) 



(165) 



expanded to order To this order the electron density is (dropping a terms which do not matter at low energies) 



?e (q) = E e iq Rej = E cx p(-*q • ( r 3 -^ ix Pi)) 



E cxp Hq-rj)(l- yqxpj) 



(166) 
(167) 



p(q) - y q x g(q)- 



(168) 



Recall the constraint generates gauge transformations. To zeroth order the constraint operator is just J^j e lq rj = 
and its action is 



qe 



-iqrj 



(169) 



Thus respecting the constraint ensures if-invariance. 

Stern et al (1999) now perform an RPA calculation of the electronic density-density correlation. RPA works because 
at a given q, the gauge field (of that q) enters only in the wiggles connecting irreducible bubbles in Figure (3) and 
nowhere inside these bubbles: every internal exchange brings with it a sum over q which introduces a small parameter 
Q. Since in the limit Q — ->■ 0, RPA diagrams are all we have, RPA respects the symmetry of H aug . We will now argue 
that if the constraint is respected, compressibility follows. For the electron density operator we have many choices 
starting with Eqn. I|168|) and using the constraint Eqn. 1164f) . In particular we can write it as 



p e = 2p 



(170) 



Since the (CS) constraint commutes with H prior to the action of U , it does so after the action of U, but of course only to leading order 
in ql since U was not implemented exactly. 
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or as 

P e = -2yq x g(q). (171) 

Written the first way, a nonzero compressibility is not surprising since there are no powers of q in the operator. In 
the second way, we see the possible paradox, since there is a q up front in each of the two factors of p e . It is here that 
the overdamped mode appears in the transverse sector, couples to g(q) and saves the day with a factor q 2 v{q) in its 
static propagator, yielding results that coincide with HLR. 

The details of this formidable calculation are not shown here in the interest of brevity. Suffice it to say that five 
different operators (p, and the two components each of g and C, where C is another vector operator) get coupled 
and K^ 1 = Kq + U is a relation among 5x5 matrices. It is remarkable that in the end all the physics of HLR, 
including the overdamped mode, finite compressibility and mass divergences at the Fermi surface, all emerge from 
the more physical quasiparticlcs, though after a lot of work. 

However this is not quite the end since we have only discussed the Q — > limit, while the actual theory has no such 
limit on Q. Stern at al (1999) finesse this question with a two-step argument: (i) First they show that X-invariance 
guarantees that the Landau parameter Fx has to be —1. (ii) Then they show that this condition generally produces 
very soft modes which restore compressibility to the dipole gas. Let us elaborate a little. Suppose we want define 
an effective low-energy theory for the composite fermions. How do we ensure that it has X-invariance? We cannot 
possibly find all the higher order correction in the Q expansion. This is also the typical situation in Landau theory, 
that although in principle the Landau parameters can be calculated given the microscopic interaction, there is no way 
to do this reliably in practice. However, there are exceptions where a symmetry is involved. For example, Galilean 
invariance can be used to relate the bare mass, the physical mass, and the Landau parameter F\. It turns out that 
here too, K -invariance can be assured for the actual problem if we choose F\ = —1. The reason is that the energy 
cost of boosting the Fermi sea is measured by (1 + F\). Choosing F\ = —1, makes the boost cost-free, i.e., implies 
if-invariance. In other words the dipolar fermion cannot have an arbitrary Fermi liquid interaction: its F\ must equal 
— 1. At the level of diagrams, if F\ = — 1 is included as an interaction, the correlation function of two dipolar densities 
will go as q° and not q 2 . 

The compressibility of the v = 1/2 system was also established by Read (1998) using a conserving approximation. 
He begins with the problem of bosons at unit filling, a problem first studied by Pasquier and Haldane (1998), though 
not in a conserving approximation. The bosons can be traded for fermions in zero average field by attaching a flux 
quantum. The role of fluctuations is the same as in the HLR problem except for the strength of the gauge field-fermion 
coupling. The hamiltonian is just the electronic interaction written in the CF basis with 

I 2 

R e = r -zxp (172) 
and the constraint is the density corresponding to 

I 2 

R,=r|-zxp. (173) 

The advantage of this starting point is that the constraint exactly commutes with H. Thus one can look for 
approximation schemes in which the constraint is respected at the level of Green's functions, i.e., Green's functions 
with any number of x(q)'s in them vanish. Read's (1998) calculation begins with the filled Fermi sea (which does 
not respect the constraint) and embellishes it with diagrammatic corrections. An infinite sum of ladders (whose legs 
contain fermions in zero field) leads to the collective mode. The density-density correlator has the appearance of 
dipolar objects that exchange a transverse gauge field. The transverse propagator (which is just the overdamped 
mode) introduces a v{q)q 2 in the denominator that offsets the q 2 upstairs. It is found that correlators containing 
x(q)'s vanish and Ward identities are satisfied. 

It is instructive to look at Read's (1998) density-density correlator which takes the form: 



K irr (q,u) = j d 2 fcd 2 fc'(e" 2kxq -l)M(k,k',q,w)(e- 4z2k ' xq -l) (174) 

where the factors at each end reduce to dipoles at small ql and M = M + Mj- is a sum of two terms, one from the 
free Fermi sea and one from the exchange of the transverse collective mode which peaks at iui ~ q 3 v(q). At high 
frequencies, Mr can be ignored and the dipoles emerge as free objects, while at low frequencies, they are coupled 
by the overdamped mode and do not behave like classical dipoles. We begin to see how the dipoles appear in the 
correlators (because we can tune the frequency and wavelength in K to expose them) but not the wavefunctions, 
which describe equal-time correlations and thus involve an integral over all frequencies. 

Arguments for the compressibility of the dipolar system were also given by D.H. Lee (1998). 

The resolution of the compressibility paradox did much to assure the community that various descriptions of the 
quasiparticle, each tailor-made for a different occasion, were mutually compatible and consistent. 
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E. Higher Landau Levels 



The advantage of the extended hamiltonian is that it keeps track of the electronic cyclotron coordinate and does 
not go to the LLL prematurely. This has many benefits. If we want to compute the Hall conductance, we can couple 
the system to an external potential and find the response. As stated previously, the response involves higher LL's so 
that the presence of T] e is crucial. The details are shown in Appendix IB! We can also study the effects of LL mixing. 
Instead of using just Hlll we can try to get an effective theory within the LLL which subsumes the effects of virtual 
transitions to higher LL's. 

To extract the leading correction due to higher LL's we write the Schrodinger equation in schematic form as as 



n 



E 



(175) 



where <f> is restricted to the space spanned by Fock states composed of just the LLL states and £ stands for everything 
else. Likewise the subscripts and n' stand for collective labels in the LLL and above the LLL respectively. The 
exact equation obeyed by (f> is 



Ha 



1 



-H n 



= E4> 



E — H n ' n 

which is not an eigenvalue problem since E appears on both sides. However we may approximate as follows 

1 1 



(176) 



E — H n ,i, 



H n .lr 



O(v/uo) 



since the eigenvalue E we are interested in is of order v and the eigenvalues of H n 
accuracy in k — v/lu we can also replace 



a' a a a n a UJ 



(177) 

are of order u>q. To the same 
(178) 



which leads to 



h; 



00 ^ — -^oo — H 0n / — H n '0 

•<- — J n't, in 



n'ujQ 



(179) 



Once we have the effective theory in the LLL, we can switch to the CF formalism: Introduce R„, exchange R e and 
R„ for R and r\ and proceed as before by a Hartree-Fock calculation. The results are in accord with earlier works 
that showed that LL-mixing reduces the transport gap, but that finite thickness reduces this effect (Yoshioka, 1986, 
Price, Platzman, and He, 1993, Melik-Alaverdian and Bonesteel, 1995, Price and Das Sarma, 1996, Melik-Alaverdian, 
Bonesteel, and Ortiz, 1997). Full details can be found in Murthy and Shankar (2002). 



VI. CORRELATION FUNCTIONS IN THE CONSERVING APPROXIMATION 

This section illustrates how one does the conserving calculation within the hamiltonian approach by deriving 
density-density correlation functions (Murthy, 2001a). This calculation differs from Read (1998) in two ways. First, it 
is done for nonzero effective field [y ^ 1/2), so that the CF's are in Landau levels instead of plane wave states. Next, 
the calculation is done in the operator approach (Anderson, 1958a, b, Rickayzen, 1958) using equations of motion 
as compared to diagrams (Nambu, 1960). The second difference is only cosmetic, and introduced here to promote 
harmony with the rest of the paper 33 . 

The Hamiltonian and constraint to be solved are 

H = ~ £^(q) %)e"W 2 /2 j5(_ q) (1 80 ) 
q 

X(q) * 0. (181) 

Let us define a time-ordered pseudovortex-electron density-density Green's function G ve (q,t) as follows (with G ee 
and G vv similarly defined): 



For an illustration of calculations in the gapped fractions using the diagrammatic approach, see Green (2001). 
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which evolves as per 



G ve (q,t-t') = -*<Tx(q,t)p(-q,t , )> 

= -iQ(t-t') <x(q,*)p(-q,0 > 

- iG^-t) <^(-q,Ox(q,*) > (182) 



-i^-G^t-t') = -8(t-1?) < R(q,i),p(-q,0] > 



St 



- ;<T[if,x(q,iM-q,0 > (183) 



Since % commutes with one immediately sees that G ve is a constant. If the constraint is set to zero initially, then 
it remains zero, and all its correlators also remain zero. 

The above is true in an exact treatment of the theory. Of course, we can only do approximate calculations. A 
calculation that respects G ve — G vv — respects the symmetries of the theory at the level of correlators, and is called 
conserving. Let us see what a natural approximation scheme might be. Consider 

d 

i—G ee (q,t) = 5(i)<[£(q,t),£(-q,0)]> 

- !<T[F^(q,#(-q,0)> (184) 

Since [p , p] ~ p, a Green's function involving three densities will arise. Additional time derivatives will produce 
higher-order density correlators, leading to a hierarchy of equations for more and more complicated Green's functions. 
The natural way to truncate this hierarchy is to make a mean-field approximation at some stage that reduces a 
product of two densities to a single density. One of the simplest of such approximations (Baym and Kadanoff, 1961) 
reduces [H, p], which is a product of four fermi operators, to a product of only two, by using the averages 



A^dfadpi -> < > d l 2 d f>2 + < d l 2 d f>2 > d iAh 

- < d lM > d lM- < 4, 2 d 0i > (185) 

Here < d) a dp >— 5 a pNF{a), where Nf{(x) is the Fermi occupation of the single-particle state a. 

Using the HF states and their occupations in the above truncation is known as the time-dependent HF (TDHF) 
approximation 34 . We will use the operator approach to TDHF as expounded by Anderson (1958a, b) and Rickayzen 
(1958). We will explicitly see below that it is conserving for all principal fractions. 

The physical picture underlying our calculation is the following. When a bosonic operator such as p(q) or x(q) acts 
on the ground state, it creates a linear combination of particle-hole pairs. In the Landau gauge, each pair is labeled by 
an index v — m, rii (not to be confused with the filling factor!) that keeps track of the CF-LL indices of the particle 
and hole, and a total pair momentum q which is conserved because the particle and hole have opposite charges and 
do not bend in the magnetic field. It is clear that in order to calculate time-dependent response functions we have to 
understand the time-evolution of these pairs. In the exact theory, the Hamiltonian can scatter a pair into two pairs, 
two pairs into four pairs, etc. This is what leads to the hierarchy of equations. However, the great simplicity of the 
TDHF approximation is that in this approximation a particle-hole pair scatters only into another particle-hole pair. 
At a given q we thus have a matrix labeled by the indices (v, v') of the incident and scattered particle- hole pairs. 
The magnetoexciton spectrum comes from the eigenvalues of this matrix, which, together with its eigenvectors, will 
be seen to explicitly determine the Green's function. 

Clearly, we must begin with an operator that creates an exciton in a state of definite momentum q, starting with 
states labeled by CF-LL indices v = (m n 2 ). The following operator does the job: 



(186) 



The TDHF approximation has been the method of choice in computing magnetoexciton dispersions in the IQHE. Here are some of the 
early references: Chiu and Quinn (1974), Horing and Yildiz (1976), Theis (1980), Bychkov, Iordanskii, and Eliashbcrg (1981), Bychkov 
and Rashba (1983), Kallin and Halperin (1984), MacDonald (1985), Hawrylak and Quinn (1985), Marmorkos and Das Sarma (1992), 
Longo and Kallin (1993). 
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Why is this so? First note that in this gauge, (Eqn.||lSJ|), the wavefunctions are plane waves in y with momentum 
k and localized in x at X — kl* 2 . Thus d' ,, 2 d „ ;»2 which creates a hole at X — qy „ and a particle at 

X + Qyl 2 , creates an exciton in a state of momentum q y in the y-direction, centered at X. Multiplying by e~ lQxX 
and summing over X creates a state of momentum q x in the x-direction. 
Now we define 

G{u- z/; q; t - t') = -i < TO v (q, t)O v , (-q, t') > . (187) 

Taking the time derivative we get 

-iJ^C?0/;i/;q;t) = -i < T[H,O v {q,t)]O uf (-q,0) > 

-S(t) < [O u (q, t), O u » (-q, 0)] > (188) 

The last piece is the standard inhomogeneous "source" term. The dynamics is controlled by the commutator with 
the Hamiltonian, to which we now turn 

[H,O v (q)] = (e(ni)-e(n 2 ))0 I/ (q) + 

(N F (n 2 ) - N F (m)) £ {-^$Y e ' q2l2/2 P^ 2 ^) x ^m(-q) 

- / 7^2^( s )^ s2i2/2 ^> 1 ( s )Pn 2 ^(- s )^ (r)2sXq )^(q) (189) 

where e(v) is the Fock energy in a state of CF-LL index v, and the TDHF approximation has been made. 

Note that the action of commuting with H on O v in the TDHF approximation can be represented as the right- 
multiplication by a matrix TL{v; \ q). It follows that if a we form a linear combination of operators = 
^(v, q)0„, the column vector * will transform linearly under the action of Ti. Diagonalizing Ti will also en- 
able us to solve for the Green's function, for the eigenvectors represent linear combinations of particle-hole states that 
are normal modes of Ti.. Assume that one has found the right and left eigenvectors and corresponding eigenvalues 
E a (q) , labeled by a 



W(i/;i/";q)**(i/";q) = E a (q)^(u;q) (190) 
vP^;q)H(^";q) = E a (q)^(v»; q) (191) 

where sums over repeated indices are implicit. 

Assuming the matrix Ti has a complete set of eigenvectors the Green's function can be written as 

G(i/;i/;q;0) = £ ^(mttj; q)^^^"!; q)(iV F (ni) - N F (n' 2 )) (192) 

where the factor (Npin^) — Np(n' 2 )) comes from the source term (the last term of Eq. I|188ll '). 

An important property of Ti. from the point of view of the conserving approximation is that it always has one left 
eigenvector with zero eigenvalue for every q, namely 

*o( n i"2;q) = Xn in2 (q) (193) 

where 

Xn 1 n 2 (q) = <rii|e-^W c )|n 2 ). (194) 

The existence of this eigenvector with zero eigenvalue (for any v(q)) is shown explicitly in Murthy (2001a). This 
zero eigenvalue is related to the constraint and is one of the conditions for the TDHF approximation to be conserving. 
To see this let us go back to the exact theory. There the Hamiltonian commutes with x(q). This means that x(q) 
acting on the ground state should produce a zero energy state at every q. In the TDHF approximation the action 
of the Hamiltonian on this state is approximated by the action of the matrix Ti on the left vector q>Q (nir^; q), 
corresponding to the operator 0^( q ) = J2Xn 1 n 2 { c L)O ni n 2 (q). The fact that H admits vPq as a left eigenvector with 
zero eigenvalue means that the zero energy state which had to be present in the exact theory is also present in the 
TDHF approximation. 
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To put it in slightly different words, the TDHF approximation represents a truncation of the Hilbert space of 
neutral excitations to states having only a single particlc-holc pair above the ground state. A priori this truncation 
need not have respected the constraint, but it does. 

The above condition is necessary for TDHF to be conserving, but not sufficient. The other condition is that the 
physical sector (excitations created by p(q)) does not couple to the zero energy sector of TL. This can also be verified 
in a straightforward way (Murthy 2001a). 



A. Small-q structure factor 

We will now put the conserving approximation to another test. Recall that the density-density correlator has two 
parts; one coming from the cyclotron pole, with a residue of q 2 , which is unrenormalized as per Kohn's theorem and 
saturates the sum rule, and another coming from the dynamics in the LLL. Hence the LLL-projected structure factor 
S(q) has to vanish faster than q 2 , and if it is analytic in q 2 , it has to go like q A . We will see if this is true. Luckily, to 
obtain the leading behavior of S(q) in the TDHF approximation one needs to keep only a finite-dimensional submatrix 
of the infinite TDHF matrix. 

The results are as follows. For v = ^ we have, upon diagonalizing a 4 x 4 matrix, 

s(g) = ^ + - (195) 

We can extend this result to all the Laughlin fractions ^^px to obtain 

3(q) = l(ql)* + ... (196) 

The coefficient |, independent of s, differs from the result | obtained from Laughlin's wavefunction (Girvin, Mac- 
Donald, and Platzman, 1986). However no general principle requires that we regain this coefficient. 

One can carry out a very similar, but much more tedious, calculation for all the principal fractions (Murthy, 2001a). 
The result, upon diagonalizing a 6 x 6 matrix, is 



( q Z) 4 P A ~ 3p 3 + \y 2 + 3p + 7 - 
p 2 -l 



S(q) = 7 4 \ (197) 



This expression, with its divergence as p — ► oo or v — ► ^ is consistent with the result (Read, 1998) that for a problem 
equivalent to the v — \ problem, S(q) ~ q 3 \og(q). As p — > oo, the radius of convergence of the power series expansion 
of S(q) must go to zero (else the structure factor would diverge for a range of q according to Eqn. I|197|) l. and the 
p — > oo limit does not commute with the q — > limit. The formula also does not hold for p = 1, the Laughlin fractions, 
for which we must use the results quoted earlier. 

Surprisingly, the results are independent of the form of the potential, which does enter the intermediate stages of 
the calculation. 



B. Magnetoexciton dispersions for | and | 

We saw that the small-g behavior of S(q) can be satisfactorily addressed using the TDHF approximation. It turns 
out that this approximation also works well for computing the dispersion of magnetoexcitons, by which we mean the 
lowest energy physical eigenstate of TL at each q. 

As shown in the previous section, for very small q, the naive magnetoexcitons (bare particle -hole states created 
by O v ) do not mix with others and become the true eigenstates of TL. As q increases they become increasingly 
coupled, and both level repulsion (between positive energy states) and level attraction (between positive and negative 
energy states) manifest themselves, giving rise to complicated magnetoexciton dispersions. The matrix TL is infinite- 
dimensional, but in any numerical calculation, only a finite matrix can be diagonalized. We saw that the lowest 
nontrivial result for S(q) could be obtained by keeping at most a 6 x 6 matrix. As q increases more and more CF-LLs 
have to be kept to obtain accurate results. The accuracy of the truncation was checked by two different methods. 
Firstly, the number of CF-LLs kept was increased until the energy of the magnetoexciton was stable. Secondly, we 
knew that at every q there should be two zero eigenvalues corresponding to the unphysical sector. The number of 
CF-LLs kept in the calculation was increased until these null eigenvalues were at least four orders of magnitude smaller 
than the smallest physical eigenvalue. 

Figures 4 and 5 are for the potential v(q) = ^e"? 13 I 2 for various values of (3. It is worth noting that though 
the depths of the minima are /3-dependent, and different from what are found in numerical diagonalizations or from 
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FIG. 4 The lowest-energy spin-polarized magnetoexciton for v — 1/3 for two values of (5. Note the " magnetoroton" at ql = 1.4. 
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FIG. 5 The lowest-energy spin-polarized magnetoexciton for v — 2/5 for two values f3. 
ql = 1.5. 



Note the minima at ql = 0.8 and 



CF-wavefunctions, the positions are correct 35 . This indicates that the TDHF approximation to the EHT does capture 
the important physics of the electronic problem even at fairly large q. 

In the case of |, the magntoexciton energy has stabilized at large ql enabling us to read off the activation gap, 
A a , defined as the minimum energy needed to produce a widely separated particle-hole pair. However, the situation 
deteriorates rapidly as p increases. It becomes prohibitively hard, even at p = 2, to get the large ql limit of the 
magnetoexciton spectrum, as is evident from Figure 5. This is because the decoupling of the naive magnetoexcitons 
is only asymptotic for large q. However, as p increases, the shortcut gives very good answers (Figures 6 and 7) with 
hardly any additional work. We now turn to it. 



35 Haldane and Rezayi (1985), Su and Wu (1987), Morf and Halperin (1987), d'Ambrumenil and Morf (1989), He, Simon, and Halperin 
(1994), He and Platzman (1996), Kamilla, Wu, and Jain (1996a,b), Jain and Kamilla (1998). 
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VII. GAPS 



Having developed the EHT at length and having established several qualitative results and matters of principle, 
we now turn to mainly quantitative issues. As mentioned at the outset, one can compute just about anything within 
this approach, some more accurately than others. In the next few sections we will discuss a variety of such quantities. 
Rather than run through an endless list we will focus on a few that give the flavor of the method and expose its 
strengths and shortcomings, and point out references containing more details or examples. An invaluable benchmark 
will be provided by results from trial wavefunctions and exact diagonalization. Comparison with experiment will be 
done with the clear understanding that no systematic attempt is made to incorporate disorder. 

In comparing to the results of exact diagonalization and trial wavefunctions one must bear in mind their limitations. 
Exact diagonalization (ignoring machine errors) suffers from the fact that the systems are necessarily small. In 
comparing these results to ours (valid for infinite systems) we must examine the approach to the thermodynamic limit 
and the scaling of gaps. We refer the reader to Morf, d'Ambrumenil, and Das Sarma (2002) for the most recent study 
of this kind. 

While trial wavefunctions can be written down for any number of particles, the evaluation of gaps requires once 
again that the number of particles be small (though not as small as in exact diagonalization). The details may be 
found in Park, Meskini and Jain (1999) to whom we compare our numbers. It must also be borne in mind that the 
correct wavefunctions may not be of the form being tried, in which case even a flawless evaluation of gaps is irrelevant. 
So far there is good reason to believe that for 1/3 < v < 1/2 the Jain- like functions are good. At very low densities 
the correct state is believed to be a Wigner crystal and when higher LL's are involved, more complicated functions 
like Pfaffians (the Moore-Read, 1991, state) need to be explored. With these caveats in mind we proceed with the 
comparison of our theory with numerics. 

In the remaining sections, we will use the preferred density and pay no further regard to constraints. As stated 
earlier, this is the most efficient way to work on problems that do not depend crucially on the deep infrared region 
(uj ~ q 3 v(q)). In the cases discussed this is ensured by either a gap, a nonzero temperature, or both. 

Let us begin with the activation or transport gap in a fully polarized sample, defined as the minimum energy needed 
to produce a widely separated particle hole pair: 



where |p) stands for the Hartree-Fock ground state with p-filled LL's and PH stands for a widely separated particle- 
hole pair. We will use a boldface symbol such as p to label a Slater determinant with p occupied Landau Levels. 
Non-boldface symbols will label single-particle states. Note also that the highest occupied CF-LL index n for the state 
labeled by p is n = p — 1 since the CF-LLL has index n — 0. As shown in Appendix [U] the particle-hole excitations 
of |p) are HF states of our H. 

The expression for the gap written above is formally the same in the wavefunction based approach of Park, Meskini 
and Jain (1999) to whose results we will compare ours. However the notation hides a big difference. They work 
in the electronic basis where p(q) = cxp(— zq • r^), and the states are the simple wavefunctions \ v , multiplied 
by the Jastrow factor and then projected to the LLL. (Projection leads to a very complicated expression for the 
wavefunctions.) In the present approach we have tried to incorporate these effects by going in the reverse direction, 
from electrons to CF's, and obtaining complicated expressions for the charge and other operators, but with simple 
expressions for the wavefunctions. While these operator expressions are unusual in form, they are simple to evaluate 
within the HF calculation. 

In all our calculations we shall use the Zhang-Das Sarma potential (Zhang and Das Sarma, 1986) 



which is a crude model for the electron-electron interaction in a sample of finite thickness. We simply take it to be 
a one parameter family of potentials. In terms of the Haldane pseudopotentials V m (which gives the interaction in a 
state of relative angular momentum m, see Haldane, 1987) we know that just one (typically Vi) dominates. We can 
think of A as controlling the operative pseudopotential. 

Rather than work with a widely separated particle-hole pair, we find the energy in a state with just the particle 
and add to it the energy of a state with just the hole and subtract double the ground state energy. Relegating the 
details to Appendix El we present the central idea. 

We begin with the second quantized expression for the preferred charge operator p p (q): 



A = (p + PH\H\p + PH) - (p\H\p) 



(198) 



27re 2 e-?' A 



(199) 



vzds = 





(q) 



(200) 
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FIG. 6 Comparison of dimensionless activation gaps S a to the work of Park, Meskini and Jain (1999) for the fractions |, |, |, | 
(p = 1, 2, 3, 4, and s=l) as a function of A, the thickness parameter in the ZDS potential. 

where d\ nn creates a particle in the state \m n) where m is the angular momentum and n is the LL index of CF in 
the weakened field A* = A/(2ps + 1) with a magnetic length 



I* = ly / 2ps + 1. (201) 

The key ingredient in the HF calculation is the matrix element p m2 „ 2 - raini which factorizes (as shown in Appendix 
EI): 



Pm2ii2;miiii Pm 2 mi ® Pnin^ (202) 

The gaps depend only on n , the superscript on which will be generally dropped. Often we will use the 
dimensionless activation gap 5 a defined by 

A Q = 4 S a - (203) 

Figure 6 shows the gaps computed for §, §, § and | for the ZDS potential and compared to the work of Park, 
Meskini, and Jain (1999) in the region < A < 3. The following features are noteworthy. 

• At A = 0, the Coulomb case, the gaps are finite in contrast to the small-g theory (Murthy and Shankar, 1999). 
This is due to the gaussian factor e^ 2 ' 2 / 2 in H(p) ((Eqn. (EJ) which was absent in the small-q theory. The 
slope of the graphs in the present theory is nonzero at this point. It is readily verified that dA/dX at A = 
is the gap due to a delta-function potential, and should vanish for spinless fermions. The description of CF 
in terms of p p does not give good answers for potentials as short ranged as the delta function, as anticipated 
earlier. Indeed even the Coulomb interaction is too singular, and the theory begins to work well only beyond 
A~l. 

• Beyond A ~ 1 the agreement is quite fair in general and best for | . However we cannot go to too large a A since 
in this range the system may not be an FQHE state. 

• The gaps which do not vanish for any fraction and any finite A, exceed the Park, Meskini, and Jain (1999) 
values for | and | and lie below them for | and | . This result is at odds with the general belief that HF always 
overestimates the gaps by neglecting fluctuations. 
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FIG. 7 Comparison of hamiltonian theory to the exact diagonalization results of Morf et al for p = 2 and 3 . 

The situation is different when we compare with the exact diagonalization results of Morf, d'Ambrumenil, and Das 
Sarma (2002) who used a potential 

v fq) = ^- e (?' fc ) 2 ) Erfc (qlb) (204) 

q 

where b is the analog of A 36 . Our numbers are compared in Figure 7. The calculated gaps always lie above the exact 
diagonalization results for the two fractions shown (as well as for the i case, not shown). 

The general disagreement with our theory is worse for this potential than for the ZDS case because at large q this 
potential goes as l/q while the ZDS potential falls exponentially. 



A. Activation gaps 

We have computed gaps for many other fractions, including for s = 2, when four vortices are attached to form 
CF's. Rather than show more plots, we will now analyze the theory in terms of the activation mass m a defined by 

A - eB * eB (2 5) 

m a (2ps + l)m a ' 

Comparison to Eqn. (|203[l shows that 

1 e 2 l e 2 l 
— = —6 a (2ps + 1) = — C a . (206) 

m a e e 



Thus 



C a =5 a (2ps+ 1). (207) 



Based on the exact diagonalization results of d'Ambrumenil and Morf (1989) Halperin, Lee, and Read (1993) 
pointed out that C a approaches a limit as we approach v — \ or p — > oo. The Halperin, Lee, Read theory expects C a 
to be modified by logarithms. We too expect these logarithms once we approach v = \ and are forced to include the 
overdamped mode. However, we find as did Halperin, Lee, and Read (1993), a good fit to the calculated gaps without 
invoking the logarithms, which are operative in a very tiny region near v = i (Morf, d'Ambrumenil, and Das Sarma, 
2002). The only difference here is that C a approaches a limit that depends on A, a parameter they set equal to zero 
(Coulomb case). 

We next turn to Pan et al (2000) whose experiments detected that the normalized mass defined by 

m 

m nor = _a (208) 

where m e is the electron mass and B(T) is the field in Tesla, is nearly the same for s = 1 and s = 2, i.e., two and four 
flux tubes. The theory predicts that m nor are comparable for s — 1 and s = 2, but makes it clear that no fundamental 



This potential is not very different in form from the ZDS potential except at very short distances. 
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TABLE II Activation gaps as a function of A for 1 < A < 2 according to the hamiltonian theory. Note the convergence of Ca' 
as p — > oo. The superscript (2) refers the number of flux tubes attached. 



p 




S {2) 


cf 


1 


5.31^B(T)/A 


.106/A 


.32/A 


2 


2.0Sy / B(T)/\ 


.042/A 


.21/A 


3 


1.23v / B(T)/A 


.025/A 


.17/A 


4 


0.87 1 [/B(r)/A 


.017/A 


.16/A 



TABLE III Approximate numbers used in this paper, with hs the Boltzmann's constant, m e the electron mass, and B(T) the 
field in Tesla. 



l.UB(T)°K 



50y/B(T)°K 
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significance can be attached to this result since it depends on A, or more generally, the potential. Further details may 
be found in Shankar (2001). 

Consider next the experiments of Du et al (1993), who have extensive data on activation gaps. We will limit 
ourselves to v < i, to which states with 1 > v > i are related by particle-hole symmetry if full polarization is 
assumed. Given that the experiments, unlike Park, Meskini, and Jain (1999), have an unknown contribution from LL 
mixing and impurities, it is not clear how to apply the theory. There is no ab initio calculation that includes these 
effects. (There is however reliable evidence that LL mixing is a very small effect at the values of A under consideration. 
Recall also our results from Section [V]) 

We will compute gaps using the ZDS potential with A as a free parameter, and ask what A fits the data, just to 
get a feel for its size. The results are summarized in Table llVl 

Comparing the above values of A extracted from data to the Local Density Approximation (LDA) (Price and Das 
Sarma, 1996, and references therein) and exact diagonalization calculations (Park, Meskini, and Jain, 1999, Morf, 
d'Ambrumenil, and Das Sarma, 2002), which suggest A ~ 1, we see that disorder has a substantial effect on activation 
gaps. 

It is possible to compute the charge density in a state with a widely separated particle-hole pair in some gapped 
fractions. The details, and a comparison to the unpublished work of Park and Jain may be found in Shankar (2001). 



B. Other potentials 

Figure 8 shows a comparison to the Park, Meskini, and Jain (1999) results for a gaussian potential 

v(q) = 2ire 2 l e~ qH2/2 (209) 

Note that except for v = | the agreement is exceptional. This is the kind of ultraviolet-soft potential for which the 
present theory works best, though unfortunately we do not presently know of any system where it is operative. At 
present, its study serves to instruct us on the domain of validity of our approach. Not surprisingly, work on the 
potential, 

v(r) = - (210) 

r 

shows an agreement that is worse than for the Coulomb case since this potential is just as bad as r — * and does not 
give the large r values a chance. Likewise 1/r 2 fares worse than 1/r. 

Note that if A — > oo we may be in trouble since the potential may not lead to FQHE states. 



TABLE IV Comparison of activation masses to Du et al, sample A, which has a density n = 1.12 • 10 11 cm 2 . The last column 
gives the best fit to A. 
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FIG. 8 Comparison of gaps for the gaussian potential. For p > 1 it is hard to separate the results from the Hamiltonian theory 
and the wavefunction results of Park et al. 

In summary, it appears as if our approach works well in any problem which is not sensitive to distances smaller 
than about a magnetic length. Since the extended formalism is mathematically equivalent to the original electronic 
problem in the LLL; it is the twin approximations the use of p p to deal with constraints and Hartee-Fock - that are 
responsible for deviations from our benchmarks. 

While our approach can reproduce the numbers of the wavefunction and exact diagonalization approaches to within 
about 10%, (if A > 1), with a lot less work (a few seconds on a PC), without these benchmarks, we would have known 
neither its range of validity nor its degree of accuracy. Without such feedback, we would not have applied the method 
with the such confidence to other phenomena not treated by wavefunctions or exact diagonalization such as relaxation 
rates and polarization at T > 0. 

VIII. MAGNETIC TRANSITIONS AT T = 0. 

Now we turn to the behavior of the spin of the system, which so far was assumed to be frozen along the applied 
field. The reader might wonder why we bother, since the magnetic fields are of order at least a few Tesla. The reason 
(originally pointed out by Halpcrin (1983)) is that the Zeeman energy Ez = g* (^)B is the smallest energy scale in 
the problem. This is because of a conspiracy of two factors: First, due to band structure effects, the band mass of the 
electron in GaAs is 0.068m e , with m e the electron's mass in vaccuum. This makes the cyclotron frequency about 14 
times what it would have been in empty space. Next, due to spin-orbit coupling, the effective g-factor of the electron 
in GaAs is g = 0.44 instead of the empty space value of 2. Thus, the Zeeman energy is about 64 times smaller than 
the cyclotron energy ojq. Typically the interaction energy is of the same order as ojq, leading to Ez being the smallest 
energy. Since the original realization by Halpcrin (1983), spin transitions have been seen in experiment 37 , and been 
explored by exact diagonalization 3S . We will compare our results with the recent work of Park and Jain (1998, 1999) 
in the CF-wavefunction approach. 

The coupling of electron spin to the applied field is given by the Zeeman term 



where n is the density and P is the polarization, to which each electron contributes ±1. 

Since the uniform external field couples to the q — component of the spin density which is unaffected by the 
canonical transformations, Hz will have the same form in the final CF representation. 

A. Magnetic transitions in gapped fractions 

When Hz is large, we expect the system to be fully polarized (P = 1). As we lower Hz, we expect P to drop. If 
translationally invariant CF states are considered for the gapped fractions, there is a discrete set of allowed values of 



Clark, Haynes, Suckling, Mallett, Wright, Harris, and Foxon (1989), Eisenstein, Stormer, Pfeiffer, and West (1989), Furneaux, Syphcrs, 
and Swanson (1989), Buckthought, Boulct, Sachrajda, Wasilewski, Zawadski, and Guillon (1991), Du, Yeh, Stormer, Tsui, Pfeiffer, and 
West, (1995,1997). 

Chakraborty and Zhang (1984a,b), Rasolt, Perrot, and MacDonald (1985), Chakraborty, Pietilainen, and Zhang (1986). 




(211) 



where g = .44, m e is the electron mass in free space, S is given by 



S = nP 



(212) 
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P. At v = pj (2ps + 1), these correspond to states of the form |p — r. r) in which p — r CF-LL's are occupied by up 
spins and r CF-LL's by down spins. 39 Thus the allowed values of polarization are given by 

V — 2r , 

P=- . (213) 

P 

For example, when p = 4, the allowed values are P = 1, .5, and corresponding to |4, 0) ,|3, 1) and |2, 2). 

Our goal is to calculate the critical fields at which the system will jump from one value of r to the next as Hz is 
varied. Let 

E(p-r,r) = (p-r,r| J ff|p-r,r) (214) 

where H does not contain the energy due to Hz- This will be the case for the single-particle and ground state energies, 
with one exception which will be clearly pointed out. Since Hz is diagonal in the HF states which have definite spin, 
its effects can be trivially incorporated. 

The HF calculation of E(p — r, r) is detailed in the Appendix IFl 

The critical field B c for the transition from r to r + 1 is given by : 

eP> c ti 

E(p-r,r)-E(p-r-l,r+l)=g- (215) 

2m e p 

where the right hand side denotes the Zeeman cost of flipping the n/p spins in the LL that switched its spin. This 
discussion assumes that B is perpendicular to the sample. If there is a tilt 9, we write 

E(p-r,r)-E(p-r-l,r + l)=g &B1 - (216) 

2m e cos u p 

When these energy differences are calculated, the same remarkable regularity first noted by Park and Jain (1998, 
1999) emerges: they can be fit by a theory of free fermions of mass m p the polarization mass, that occupy LL's with 
a gap A p = eB* /m p . By this we mean the following. In a free theory of gap A p , we would have 

E{p-r,r) - E(p-r-l,r+l) = n( - p ~ 2r Z 12 a (217) 

P 

since (n/p) spin-up fermions of energy (p — r — 1 + ^)A p drop to the spin-down level with energy (r + |) A p . Suppose 
we evaluate the left-hand-side of Eqn. 1217|) in the HF approximation and define 

A{r)def = ££?( P -r,r)-E(p-r-l,r+l)_ 
n p — 2r — 1 

Given that H is not free, there is no reason why A p (r) de f should be r-independent. But it is very nearly so. For 
example at p = 6, A = 1, 

e 2 

A p (0,l,2) rfe/ = —(0.00660,0.00649,0.00641) (219) 

which describe the transitions |6, 0) — > |5, 1), |5, 1) — > |4, 2), and |4, 2) — > |3, 3). This r-independencc of the gaps was 
true for every fraction and every value of A we looked at. We will soon demystify this apparent free- field behavior, 
which has a counterpart in the gapless case as well. 

The transition |p — r, r) — > |p — r — 1, r + 1) occurs when the spin-flip energy equals the energy difference of the 
two competing ground states: 

<?^^ = (p-2r-l)A p . (220) 
2m e cos 



39 In Jain's approach, the actual wavefunction will be such a state times the Jastrow factor, followed by projection to the LLL. In the 
present approach, p — r, r) is literally the state, but the operators for charge and spin are obtained by transformations to the CF basis. 
For the interested reader we mention that these transformations are spin-independent. 
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B. Magnetic properties of gapless fractions 



Let us now turn to the gapless fractions \ and \ . The discrete labels p — r and r of the HF states that count the 
spin-up and down LL's are now replaced by continuous variables k±p which label the Fermi momenta of the spin-up 
and down seas. These momenta are such that the total number of particles equals n: 



k\ F + k 2 _ F = kp = 47m 



(221) 



where kp denotes the Fermi momentum of a fully polarized sea. 

In the gapped case there were several critical fields B c , each corresponding to one more CF-LL flipping its spin, 
each describing one more jump in the allowed values of P. In the gapless case the situation is different. For very large 
Zeeman energy, the sea will be fully polarized. It will not be worth including even one fermion of the opposite spin 
since the Zeeman energy cost alone will exceed the Fermi energy of the polarized sea. As we lower the Zeeman term, 
we will reach a critical field at which it will be worth introducing one fermion of the other spin with zero (effective) 
kinetic energy. At this point the energy of a particle on top of the spin-up sea obeys 



£ + (k +F ) = g 



2m„ cos 6 



(222) 



If we lower the Zeeman term further, the polarization will fall continuously and be determined by £±(k± F ) 1 the 
energies of the particles on top of these two seas according to 



£+{k +F ) - £-(k- F ) = g 



e B ± 
2m e cos 9 



(223) 



This equation states that the system is indifferent to the transfer of a particle from one sea to another, i.e., has 
minimized its energy with respect to polarization. 

Since the effective magnetic field vanishes at the gapless fractions, we deal with a very simple expression for p(q) p : 



It is easy to do a HF calculation and obtain 



(224) 



£ ± {k)) = 

~ [ d 2 q , , . 



sin 



kxql 2 



-4 



d 2 k' 



4tt 2 
£o+£i 



n£(|fc'|M|k-k'|)sin 2 



k'xk! 2 



where the Zeeman energy is not included, n± is the (step) Fermi function for the two species, £ and £j represent 
single particle energy (due to what was called H earlier) and the energy of interaction of this particle at the Fermi 
surface with those inside the sea, and 



v{k) = v{k)e- kH2 12 . 



(225) 



When this result is used to compute £+(k+ F ) — £-(k- F ), we find once again that the numbers fit a free theory in 
the following sense. Imagine that CF were free and had a mass m p . We would then have 



£+(k +F ) - £_(k- F ) 



k 2 
+F 



k 2 



2m E 



(226) 



What we find is that the HF number for £{k + p) — £(k-p) may be fit very well to the above form with an m p that is 
essentially constant as we vary k±p i.e., the relative sizes of the up and down seas (which is analogous to an m p that 
docs not depend on the index r in the gapped case). This constant m p defined by 



1 



= 2 



£+(k +F ) -£-{k_ F ) 



k 2 

H + p 



k 2 



(227) 



matches smoothly with that defined for the nearby gapped fractions. 
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This free-field behavior is surprising because we will shortly see that there are many reasons to believe that the 
CF's are not free. For now, note that the HF energies are not even quadratic in momenta: for example at v — i and 
A = 1 there is a hefty quartic term: 

£{k±F) =a(^)\b('^) i (228) 



(e 2 jel) \ ftF J \ kF 
where a = .075, b = - .030. 

There is no reason the CF kinetic energy should be quadratic in momentum. These particles owe their kinetic energy 
to electron-electron interactions, and given this fact, all we can say is that their energy must be an even function 
of k, starting out as k 2 at small k. What constitutes small k is an open question that is answered unambiguously 
here: our expression of the energy has substantial fc 4 terms for momenta of interest, but not from higher powers. An 
analogous result holds for the gapped fractions, where the CF-LL's are not equally spaced (Murthy 1999, Mandal and 
Jain, 2001a) 

The proper interpretation of this free- field behaviour will be taken up next. 

C. Composite Fermions: Free at last? 

The fact that magnetic phenomena at T = can be described (to excellent accuracy) by free fermions of mass m p 
(Park and Jain, 1998, 1999, Shankar, 2000, 2001) needs to be properly understood and interpreted. In particular, 
one must resist the thought that perhaps by some further change of variables one could take the present hamiltonian 
and convert it to a free one. This is because if there were really an underlying free theory, it would have a single 
mass mcF for both activation and polarization phenomena, with the CF's forming LL's of spacing eE>* jnicF- But 
we know from the extended hamiltonian theory, Jain's approach, or experiment, that there are two masses m a and 
m p that differ by at least a factor of two. Furthermore the shape of the magnetoexciton dispersions (Kamilla and 
Jain, 1998, and Figure 5 also points to sizeable CF interactions: as ql 2 , the distance between the particle and hole 
varies, the (binding) energy varies by an amount comparable to their individual energies, whose sum is given by the 
value at large ql. 

It has been shown (Shankar, 2000, 2001) that a single assumption about the form of the ground state energy, 
an assumption that is not equivalent to the free-field assumption or even to a quadratic dispersion relation in the 
gapless cases, will explain this behavior for gapped and gapless fractions. Consider E(S), the ground state energy as 
a function of S = nP, where P is the polarization. By rotational invariance it must have only even powers of S in its 
series. Assume the series is dominated by the first two terms: 

E{S) = E(0) + ^S 2 (229) 

where a is the inverse linear static susceptibility. 

Consider first the gapless case. When dn particles go from spin-down to spin-up, 

dE = a S dS = a S {2 dn) (230) 



Ot- 



ic- 2 — h 2 
+F ^—F 



47T 



(2 dn) (231) 



using the volumes of the Fermi seas. We see that dE has precisely the form of the kinetic energy difference of particles 
of mass nrip given by 

— = -. (232) 
m. p 7T 

Thus nip is essentially the static susceptibility, which happens to have dimensions of mass in d = 2. The statement 
that m p has no r-dependence in the gapped case or no spin dependence in the gapless case is the same as saying that 
the full nonlinear susceptibility does not depend on the spin S, which in turn means E(S) is quadratic in S. 

Note that the free-field form of dE comes from E ~ S 2 and d = 2: in d = 3, we would have dE/dn ~ S ~ 
(fc+F — ^-f) which no one would interpret as a difference of kinetic energies. 

Let us see how this general argument applies to the specific example we have been working on. 

Consider the HF energies quoted earlier 
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The quartic terms miraculously drop out in the energy cost of transferring a particle from the top of the spin-down 
sea to the top of the spin-up sea: 

a 3 + b -. (234) 



(e 2 /el) k% k F 



k+F k_ F ,{k 2 + F k 2 _p){k 2 +F + k 2 _ F ) coor-N 

a p + b p (235) 

n,p h,p 

{a + b) -(k 2 +F -k 2 _p) (236) 



k 



2 



using 



kj_p -\- k_p — kp. (237) 



Note how d = 2 was essential to this argument: in d = 3 we would have k^_ F + k^ F = kp. 

Thus the k A terms in £(k±) are not the cause of the S 4 term. However, a k 6 term in £ (fc±), can be shown to 
produce an S 4 term in E(S). 

Thus the apparent free-field behavior is tied to the smallness of terms of order k e and higher. To understand why 
the k e term is so small, we turn to Eqn. (|154|l for H . Expanding the sine 2 in a series, we find the k e term is down 
by a factor of at least 15 (50) relative to the k 2 term, at A = (A = 1), all the way up to k = kp. Presumably this 
feature (and its counterpart in the gapped case) persists in the HF approximation to H and keeps E(S) essentially 
quadratic, which in turn mimics free-field behavior. 

The reader is referred the original work (Shankar, 2001) for a proof that E(S) = E(0) + ^S 2 implies that A(r) 
will be r-independent in the gapped case as well. 

Composite Fermions are not free fcrmions but are like Landau quasiparticles in a Fermi liquid . These objects 
too are labeled by free-particle quantum numbers and long-lived. They do have fairly strong interactions: the 
dimensionless Landau parameters that describe these interactions are not small and produce effects like zero sound. 
They are adiabatically connected to free fermions in zero field just as the CF's are adiabatically connected to free 
fermions in a reduced field B*. 



D. Effective potentials for experimental systems with disorder 

In comparing to experimental results one cannot neglect disorder. Our theory ignores disorder and our results are 
completely determined by the electron-electron interaction. Here we ask if it is possible that a ZDS potential with 
some effective A can describe a dirty system. First of all, we realize this cannot be true with respect to all observables, 
if at all it is true for any. For example, if one were considering conductance, one knows the electron in a disordered 
potential will typically get localized whereas no ZDS interaction will predict this. As for transport gaps, the present 
day samples, with a disorder broadening of the same order as the gaps, again preclude this possibility. Magnetic 
transitions, on the other hand, are controlled by total energies and one may expect that disorder will have a rather 
innocuous effect and can be represented in an average way by some translationally invariant interaction. We raise 
this issue because in several magnetic phenomena to be described shortly, it appears that a single A characterizes a 
sample. Specifically, A extracted from one data point can be used to explain the rest of the data from that sample. If 
the other data points differ only in the temperature T, the same A is used. If it differs in B or n or v, the following 
scaling argument applies (Ando, Fowler, and Stern, 1982): In a heterojunction, the donors of density n produce a 
confining linear potential of slope that goes as n. If one considers a variational wavefunction of the Fang-Howard 
(1966) form ip(z) = A(w)zexp(~z/w) in the transverse direction, then the optimal w (to which A, the well width 
must be proportional), varies as w ~ n~~3. Consequently the dimensionless width, A = A/ 1 varies as 

A ~ n~%B* ~ B 1/6 ^^ ~ n 1/6 u-i. (238) 

Arguments can be given (Shankar 2001) for why in certain limiting conditions, not realized in today's experi- 
ments, an effective potential exists. The question of why this works in realistic situations far from this limit remains 
unanswered. 

With this preamble, let us turn comparison with experimental results. Kukushkin, von Klitzing, and Eberl (1999) 
vary both n and B and drive the system through various transitions at T — (by extrapolation). The field B is 



This was already suggested by Halperin, Lee, and Read (1993). Also in this context, see the work of Mandal and Jain (2001a). 
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TABLE V Critical fields based on a fit at f .The rows are ordered by the last column which measures density. 
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always perpendicular to the sample. We will compare the hamiltonian theory to these experiments by calculating the 
critical fields at which the v = | and v = i systems saturate (P — 1) and the gapped fractions undergo transitions 
from one quantized value of P to the next. 

Let us recall that as far as these transitions go, the systems behave like free fcmions of mass m p which is independent 
of the index r which labels how many LL's have reversed their spins in the gapped case or the size of the up and down 
Fermi circles in the gapless cases. 

We consider B c, s at which the systems at 1/4, |, |, |, and \ lose full polarization (r = for gapped cases, saturation 
for the gapless cases) and, for |, also the r = 1 transition, |3, 1) — |2, 2). 

An experimental complication needs to be addressed first. Each of these transitions seems to take place via a 
narrow intermediate step (Kukushkin, von Klitzing, and Ebcrl, 1999) with a polarization half-way between the ones 
allowed by CF theory based on spatially homogeneous states. We use the center of these narrow steps as the transition 
points for comparison to the present theory. The physics of these intermediate steps will be addressed in Section Ixl 

In accordance with our strategy we find A by fitting the theory to the experiment for a particular transition. We 
will use this value (or for samples with changing field and density, Eqn. (J23EJ0 to predict B c for other transitions 
using Eqn. (|25U|) . We obtain As = 1.42 from the transition |3, 0) — |2, 1) at B c = 4.5T. 

For the gapless cases, there are two equivalent approaches. First, at the critical field the Fermi energy of the up 
spins equals the Zeeman energy of the down spins: 



eB c 



2m e 



kp 2im eBv 



2m E 



(239) 



Equivalently we can write for the total ground state energy density E Z (S), (where the superscript indicates that 
the Zeeman energy is included), 

E z (S) = ^S 2 -g n e g± f (240) 
K ! 2 y 2m e cos ' 

where a = Tr/m p . This expression is minimized (forP < 1) to give P. Setting P = 1 gives the critical fields. 

The comparison to experiment is made in Table (|Vl Note that in rows above (below) |, where we fit A, the 
predicted B c, s are lower (higher) than the observed values, i.e., the actual A's are less (more) than what Eqn. I|238|) 
gives. This is consistent with the expectation that interactions will increase the effective thickness with increased 
density. If we fit to the | point, we obtain similar numbers, with the agreement worsening as we move off in density 
from |. Thus | was chosen as the fitting point since its density was somewhere in the middle of all the densities 
considered. 



IX. PHYSICS AT NONZERO TEMPERATURES T > 0. 

So far we have seen that the EHT may be used to compute quantities such as gaps, particle-hole profiles, critical 
fields for magnetic transitions and so on to 10-20% accuracy. All such quantities have been readily computed using 
trial wavefunctions, giving numbers that are superior to ours. Our main emphasis has been to expose the underlying 
physics as transparently as possible and to resolve questions such as why CF's behave like free particles on some 
occasions. 

We turn to physics at finite T (Murthy, 2000c, Shankar, 2000, 2001) where the hamiltonian method has few rivals. 
Exact diagonalization (Chakraborty and Pietilalnen, 1996, Chakraborty, Niemela, and Pictilaincn, 1998) is limited 
to very small systems and trial wavefunctions typically cover the ground state and very low-energy excitations. The 
hamiltonian approach is able to yield the temperature dependence of polarization P and the relaxation rate 1 /T± for 
the gapless states in the thermodynamic limit. If A is treated as before (fit to one data point per sample) we will see 
it is possible to give a very satisfactory account of experiments in gapless systems up to about l°K, which is of the 
order of the Fermi energy (Shankar, 2000, 2001). 
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FIG. 9 Hartree Fock energies at v = i for up and down spins (upper and lower curves) at T — .3°K at B — 5.52 and zero 
tilt. Note that they are not simply quadratic in momenta, and that at the chemical potential, indicated by the horizontal line, 
the two graphs have very different slopes, i.e., density of states. 



We then address finite-temperature polarization in gapped states, which is complicated by a nonzero spontaneous 
polarization and the attendant spin waves. It turns out to be essential to take the finitc-T behavior of these spin- 
waves into account. Once this is done, the theoretical predictions are in excellent agreement with the experiment up 
to several Kelvin (Murthy, 2000c). 

We start with the gapless case since it is simpler. 



A. Polarization and relaxation in gapless states 

The polarization P is computed as follows. First we compute the HF energy of a particle including the Zeeman 
energy, which is the self-consistent solution to 
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where the superscript on £± reminds us it is the total energy including the Zeeman part, the Fermi functions 



;a*i) 



exp [(£f(fc) - n)/kT] +1 



(241) 



depend on the energies £±{k) and the chemical potential [i. Fig. 9 shows some typical results. At each T, one must 
choose a /j,, solve for £±(k) till a self-consistent answer with the right total particle density n is obtained. From this 
one may obtain the polarization by taking the difference of up and down densities. As usual we use the ZDS potential 
for which 



v(q) 



2 li 



(242) 



The computation of the longitudinal nuclear relaxation rate 1/Ti is more involved (Shankar, 2000, 2001). The 
question we ask is the following. The fermions are in a quantum well, with their density varying across the width. 
So the nuclear relaxation rate will be a function of position. Consider a nucleus at the center of the quantum well, 
(as well as the x — y plane) where the density is the largest. Let us call this point the origin and let 1/Ti be the 
relaxation rate here. The theory predicts 
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FIG. 10 Comparison to the work of Dementyev et al. The value of A is fit to P at 300 mK, B± — 5.52 T and the rest follows 
from the theory. Notice the correlation between the curvature of 1/Ti and the limit of P as T — > 0°K. 



where k± are solutions to £±(k±) = E, Iq is the Bessel function, Eq is the lowest possible energy for up spin fermions, 
and K™ ax is the measured maximum Knight shift (at the center of the sample) for the fraction v = | or |. 

Here is a rough description of the derivation, the details of which may be found in Shankar (2000). Suppose for 
a moment we were dealing with electrons and not CF's. The Knight shift at the chosen point, the origin, will be 
determined by the spin density there. The same parameter enters the 1/Ti calculation quadratically. This is why 
K™ ax enters the answer. Thus K™ ax is not calculated ab initio but taken from the same experiment. The density of 
states and Fermi factor are standard. The only new feature here is the presence of F(h+, fc_) which reflects the fact 
that the spin density has to be projected into the LLL when going to the CF basis. The effect of this factor (which is 
none other than the e~ q 1 I 2 which appeared on the projected charge density) is to suppress processes with momenta 
much larger than l/l, as these have no place within the LLL. 

We now compare to some experiments at v = \ and T > 0. Consider first Dementyev et al (1999). From their 
data point P = .75 for B = B ± = 5.52T at 300 mK We deduce 

A = 1.75. (246) 

We have once again chosen to see to what extent a sole parameter A, can describe P and 1/Ti for the given sample 
at a given B±, but various temperatures and tilts. 

Since there does not exist a model, including disorder, that describes how A should vary with tilt we include no 
such variation. 

Dementyev et al (1999) find K™ ax = 4.856 • 10~ 7 °K, which is believed to describe a saturated system at P = 1. 

3 

They estimate that K™ ax = . 953TsT7 ltt:r , which is what we need here. Given this information, 1/Ti follows. 

2 3 

The top and bottom halves of Figure 10 compare the HF calculation of 1/Ti and P respectively, to the data. (The 
graphs for 1/Ti differ slightly from those in Shankar (2000) since the present calculation treats the spin of the CF 
more carefully. The 1/Ti graph at 5.52 T appears a little jagged since it was computed at just six points which were 
then connected. This is not apparent in the tilted case since the points lie on a straight line. 

Dementyev et al (1999) had pointed out that a two parameter fit (using a mass m and interaction J), led to 
disjoint sets of values (in the m — J plane) for these four curves. Given that H is neither free nor of the standard 
form (p 2 /2m + V(x)) this is to be expected. By contrast, a single A is able to describe the data here rather well since 
H has the right functional form. Given how the theory fits the polarization data up to the Fermi energy of ~ 1°K, 
it is clear that changing the data point used to fix A will be inconsequential. 

The present work establishes a phenomenological, nontrivial, and nonobvious fact that a single A parameter, 
determined from one data point, can describe both P and 1/Ti for the given sample under a variety of conditions. 
That the fitted A is larger than the Local Density Approximation value makes sense, as both disorder and LL mixing 
will lower the gap and raise A. 
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FIG. 11 Comparison to Melinte et al at = and 60° with A = 1.6 fit to P at 60 mK, B± = 7.1 T. 



Consider next sample M280 of Melinte et al (2000) which had P = .76 at .06° AT and B = B± = 7.1T, from which 
we deduced A = 1.6. Figure 11 compares the theoretical T -dependence of P(T) with data. There is a factor of 2 
between theory and experiment for 1/Ti, not shown. 

The fair agreement for the tilted cases is unexpected in both the Dementyev et al (1999) and Melinte et al (2000) 
data . First of all, orbital effects have to be considered due to the tilt. The thickness parameter A can be affected by 
it. As pointed out by Jungwirth in a private communication, once there is an in-plane component of B, the problem 
is no longer rotationally invariant. This means that our states are no longer HF states and get scattered into each 
other by the potential. No attempt is made here to take into account all the effects of the tilt. Instead we included 
just the increased Zeeman coupling and hoped for the best. 

For the benefit of others who measure 1/Ti at v — | in the future on similar samples, we give some very approximate 
formulae (to be used for zero or small tilts). From Figure 10, we note that in general, the graphs of 1/Ti become 
linear and parallel for temperatures above .3 °K. In this region we can write 



dT 



K 



lO-'V 1 TAT 1 



for T > .3°K 



(247) 



with K the Knight shift in KHz and 



10 10 /cm 2 



(248) 



(In this approximate formula, we ignore the A dependence of Eqns. (|243I245[1 . and the distinction between the 
average and maximum Knight shift.) 

The graphs do not generally obey the Korringa-like law because as T — > they are sublinear (supcrlinear)for 
saturated (unsaturated) cases. Only the critical case with P(0) — > 1 as T — > is linear. For T > .3 ° K, (which in 
general must be replaced by either the energy gap or energy overlap between the up and down Fermi energies) we 
have the approximate result 



i ~ [3 T(°K) + C] \§ 
J i \_n 

C = (critical) 

= > (unsaturated) 

= < (saturated) 



10 



(249) 

(250) 
(251) 
(252) 



For the critical case (only) , we have a Korringa law 



Ti T(°K) 



[3 



l 2 



(253) 



For Dementyev et al (1999) C ~ 1. This value may be used as a first approximation. For more accurate results 
one must solve Eqns. (|243I245|I . 
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B. Polarization in gapped states 

Let us now turn to the theoretical description of the finite-temperature polarization of the v = ^ state, detailed 
measurements of which have recently been carried out (Khandelwal et al, 1998, Melinte et al, 2000), Freytag (2001). 
Theoretical details can be found in Murthy (2000c). 

While the CF-HF approximation provided an adequate description of v = ^ it is qualitatively incorrect for a 
spontaneously polarized state like v — \. The reason is that it underestimates the effects of the excitations that 
destroy order. For example, it predicts a nonzero spontaneous polarization at E z = gpB tot = and T > 0, in 
violation of the Hohenberg-Mermin- Wagner theorem (Mermin and Wagner, 1966, Hohcnbcrg, 1967) that forbids 
spontaneous breaking of a continuous symmetry in two dimensions, except at T = 0. The T = spontaneous 
polarization is driven by the fact that fermions of the same spin will avoid each other due to the Pauli principle 
and thus have a lower interaction energy, just as in the v = 1 spontaneous quantum ferromagnet 41 , which has been 
extensively studied theoretically 42 and experimentally 43 . 

Thus we must begin by understanding the disordering mechanism and seeking a proper description of it. Consider 
a fully polarized v = i state at T = with all the fermions in the CF-LLL state. As the system is heated, some 
fermions will go to the n = 1 CF-LL with spin up (which does not change the polarization) and some will go to spin 
down n = CF-LL, which reduces the polarization and costs an energy per spin-flip of A$r, the spin-reversed gap. 

This description completely misses the spin- waves which are related to the particle-hole excitations as follows. Just 
as in the case of the magnetoexciton (Figure 4) where no spin was flipped, a spin-flip particle-hole excitation has an 
energy that varies with q 44 . As q — > oo, the dispersion settles down at A$r, the energy to create a widely separated 
pair. However as q — > the pair energy vanishes, for Ez = 0, by Goldstone's theorem. For nonzero Ez the long 
wavelength limit for the spin-wave energy is Ez by Larmor's theorem. These are the modes to reckon with for they 
are very low in energy and plentiful at low temperatures. 

The simplest way to describe these spin-waves (including their self-interaction) is the continuum quantum ferro- 
magnet (CQFM)(Read and Sachdcv, 1996). This model assumes all high-energy modes (at the electronic and CF 
cyclotronic scales) have been integrated out. The only modes left are slow and long- wavelength fluctuations of the 
spin polarization, which have the action 

l/T 

S = J d d x J dr{iM A(n) • V r n + f (V x n) 2 

o 

-M H • n + • • •) (254) 

where Mo is the magnetization density, n is a local vector of unit length pointing in the direction of the magnetization, 
A(n) is the field that implements the Berry's phase needed to obtain the correct quantum commutation relations 
between the spin components, p s is the spin stiffness, and H = <7*/zbB is the Zeeman field (|H| = Ez)- 

This model is still nontrivial (because of the condition |n| 2 = 1). However it can be solved (Read and Sachdev, 
1996) in the limit when N — > oo, where N is the number of components of n. The limit N — > oo appears to describe 
the actual case of N = 3 in the case v = 1. 

Our strategy then will be deduce reasonable values for Mq, the magnetization density, and p s , the spin stiffness, 
and then plug them into the known results from the large-N limit for the magnetization as a function of T is given by 

P(T) =M^ M {r,h) (255) 

where r = p s /T and h = Ez /T are scaling variables, and $m is a known scaling function (Read and Sachdev, 1996) 45 . 

To find the values of the parameters M and p s corresponding to the v = ^ state, we will use some Hartree-Fock 
results. Since the underlying fermionic theory responds to temperature by self-consistently modifying occupations 
and energies, we expect to obtain temperature-dependent values 46 M (T) and p s (T). 



The classic study of nontrivial interaction effects (beyond ferromagnctism) at v = 1 is Sondhi, Karlhede, Kivelson, and Rczayi (1993). 
Kopictz and Chakravarty (1989), Read and Sachdev (1996), Kasner and MacDonald (1996), Haussmann (1997), Timm, Girvin, Hcnelius, 
and Sandvik (1998), Kasner, Palacios, and MacDonald (2000), 

Tycko, Barrett, Dabbagh, Pfeiffer, and West, (1995), Barrett, Dabbagh, Pfeiffer, West, and Tycko, (1995), Aifer, Goldberg, and Broido 
(1996), Manfra, Aifer, Goldberg, Broido, Pfeiffer, and West (1996). 
For an example of this calculation see Murthy (1999) 

There are actually two different large- ./V approximations corresponding to the fact that the symmetry group can be viewed as an example 
of O(N) with N = 3, or as an example of SU(N) with N = 2. Both are considered by Read and Sachdev (1996) and scaling functions 
given. 

This should be contrasted with the case v = 1 where due to the huge exchange gap, particlc-holc excitations are frozen at all temperatures 
of interest, and the parameters of the CQFM are T-independent. 
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FIG. 12 Polarization P(T) as a lunction ol T for | for the data of Khandelwal et aVs (1998) 10W sample. The thickness 
parameter has been set to A = 1.5. Note the important role played by the spin-waves in bringing the HF value of the polarization 
down to the nearly experimental value. 

First consider the spin stiffness p s (T). At a given temperature T the self-consistent occupations Nf,gs{&,ti) and 
energies e(<7, n) in the ground state are computed using the procedure described in the gapless case. Now one creates 
a twisted spin state and computes the HF energy of the twisted ground state (Murthy, 2000c) , and thence the excess 
energy to order q 2 . Comparing to the energy cost of a twist in the CQFM, which is (p s /2)L 2 q 2 , one finds the spin 
stiffness 

P' = ifc / t€f W ( S ) E \Pn in2 {s)\ 2 X 

Tlx ,™2 

(N F (\, m) - N F (i, m))(JV>(T, n 2 ) - N F (l, n 2 )) (256) 

where L 2 is the area of the system, and p ni n 2 is the matrix element of equation (jDlfl . The above should be regarded 
as an estimate for the twist rather than a rigorous calculation (even in HF), since ideally one should compute the free 
energy cost of a twist, rather than just the internal energy cost, as we have done. 

To find Mo(T) we need a more devious approach (Murthy, 2000c), based on Eqn. 1255|) . We already know how to 
compute the CF-HF magnetization Phf(T). We set 

Phf{T) = M (T)$ M (r = 0,h= A SR /T) (257) 

and justify it as follows. In the CF-HF theory the particles and holes are treated as independent, or noninteracting, 
with a gap equal to Asr independent of the distance between them. This corresponds to a collective mode dispersion 
that is completely fiat, w(q) = Asr. The CQFM description that corresponds most closely to the CF-HF is the one 
that has the same spin- flip excitation spectrum, namely one with no spin stiffness (r = 0), and an effective Zeeman 
field E e J f = A S r. 

Armed with Mq(T) and p s (T), we evaluate the large- N scaling functions and hence P{T). The predictions are 
shown in Figure 12 for parameters corresponding to the experiment of Khandelwal et al (1998). As can be seen, CF-HF 
overestimates the spin magnetization considerably at low temperatures, but the inclusion of interacting spin-waves 
gives a prediction in almost perfect agreement with the data. The results are quite insensitive to A and a typical value 
of A = 1.5 was used. The results are also insensitive to which of the two large- N approximations is used. Similar 
good agreement over a wide range of temperatures is found in the comparison to the data of Mclintc et al's (2000) 
M242 sample, shown in Figure 13. Here we have plotted the Knight shift vs. the temperature. The Knight shift at 
very low temperature shows considerable scatter, and an appropriate intermediate value has been used to fit to the 
theory. 

It is somewhat surprising that the theory agrees so well with the data up to temperatures much higher than those 
for which agreement was found for the gapless fractions. We expect the theory to work only when CF's are well- 
defined, which should be true up to temperatures of the order of the ^ gap. In fact, the agreement persists to about 
10K, which is higher than a typical activation gap for | in such samples. For more details see Murthy (2000c). 
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FIG. 13 Polarization P(T) as a function of T for | for the data of Melinte et al's (2000) M242 sample. The thickness parameter 
has again been set to A = 1.5. Once again spin fluctuations are seen to be crucial in reproducing the experimental polarization. 
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The same considerations can be applied to the v = | state. Since it is unpolarizcd for E z = we can stop at the 
CF-HF stage. Figure 14 shows the P(T) curves for v = | for A = 1.51 for a range of Zccman couplings. The rise 
and fall of P(T) at small Ez can be understood as follows. At T = 0, the system is in a singlet state with the up 
and down spins occupying the CF-LLL. As the system is heated, some particles go the next LL. Here the ones with 
spin parallel to the external field are preferred due to the Zccman term and hence a nonzero polarization develops. 
At higher T, the two spin species start to get occupied with nearly equal probabilities and P starts to decline. There 
is a transition to the fully polarized state around Ez = O.OlEc- 

One may compute P(T) for arbitrary fractions by exact diagonalization (keeping all the excited states) and sub- 
sequent calculation of thermodynamic quantities (Chakraborty and Pietilaincn, 1996). Due to computational limita- 
tions, this method is restricted to fairly small systems. For example, the largest system studied by Chakraborty and 
Pietilaincn (1996) for v = | has 5 electrons, and for v = | has 4 electrons. Allowing for this, the results seem fairly 
consistent with ours. Currently experimental data are not available for comparison. 

Finally one use the CS theory to do finite T calculations, using to* as a free parameter instead of A. 



X. IN HOMOGENEOUS STATES OF CF'S 



We have repeatedly emphasized that CF's are interacting particles. The uniform liquid states studied so far were 
not contingent on CF interactions- such states could exist even for noninteracting CF's. We now turn to some 
interesting inhomogeneous states, which would not exist if the CF's were free, such as the high-field Wigncr Crystal, 
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and possibly the partial polarized states seen in the gapped fractions by Kukushkin, von Klitzing, and Eberl (1999). 

A. The high-field Wigner Crystal 

The clearest instance of an inhomogeneous state in the fractional quantum Hall regime is the high-field Wigner 
crystal (WC) (Yoshioka and Fukuyama, 1979). If we start with a free electron gas and slowly crank up the interactions, 
or lower the density, the emphasis goes from the kinetic to potential energy. In the Wigner crystal the electrons seek 
an arrangement aimed at minimizing the potential energy. In the Hall case, due the quenching of kinetic energy, the 
WC seems even more likely at low densities. It was initially proposed (Fukuyama and Platzman, 1982) as a possible 
explanation of the fractional quantum Hall effect. However, it was soon realized (Yoshioka and Lee, 1982) that this 
state is very different from the fractional quantum Hall states, because both its longitudinal and Hall conductances 
vanish at zero temperature and when the crystal is pinned by disorder (Yoshioka 1983, MacDonald 1984), leading to 
an insulating state. Further, the WC state is not tied to any particular commensurate filling. 

There are now experimental observations 47 which support the existence of the WC near v = i. In fact, experiments 
see a re-entrant transition (Jiang et al, 1990, Jiang et al, 1991) in which there is a putative WC both above and below 
i. To see how this might happen, recall that the incompressibility of the Laughlin state results in a downward cusp 
in the ground state energy as a function of filling factor in the neighborhood of (say) |. Thus, one can easily imagine 
that in a small neighborhood of -| the Laughlin liquid has a lower energy than the WC, leading to the re-entrant WC 
near |. 

Theoretical work on the WC 48 , based on a study of trial wavefunctions and collective excitations has established 
that the Laughlin state becomes unstable around v « ^ to a density wave even in the absence of disorder. In the 
wavefunction approach, Hartree-Fock and (weakly) correlated wavefunctions have also been written down, and their 
energy evaluated. By studying the excitonic instabilities of the Laughlin liquid, Jain and Kamilla (1998) showed that 
it becomes unstable to crystallization around v = |. 

Let us begin with one of the simplest HF wavefunctions for the crystal, that of Maki and Zotos (1983): 

* HF ({ri})=Al[(t>nM)- (258) 

i 

where A is the antisymmetrization operator, and 0R i is a single-particle wavefunction that is localized at R; (lattice 
site) and belongs to the LLL. It is given by 

0R s (r) = e -l r - R *l 2 / 4 'o-* rxR ^/ 2i o. (259) 

The wavefunction l|258|l can be improved by adding a correlation factor corresponding to including small fluctuations 
around the HF state. The energy of this state becomes lower than that of the liquid state at about the experimentally 
right filling fraction (y « i) (Maki and Zotos, 1983, Lam and Girvin, 1984, Levesque, Weiss, and MacDonald, 1984). 

Experiments (Jiang et al, 1990, Jiang et al, 1991) show transport gaps two orders of magnitude smaller than the 
theoretical estimate as calculated using the Hartree-Fock approximation (Fukuyama, Platzman, and Anderson, 1979, 
Yoshioka and Fukuyama, 1979). The measurements of the Hall resistivity p xy are surprising as well (Goldman, Wang, 
Su, and Shayegan, 1993, Goldys et al, 1992) . As mentioned above, the WC is expected to have a vanishing Hall 
conductance a xy = (when pinned), which implies a vanishing Hall resistance p xy — 0. However, the experiments see 
"Hall insulating" behavior (Zhang, Kivelson, and Lee 1992), that is, p xy « These problems led Yi and Fertig 

(1998) to consider crystalline states with correlation zeroes that keep electrons apart. Each electron is combined with 
2s vortices to obtain the trial wavefunction 

*(w)=^n^-^) 2s n^(^)- ( 26 °) 

i^j i 

They also find that near 1/5 the best energies are obtained by attaching 4 zeroes to each electron. In other words, 
they take the Laughlin FQHE wavefunction at ^ which is a product of a Jastrow factor with quartic zeros and a single 
filled CF Landau level xi, an d replace xi with a crystal. The Coulomb energy for this wavefunction is computed using 



Jiang, Willett, Stormer, Tsui, Pfeiffer and West (1990), Jiang, Stormer, Tsui, Pfeiffer and West (1991), Engel, Sajoto, Li, Tsui and 
Shayegan (1992), Goldman, Wang, Su, and Shayegan (1993), Goldys, Brown, Dunford, Davies, Newbury, Clark, Simmonds, Harris and 
Foxon (1992). For a recent review of the experimental situation, see Shayegan (1997) 

Maki and Zotos (1983), Lam and Girvin (1984), Levesque, Weiss and MacDonald (1984), Price, Platzman, and He (1993). For a recent 
review of the theory, see Fertig (1997). 
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Monte Carlo methods. Yi and Fertig (1998) have shown that the ground state energy of the correlated WC is lower 
than that of the usual WC at experimentally relevant filling fractions. Moreover by introducing Laughlin-Jastrow 
correlations between the interstitials and the lattice electrons the experimentally observed p xy (Zheng and Fertig, 
1994a, 1994b) can be explained. Unfortunately, the method becomes too computationally demanding to allow one to 
calculate other quantities of interest, such as the excitation spectrum. 

This is a situation tailor-made for the EHT. Before we describe the theory and its results, let us note that there 
are two features of the experiment which hint that a CF-WC is involved. The first is the nonmonotonic behavior of 
the gap near | , which shows that the | Laughlin liquid correlations are being felt in the nearby WC state. Secondly, 
the threshold electric field beyond which the WC becomes depinned and starts sliding increases near i. Below we 
will see its connection with the structure and properties of the CF-WC. 

Now let us see how to set up the EHT. Since attaching the zeros to electrons converts them into CF's, we are 
naturally led to consider a WC of CF's. Here is how one proceeds (details can be found in Narevich, Murthy, and 
Fertig, 2001). The crystal is characterized by density wave order parameters at the reciprocal lattice vectors G. 

Aw(G) = ^P-J2 e ~ iGmX < d lx-G y i^/2 d n',x + G^/2 > (261) 

X 

First, one has to assume a particular lattice structure (shape and size). An important parameter of lattice is the 
number of quanta of effective flux that penetrate each unit cell 49 . Things become simple when this number is rational, 
of the form p/ q, where p and q have no common factors. In this case each CF-LL breaks up into p sub-bands with the 
total number of states in the original CF-LL being equally divided among the sub-bands (Yoshioka and Fukuyama, 
1979). While the original CF-LL had a sharp energy, the sub-bands have a nonzero energy dispersion, which can be 
found from the HF hamiltonian. One finally closes the circle by demanding self-consistency; the ground state formed 
by filling up the sub-bands with the correct number of particles should reproduce the assumed form of Eqn. (|261|l . 
Once one has a self-consistent HF solution, the transport gap is found as the energy difference between the centers 
of the highest occupied sub-band and the lowest unoccupied sub-band. The WC is characterized by one particle per 
unit cell, and depending on the filling, this translates to different p and q. 

Here are the results: The gaps calculated from the Hamiltonian theory are within a factor of two of the experimental 
results, (Figure 15) which is to be contrasted to the more than two orders of magnitude disagreement with electronic 
HF calculations. The theory also predicts that the gap below 1/5 should be substantially more than the gap between 
1/5 and 2/9 (the next quantum Hall liquid state). Narevich, Murthy, and Fertig (2001) are not able to go very close 
to 1/5, because then p and q become large, and the problem becomes computationally prohibitive. 

Two other results emerge at a qualitative level. First, the shear modulus becomes very small as | is approached, 
which means that the crystal becomes very soft to deformations. The reason is that as | approaches the WC tends 
more and more to a Laughlin liquid. The density wave order parameter decreases, and the shape of the crystal 
matters less and less. This is connected to the depinning threshold of the WC to electric fields. This connection can 
be understood as follows (Fukuyama and Lee, 1978, Blatter et al, 1994); when the crystal is stiff it is not able to 
take advantage of all the minima in the disorder potential, and so is not strongly pinned. However, when the crystal 
becomes soft, it does not cost a lot of energy to deform and take advantage of local minima of the disorder, and the 
crystal becomes more strongly pinned. Disorder may also have a dominant role in the behavior of the gap near jr. 
Since the crystal is soft, its actual shape is determined by the local disorder. Different shapes give gaps varying by 
about a factor of two. 

The second qualitative result is that the density inhomogeniety in the CF-WC state is quite small in absolute terms 
(about 20% of the background density) for all ^, not just near |. This is at odds with the conventional view that 
electrons are localized in a WC. The correlation zeroes seem to prefer a more homogeneous state 50 . The electronic 
WC has to partially melt and become more homogeneous to accommodate Laughlin-Jastrow correlations. This has 
interesting similarities with earlier ideas concerning cooperative ring-exchanges in a WC, and the melting of the WC 
(Kivelson, Kallin, Arovas, and Schrieffer, 1986a, b). 

B. Partially polarized Hall Crystal states 

Another class of candidates for an inhomogeneous state of CF's stems from the observation of Kukushkin, von 
Klitzing and Eberl (1999), who measured the spin polarization of various principal fractions as a function of Ez, the 



Since the CF's see effective flux rather than external flux this point is crucial. 

It is easy to see that correlation zeroes become meaningless if electrons are strictly localized; then they do not come close enough for 
the zeros to be operative. 
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FIG. 15 Experimental (stars, from Jiang et al 1991) and theoretical (squares and diamonds, from Narevich, Murthy, and Fertig 
2001) gaps versus v for the Wigner crystal. The squares are theoretical gaps for the triangular lattice, while the diamonds 
represent gaps for an oblique lattice. A change of shape of the lattice just above -| could explain the nonmonotonicity of the 
experimental gap. 



Zeeman energy and found polarizations not expected in the CF theory with uniform states. As an example, consider 
v = |. There are two filled CF-LL's, and at small Ez we expect that the state is unpolarized with the n = 0, | 
and n = 0,1 CF-LL's occupied and the rest empty. At large Ez we expect a fully polarized state with n = 0, f and 
n = 1, t occupied. If these arc the only states we are allowed to think about, then at some intermediate Ez there 
should be a first order transition. Such transitions between P = and P = 1 were extensively discussed earlier. 

However, Kukushkin, von Klitzing, and Eberl (1999) observe a plateau at half the maximum polarization in a narrow 
region of Ez for |. Analogous plateaus are observed at other principal fractions 51 . Other puzzling phenomena near 
this transition are that of hysteresis over a range of fillings near 2/5 (Cho et al 1998) and very slow dynamics (Eom et 
al 2000). One of us (Murthy 2000a) has proposed that these observations might be explained by considering partially 
polarized Hall Crystals of CFs. 52 

Let us first review some terminology. A state which shows the co-existence of density wave and quantum Hall order 
(a quantized <J xy ) is called a Hall Crystal. Such states were implicit in the ring-exchange theory (Kivelson, Kallin, 
Arovas, and Schrieffer, 1986a, b), but were first postulated explicitly by Halperin, Tesanovic, and Axel (1986), and 
analyzed in great detail with explicit examples by Tesanovic, Axel, and Halperin (1989) 53 . A WC is an example of a 
Hall Crystal, as is the lattice of quasiparticles formed at a filling v — | + S, provided the density wave is pinned. The 
reason the partially polarized state at 2/5 has to be a Hall Crystal (if it is a crystal at all) is that no changes in the 
quantized Hall conductance are observed as Ez is varied. 

Now let us understand the nature of the proposed state (Murthy, 2000a). Imagine that CF's were really free. Then 
as Ez is raised from the system remains in a singlet state until the n = 0, J. CF-LL (whose energy is increasing) 
crosses the n = l,f CF-LL (whose energy is decreasing) at a special Zeeman energy E* z . This is when the system 
makes the first-order transition. Now let us "turn on" interactions between CF's, which we know are present in the 
full theory. At E* z there are twice the number of degenerate states (from the n = 0, J. and n = 1, t CF-LL's) as there 
are CF's at the Fermi level. The true many-body ground state will therefore be picked by the interactions. This state 
should be stable in a small range of Ez near E* z determined by the strength of the CF interactions. One possibility 
is a charge/spin density- wave state. 

As we saw in the previous subsection, when the system self-consistently chooses a particular density-wave order 
with p/q quanta of effective flux per unit cell, each CF-LL splits up into p subbands of equal degeneracy (Yoshioka 
and Fukuyama, 1979). The partially polarized | state has two quanta of effective flux penetrating every unit cell, 
which means p = 2, q = 1, and each CF-LL is split into two subbands. The proposed state corresponds to occupying 
3 subbands of |-CFs and 1 subband of |-CFs. Once again a CF-HF calculation allows us to calculate the gap. The 



Very recently, polarizations not allowed by translationally invariant CF states have been seen at ^ by NMR techniques as well: Freytag, 
Tokunaga, Horvatic, Berthier, Shayegan, and Levy (2001). 
An alternate explanation is provided by Mariani et al (2002). 

The topological numbers associated with these states had been noticed before by a number of authors: Wannier (1978), Thouless, 
Kohmoto, Nightingale, and den Nijs (1981), Johnson and Moser (1982), Thouless (1983), MacDonald (1983), Dana, Avron, and Zak 
(1985), Avron and Seiler (1985), Avron and Yaffe (1986), Kunz (1986). However, Tesanovic, Axel, and Halperin (1989) were the first to 
give a concrete example where nontrivial topological integers were realized, and to make an important connection between the topological 
integers and gapless collective excitations. For an extension of the Hall Crystal concept to the FQH regime, see Kol and Read (1993). 
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theoretical prediction for the range in Ez over which the state should be stable is consistent with the experiment 
(Murthy 2000a). 

This possibility can be investigated in more detail in the integer quantum Hall analogue of v — 2/5, which is v = 2. 
There it can be shown (Murthy 2000b) that there are indeed charge/spin density wave states whose energy (in HF) 
is less than both the singlet and fully polarized liquid states. It can also be shown that these states are Hall Crystals 
with a nonzero quantized Hall conductance (Murthy 2000b). 

However, there are other possibilities for the partially polarized state at v = 2/5. For example, Apalkov, 
Chakraborty, Pictilainen, and Niemela (2001), have proposed that the CF's belonging to the n — 0,1 and n = 1, j 
CF-LL's form a Halperin(lll) state (Halperin, 1983), which is a liquid. Their proposal is based on a finite-size calcu- 
lation in a spherical geometry where only the CF's belonging to the two "active" levels were kept, and their interaction 
was approximated as being a two-body interaction. While their state remains a contender, Apalkov, Chakraborty, 
Pictilainen, and Niemela, (2001) wrongly conclude, by neglecting some terms in the ground state energy (see Murthy 
2001c), that the ground state energy of the charge/spin density wave state would be higher than the Halperin(lll) 
analogue state. 

At the moment it is fair to say that the situation is not resolved. There are some key differences between the 
two proposals which could be used to distinguish them experimentally. For example, in the liquid state of Apalkov, 
Chakraborty, Pictilainen, and Niemela (2001) one would expect no gapless excitations whatsoever, and no spatial 
variation in the spin density, and thus the NMR Knight shift. On the other hand, there are generically gapless density 
excitations in the partially polarized Hall Crystal state (Tesanovic, Axel, and Halperin 1989, Murthy 2000b), and one 
would expect (at very low temperatures when the crystal is frozen) periodic spatial variation in the Knight shift, with 
the maximum Knight shift being the same as in the fully polarized state. A complete understanding of these partially 
polarized states will hinge on future experiments. 

XI. CRITICAL REVIEW OF EHT AND ITS APPROXIMATE SOLUTIONS 

The EHT was arrived at by starting with the electronic hamiltonian, introducing a new (pseudovortex) coordinate 
R„, combining it with the electronic guiding center R e to form the CF variables rj and R. A (weak) constraint 
x(q) — completed the picture. For v = 1/2 the formalism coincides with that of Pasquicr and Haldane (1998). 

As an exact restatement of the original electronic problem, the EHT is as good as any, with no loss of information, 
including that pertaining to higher LL's. Its main claim to superiority over other formulations is its amenability to 
Hartree-Fock and attendant approximations, because the problem is now expressed in terms of CF variables (rj and 
R) which have a natural HF ground state that fills an integral number of LL's. 

We discussed two approximation schemes, both of which rely on a nondegenerate HF state, but which differ in how 
the constraints are handled. Since the separation of the LLL physics (and thus the limit m — > 0) was straightforward 
in this formalism, we focused on this limit and the LLL. 

In a conserving approximation, one finds approximate Green's functions that respect the constraints. The time- 
dependent Hartree-Fock (TDHF) approximation is conserving. It was employed (Murthy 2001a) to show that in the 
LLL, S(q) ~ q 4 for gapped states, in compliance with Kohn's theorem. For v = 1/2 it was used by Read (1998) to show 
that the state is compressible, and the fermion effective mass has a logarithmic divergence at the Fermi surface, because 
the overdamped mode of HLR responsible for all this is generated by summing ladder diagrams. Thus there is no 
conflict between the HLR results based on CS fermions (that had unit charge) and those based on Composite Fermions 
with charge e* (which vanishes for v = 1/2.) However results that are easily obtained in the CS approach are obtained 
only following a careful, gauge-invariant calculation in the EHT. It may well be that without the CS results, the 
compressibility at v — 1/2 would have been overlooked in a description that succeeded in displaying dipolar expressions 
for the electron charge operator. While the primary use of the conserving approximation thus far has been to 
demonstrate that correct results in the infrared, long-wavelength limit require a proper implementation of constraints, 
it seems to numerically capture some of the short-distance physics as well, as evidenced by the magnetoexciton 
dispersions for | and | (Murthy, 2001a, section VIB). 

In the other approximation, (the shortcut) one uses the preferred charge p p (q) = p(q) — c 2 x(q) in place of p(q) 
on the grounds that in an exact calculation such a substitution is permitted. This single substitution manages to 
accomplish at tree level some of the work of the constraints: e — > e* , d — > cP, and S(q) ~ q 2 — > S(q) ~ q 4 . The 
constraints are then ignored, so that their effects at higher orders in ql are not built in, and may or may not be 
described well by the simple HF calculation that follows. Thus, although the EHT itself is exact, the small parameter 
ql is in the shortcut needed because of the way we handled the constraints. Note that it is ql and not ql* that has to 
be small, as evidenced by the success of the approach near and at v — 1/2 where I* — ► oo. 

Also missing in the shortcut is the overdamped mode and its attendant consequences. However in problems with 
a gap or T > this does not seem to matter. On this basis one expects it to work in problems where the potential 
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favors small ql, and the problem is still described by CF's 54 . 

These expectations were borne out in the computation of gaps, magnetic transitions and T-depcndcncc of polar- 
ization and nuclear spin relaxation, which were compared to results from exact diagonalization, trial wavefunctions 
and experiment. The very special form the hamiltonian that comes out of this theory does what ad hoc hamilto- 
nians parametrized by standard kinetic and interaction terms cannot: Explain a host of phenomena in any given 
sample (polarization, relaxation) with a single parameter describing electro-electron interactions. We are not aware 
of any other quantitative, analytical approximation scheme for computing these quantities which is not plagued by 
singularities as m — > 0. 

The preferred charge description also answers many qualitative issues: In what sense does an electron bind to 
its correlation hole (represented here by pseudovortices) to form the CF? How can CF's appear to be free in some 
magnetic phenomena when they are surely not free? What does the dipole picture really mean? How is an effective 
mass generated from the interactions alone? 

On the down side, the preferred charge as a way to implement constraint is peculiar to the FQHE and unprecedented. 
We do not know a priori how well this procedure will work and how to systematically improve it. We have no internal 
signal that the CF description is failing as the interaction is varied (say by sending A — > oo in the ZDS interaction) or 
by mapping v > 1 problems into effective LL problems with modified interactions, where Pfaffian states and striped 
states might be the answer. 

XII. SUMMARY AND OPEN QUESTIONS 

The goal of the hamiltonian theory was to start with the hamiltonian for two-dimensional interacting electrons in 
a magnetic field and arrive at a comprehensive description of fractional quantum Hall phenomena by a sequence of 
transformations to the final quasiparticles. The final theory was expected to encode and display the wisdom inherited 
from the study of trial wavefunctions, make the various pictures precise, and permit the computation of all physical 
quantities (to some accuracy) by the use of standard approximation techniques, such as HF. 

Thus we first surveyed the wavefunction approach to see what it could teach us and then described the two 
hamiltonian approaches: the Chern-Simons (CS) approach and our extension of it, the Extended Hamiltonian Theory 



The wavefunctions taught us the following. The FQHE ground states are described by incompressible fluids. The 
wavefunctions for the ground states and quaiparticle (or hole) excitations are independent of the mass m (in the 
m — > limit) and is built entirely of LLL wavefunctions. In the Laughlin fractions, the elementary excitation is a 
vortex of charge — l/(2s + 1) in electronic units. When an extra electron is introduced into the system, it is screened 
by 2s such vortices, leading to a composite fermion of charge e* = l/(2s + 1). These vortices sit on the electron. In 
the Jain fractions the wavefunctions are obtained by projection to the LLL. Prior to projection one can see the CF as 
made of an electron and 2s vortices sitting on them. After projection, many zeros are annihilated, and typically move 
off the electrons and are no longer organized into vortices. This holds for v = 1/2 as well, as is clear when f = 1/2 is 
seen as a limit of nearby Jain fractions. This means that the naive dipole picture based on zeros of the wavefunctions 
is not tenable. Remarkably, despite the dissolution of vortices on projection, the CF still has the charge of an electron 
and 2s vortices. The proper interpretation of this requires the hamiltonian approach, in particular the EHT, which 
in turn is an offshoot of CS theory, which was described next. 

We emphasized that our aim was not simply to get yet another exact reformulation of the problem, but one that 
lends itself to approximations. For example the problem in terms of electrons: 



contains all the answers in principle, but is flawed in practice because of the LLL degeneracy of H precludes the use 
of HF or perturbation theory. 

The degeneracy of the noninteracting problem is overcome by the flux attachment transformation of Chern-Simons 
theory. Lack of space prevented us from discussing the very first application of this idea to the Laughlin fractions 
by Zhang, Hansson, and Kivelson (1989). They converted the problem to one of composite bosons in zero average 
field, which led to a detailed analogy to superfluidity. This review focused on CF's, turning to the work of Lopez 
and Fradkin (1991,1992,1993) who implemented Jain's idea in operator form, by attached 2s flux quanta to electrons 



If one lets the "thickness" parameter of the Zhang-Das Sarma potential become very large, CF's will likely become unstable, and the 
state will not be dcscribable in the EHT. 



(EHT). 




(262) 
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converting them to CS fermions described by 

H cs = Y, {U+ ' 2 * CS:)2i +V (263) 

i 

U = p + eA* (A* = ^A_) (264) 
Vx : a cs : = ins : p : (265) 

At mean- field one now has a unique vacuum of p- filled LL's of CF's upon neglecting V and : a cs :. Their effects 
could be included perturbatively. In particular in an RPA approximation one could include both these and obtain 
the cyclotron mode at the right frequency, and in the Laughlin fraction, the wavefunction in the long-distance limit 
(Lopez and Fradkin, 1991,1992,1993). 

However, two problems remained: A singular limit as m, the bare mass of the electron, vanished, and a quasiparticle 
charge of unity (at tree level) instead of e* , owing to the fact that flux tubes and not vortices were attached. These 
features made it hard to make quantitative predictions. 

We mentioned Kalmeyer and Zhang (1992) who considered v = 1/2 and pointed out that impurity scattering would 
be stronger than naively expected since any charge inhomogeneity led to additional gauge flux (by the CS condition) 
which produced strong scattering. 

We then moved to the Halperin, Lee, and Read (1993) treatment of v ~ i, which was extensive and aimed at 
confronting theory with experiment at a quantitative level. The success of the CF theory in the region of small or 
zero gap was quite unexpected. They not only highlighted the fact that at v = 1/2 the system saw zero gauge field 
on average and therefore had a Fermi surface, they also argued that the best way to think of the region near v = 1/2 
was in terms of CF's seeing an effective field B* = B — B 1 / 2 - This meant the particles would bend with a radius 
R* = 2j& with kp = V^rrn , a result that was verified experimentally by Kang et al (1993), Goldman, Su, and Jain 



(1994)fand Smet et al (1996). 

HLR computed the electromagnetic response within the RPA. They showed that the system was as compressible 
as a traditional Fermi liquid. It had a longitudinal conductivity a xx = This was tested by the damping and 

velocity shift of surface acoustic waves by Willett et al (1990,1996) who also found that away from v — 1/2, there was 
a resonance in the velocity shift when the wavelength of the surface acoustic wave coincided with 2R*. 

HLR identified an overdamped mode in the density-density correlation with dispersion u> ~ iq 3 v(q). This mode, 
responsible for the compressibility of the system, also produces a (logarithmic) divergence in the CF mass by entering 
its self-energy. This in turn leads to a gap E ~ l/(plnp) for v — pj (2p +1) as p — > oo. It has been shown that the 
mass divergences do not affect bosonic (e.g., density-density) correlations 55 

These successes notwithstanding, the HLR approach had some room for improvement. It did not allow a clear 
separation of the LLL physics, i.e., it did not have a smooth m — > limit. The CS fermion, obtained by attaching 
two flux tubes to electrons had a charge of unity (and not e* =0). It did not make contact with the dipole picture 
(Read, 1994,1996) based on a wavefunction analysis. 

The Extended Hamiltonian Theory (EHT) was devised by us to address some of these issues. The key idea is that 
in order to discuss the correlation hole that accompanies the electron to form the CF in tractable form, one must 
enlarge the Hilbert space to describe collective charge degrees of freedom and to place a suitable number of constraints. 
This idea, along with expressions for the new coordinates, and a certain change of variables, all arose in our original 
tortuous route (Shankar and Murthy, 1997, Murthy and Shankar, 1998a) from electrons to CF's. In this review, we 
spared the reader details of the historical route and gave an axiomatic description (Murthy and Shankar, 2002) of the 
EHT. 

Here we begin with the electronic hamiltonian Eqn. (|262|l . and add for each electron a new pseudovortex guiding 
center coordinate R„ (whose algebra corresponding to charge — c 2 ). Thus the Hilbert space is bigger now and R„ has 
no dynamics. 

If we are interested in physics at the cyclotron scale, we can focus on the kinetic term. The mass m appears here as 
it should, in ojq. The collective coordinate formed out of rj e carries the entire Hall current. (This answers the question 
of who carries the Hall current at v = | when e* =0.) In fact they carry the Hall current even at gapped fractions 
where e* ^ 0, thus substantiating the general belief that the Hall current is not affected by disorder because it is 
carried by collective coordinates. Having identified them as our oscillators, (for which we have a specific hamiltonian) 
we have paved the way for a detailed analysis of the Hall response. We also showed how extract an effective LLL 
theory by approximately integrating out the coordinate t] e in our discussion of LL mixing. 



55 Kim, Purusaki, Wen, and Lee (1994), Kim, Lee, Wen, and Stamp (1994), Kim, Lee, and Wen (1995), Stern and Halperin (1995). 
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Let us now proceed to drop higher LL's and work with 

H{P) = V = \Y1 KQ) v(q)e- {ql)2/2 K9) (266) 

= ex P H<1 • (267) 

i 

The densities p(q) obey the magnetic translation algebra since R e is the electron guiding center coordinate in the CF 
basis. The advantage of this basis is that the velocity operator II sees a weaker field B* , leading to a nondegenerate 
HF ground state. 

The hamiltonian commutes with the operators 

x(q) = Xl exp 

3 

The pseudovortex densities x(q) thus form a closed algebra, the symmetry algebra of H . 

In the first derivation of the theory (Shankar, 1999), the following constraint naturally emerged: 

X(q)|physical state) = 0. (269) 

In the extended approach (Murthy and Shankar, 2002), H v is a cyclic coordinate with no dynamics. One is free to 
supplement the theory with the above constraint since nothing physical depends on the auxiliary coordinate R v . We 
made this choice so that there would be just one set of equations to deal with. 

We described two ways to proceed, the choice being dictated by what we want to calculate; the conserving approx- 
imation and the shortcut that used the preferred charge. 

For situations where the symmetries of H are important one uses the conserving approximation, i.e., one in which 
the constraint is respected at the level of Greens functions. This amounts to ensuring gauge invariance. We reviewed 
the compressibility paradox where gauge invariance made all the difference. The paradox concerned the system at 
v = h which Halperin and Stern (1998) argued must be compressible (as had been predicted by HLR). This result 
seemed to be at odds with our operator description in which the charge was dipolar. Halperin and Stern gave heuristic 
arguments for how a system of dipoles could still be compressible if their hamiltonian had X-invariance- invariance 
of the energy under the shift all momenta, a symmetry first noted by Haldane, and which appeared in our work as 
part of a gauge symmetry. The more detailed analysis of Stern, Halperin, von Oppen, and Simon (1999) drove the 
point home. They considered first a model where the flux tubes were spread over a distance 1/Q and the hamiltonian 
had if-invariance to sufficiently high order in Q to examine the question of compressibility. The density-density 
correlation (now of the dipole-dipole form) was computed in the RPA (which was exact in the limit Q — > 0) and the 
factors of q 2 in the numerators from the dipoles were cancelled by the exchange of the overdamped mode. They then 
showed that in the actual FQHE problem, if -invariance could be built in if the Landau parameter F\ = — 1. The 
conserving approximation was also employed by Read (1998) within the Pasquier-Haldane (1998) formalism for v = 1 
bosons (to which our theory reduces if we set c = 1 or v = |). Read (1998) established that the fermions interacted 
with each other as dipoles by exchanging a transverse collective mode, but that as q — ► 0, the propagator for this 
mode produced enough negative powers of q to overturn the q 2 coming from the dipolar factors from the two ends. 
These calculations are important in showing that if constraints are taken into account, the correct charge of the CF 
will emerge and that despite the weaknesses in the wavefunction based arguments leading to the dipole picture (see 
Section li.D.3(l . the predicted dipole moment is correct in the above sense. Murthy (2001a) was able to show that in 
the gapped fractions one could obtain structure factors in accord with Kohn's (1961) theorem (S ~ q 4 ) in a conserving 
calculation. He also showed how to compute magnetoexciton dispersions. 

In all other situations (where a gap or temperature or both suppress the deep infrared region u> ~ q 3 v(q)) the 
shortcut using /O p (q) was the weapon of choice. Here we make the replacement p — > p p = p — c 2 x, which is allowed 
in the exact theory. This choice allows one to employ naive Hartree-Fock calculations that respect Kohn's Theorem 
(S ~ g 4 ). As a bonus the CF charge and dipole moment emerge in the power series expansion of ff . It is only in this 
sense that the CF can be viewed as the union of an electron and a correlation hole of charge — c 2 . With these features 
built in at tree level, the usual approximations such as Hartree-Fock are applicable as long as the large ql region is 
avoided by the potential. 

The EHT gives a uniform and precise description of the internal structure of the CF. Whereas in the wavefunction 
based description, the Laughlin fractions allowed for a simple picture of the CF (an electron bound to 2s vortices), and 
the rest of the Jain series (upon projection) did not, in the EHT the CF is viewed as an electron plus a pseudovortex. 
Especially interesting is the case of v = h. The expansion for ]? in a power series shows that it begins with a term 
that couples to an external electric field exactly as a dipole moment of strength l 2 z x p would. This formula, in 
operator form, makes no reference to zeros of the wavefunction or vortices, neither of which is robust. It does not 
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have the problems of the wavefunction-based dipole picture. These problems arose upon antisymmtrization, which we 
carry out by expressing p p in second quantized form in terms of fermionic operators, we saw that the place look for 
dipoles is not in the wavefunction, but in correlation functions at high frequencies at and above the CF Fermi energy. 

The operator approach gives a concrete realization of one of the primary expectations in FQHE: Once the kinetic 
energy is quenched by restricting electrons to the LLL, it will be resurrected by interactions and that this low-energy 
problem will be characterized by one common scale, the electron-electron interaction. Eqn. Ij61|l which gives H as a 
quadratic function of p, (or p p if we use the preferred combination) embodies all these expectations. 

The HF approximation to H(p p ) was used to compute transport gaps (Murthy and Shankar, 1999, Shankar, 2001) 
good to within 10 — 20% for potentials that vanish rapidly with ql. It could be used to explore magnetic transitions 
from one quantized value of magnetization to the next (Shankar, 2000, 2001). In the absence of disorder it is clear 
that all physical quantities pertaining to the FQHE (restricted to the LLL) are functionals of the potential v(q). For 
the Zhang-Das Sarma potential which we use for most of our work, this means a function of A. We saw that given a 
value of A from experiment, we did not need several masses for several phenomena, they all came from one potential 
(Shankar, 2000, 2001). This is because the hamiltonian for the CF's is rather unusual 56 , and contains the kinetic and 
potential terms in a monolithic form. Hidden in it are the various mass scales m a and m p appropriate to various 
phenomena (activation, polarization), all functionals of the interaction v(q). 

The operator approach also clarified the question of whether or not CF's are free. Given the prominent variations 
in the magnetoexciton dispersions and the fact that it takes two very different masses to describe polarization and 
activation, it is clear that they are not. But why do they appear to be free for some magnetic phenomena? Our 
theory (Shankar, 2000, 2001) shows that it is an accident coming from rotational invariance and d = 2. 

The EHT allows us to compute physical quantities at T > 0, such as polarization and relaxation rates for gapless 
states (Shankar, 2000, 2001). In the experiments of Dementyev at al (1999), a A determined from one polarization 
data point gives the polarization and relaxation curves for two tilts and a range of temperatures. This is to be 
contrasted with attempts to fit the data with a mass and interaction pair (m, J), where the four curves require four 
nonoverlapping values of these pairs. In the case of Melinte et al (2000), Freytag (2001), the 1/Ti predictions (that 
vary over orders of magnitude) are off by a factor of 2, but the polarization data are well described by a single A. 
Once again the success can be traced back to the fact that the CF hamiltonian we use is of a nonstandard form 
parametrized by a A, which in turn determines all the mass and energy scales relevant to each given process. 

When calculating the polarization of gapped states, it turns out to be essential to take into account spin-waves for 
spontaneously polarized cases. This is done (Murthy, 2000c) by mapping the low-energy dynamics of the problem 
on to the continuum quantum ferromagnet treated in the large-iV approximation (Read and Sachdev, 1996), the 
parameters of which are extracted from the HF treatment of the problem in the hamiltonian theory. The results 
(Murthy, 2000c) agree extremely well with experiments (Khandelwal et al, 1998, Melinte et al, 2000) up to very high 
temperatures. This calculation, like all the others in the Hamiltonian theory, is in the thermodynamic limit and free 
from finite-size effects. 

The EHT allows us to calculate the gaps of the Wigner crystal in terms of CF's (Narevich, Murthy, and Fertig, 
2001). These gaps are off by a factor of 2 when compared to experiment (Jiang et al, 1990, Jiang et al, 1991). This 
should be contrasted to previous approaches in which the gaps were off by two orders of magnitude. This approach 
also allows us to consider inhomogeneous states (Murthy, 2000a) with polarizations not allowed by the CF theory of 
homogeneous ground states. 

Hopefully we have succeeded in establishing that the hamiltonian theory of the FQHE is a comprehensive scheme 
for addressing and answering a variety of qualitative and quantitative questions, for gapped and gapless states, at 
zero and nonzero temperatures. 

We have given the reader a taste of what can be done using this formalism. While many things have been clarified, 
there are many open problems where we believe this approach may be fruitfully applied. 

One outstanding open problem is that of computing transport coefficients in the quantum Hall regime from a 
microscopic theory. By identifying the collective mode which carries the Hall current, we have set the stage for a 
study of transport in the presence of disorder. 

The FQH edge 57 is another open problem where the hamiltonian approach is applicable. After Wen's ( 1990a, b, 
1992) description of the edge as a chiral Luttinger liquid, other descriptions have appeared based on wavefunction and 
field-theoretic approaches 58 . The field-theoretic descriptions are all effective theories whose connection to the electron 
problem has not been rigorously established. While the wavefunction approaches are microscopic, it is impractical 
to calculate time-dependent response functions in them. As we have seen in the extended formalism, we have an 
exact rewriting of the microscopic electron problem, but with the added advantage of a nondegenerate starting point. 



Consider for example Eqn. 11541 for v = 1. 

For an excellent review, see Kane and Fisher (1997). 

For wavefunction-based approaches, see: Ziilicke, U., and A.H.MacDonald (1999), Goldman and Tsiper (2001), Mandal and Jain (2001b). 
For field-theoretic descriptions, see: Lee and Wen (1998), Lopez and Fradkin (1999), Levitov, Shytov, and Halperin (2001). 
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The calculation of edge reconstructions (MacDonald, Yang and Johnson, 1983, Chamon and Wen, 1994) in CF-HF 
and the number and dispersions of the edge collective modes in TDHF seem to be very accessible in the hamiltonian 
approach. One of the interesting things to consider in the edge problem is tunneling 59 . For this one needs a description 
of electron creation and destruction operators within the CF basis, which is an open problem. 

Since there is a gap, the effects of disorder on the Jain series ought to be describable in some simple approximation, 
such as the self-consistent Born approximation for single-particle properties. Such a treatment seems to capture many 
of the experimental facts at v = 1 (see Murthy, 2001b), such as the reduction of the transport gap due to disorder, 
the variation of the transport gap as a function of Ez (see, for example, Schmeller et al (1995)), the polarization at 
v = 1 60 , etc. Such an approach ought to be applicable to the gapped fractions. The effect of disorder on excitons is 
also interesting, and can be treated in a "self-consistent exciton approximation" (Kallin and Halperin 1985). 

A more ambitious problem that might benefit from the hamiltonian approach is the question of whether the 
v = \ state remains metallic at T = 0, when disorder is included. The CS formalism of HLR (with its logarithmic 
corrections) suggests that metallic behavior disappears. It is known that noninteracting fcrmions are always localized 
in two dimensions, regardless of how weak disorder is (Abrahams, Anderson, Licciardello, and Ramakrishnan, 1979). 
By extension, it seems plausible that a Fermi liquid, which is adiabatically connected to noninteracting fermions, 
should also be an insulator on the longest length scales. The v = \ system is rather unique, in that the state is 
produced by interactions. The Extended hamiltonian has additional symmetries absent in the zero-field problem. 

We have excluded from this review many interesting topics to which the hamiltonian is applicable, such as double- 
layer systems 61 , paired states 62 , etc. These omissions reflect the double constraints of limited space and our own lack 
of expertise. 
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APPENDIX A: Matrix elements 

Many of the calculations performed in this paper deal with the preferred density p p . In second quantization we 
write it as 

W{*l) — ^ ' d\ n2n2 d mini p m2n2 .. rnini (Al) 

m2ii2;mi)ii 

where d) m2n2 creates a particle in the state \m,i n-z) where m is the angular momentum and n is the LL index. They 
are related to the CF cyclotron and guiding center coordinates, R and r] as follows. Let 

_ R x — iR y | _ R x + iR y CA2) 

where I* = lj\J\ — c? is the CF magnetic length. These obey the oscillator algebra 

M f ]=l (A3) 



The theory was described in a series of beautiful papers by Kane and Fisher (1992a, b, 1994). Recently there have appeared some truly 
remarkable tunneling experiments where the data span several orders of magnitude of voltage: Grayson et al (1998), Chang et al (2001) 
Tycko et al (1995), Barrett et al (1995), Aifer et al (1996), Manfra et al (1996) 

For some of the early theoretical references see: Fertig (1989), Wen and Zee (1992, 1993), Ezawa and Iwazaki (1993), Yang et al (1994), 
Moon et al (1995). For a recent review, see Girvin and MacDonald (1997): For recent experiments in double-layer systems see: Spielman 
et al (2001), Kellogg et al (2002) 

The experiment that stimulated this subfield was the observation of a plateau in the Hall resistance at v = |: Willett, Eisenstein, 
Stormer, Tsui, Gossard, and English (1987). For some of the theoretical references, see: Haldane and Rezayi (1988), Moore and Read 
(1991), Greiter, Wen, and Wilczek (1992), Ho (1995), Park, Mclik-Alavcrdian, Boncstccl, and Jain (1998), Morf (1998), Bonesteel 
(1999), Read and Green (2000). 



given 

[R x , R y ] = —il* 2 . 

Similarly we define, in terms of the cyclotron coordinates, 

Vx + i-Vy + '/• _ "i; 
a = — ■ a' = ; 

V2/* 5 VW 1 



which obey the oscillator algebra 



given 



The states \mn) are just the tensor products 



where |00) is annihilated by both a and b. 
We will now show that 



[a , a 1 ] = 1 
(tf) m (at)" 



ml vn! 



100) 



G4 

(A4) 

(A5) 

(A6) 
(A7) 

(A8) 



(m 2 |e- ic »- R |mi) 



m 2 ! _-x/2 f- i( l+ l 



mi'. 



* \ mi—m-2 



V2 



where 

x = q 2 l* 2 /2, q± = q x ± iq y 
L is the associated Laguerre polynomial, and m\ > ni2. If mi < m 2 one may invoke the relation 

(m 2 \e-^ R \ mi ) = (m 1 \e+^' R \m 2 )* . 
Likewise to establish Eqn. (|A9|) . consider the coherent states 



: |o> = E 



with the inner product 



First we write from the definitions given above 

(^|e- iq - R |z) 



oo oo 

E E- 

mi— m2— 

R(z,z,q). 



zm 2 mi 



lm, 2 \ \Jm\\ 



{male-^lmi) 



On the other hand 



(z\e-^ R \z) = (z\ 



il* 



exp(--y=(q+b i +q-b))\z) 



il* 



il* 



exp 
R(z, z, q) 



zz- -^( z i- +1+ z ) 



-g 2 r 2 /4 



(A9) 
(A10) 
(All) 

(A12) 

(A13) 



(A14) 
(A15) 

(A16) 

(A17) 

(A18) 
(A19) 
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Comparing Eqns. (|A14IA18(I and matching powers of z a z b we obtain Eqn. I|A9(1 if we recall 



nil • 



(-1)* 



t ^(m 2 -ty.(m l -m 2 + t)\ t\ 



Likewise, to establish 



{n 2 \e 



(A20) 



(A21) 



(again for n\ > n 2 ) we just need to remember that the commutation rules of the components of r\ have a minus sign 
relative to those of R, which exchanges the roles of creation and destruction operators and hence q+ and q- . 

Now we consider matrix elements of p,x, P P ■ As a first step, let us express the operators R e and R v in terms 
of CF guiding center and vortex coordinates R and r\. We have seen that in the CF representation 



R P = r-l 



,z x n 

l + c 



R + 77C 



(A22) 



if we recall I 2 = l* 2 (l -c 2 ). 
It can similarly be shown that 



Thus in first quantization 



R„ = r + r 



z x n 

c(l + c) 



= R + r]/c. 



= 2 = 
P = P ~ c x 

P = ex P(- ?; q ■ R i) exp(-iq ■ T] lC ) 

i 

X = Ycxpi-iq-R^expi-iq-rii/c) 



(A23) 

(A24) 
(A25) 

(A26) 
(A27) 



Armed with Eqns. l|A9l and lA21l ) we may finally write for the matrix elements of p(q) defined in Eqn. (|A1() . 



ni2tL2\mini 



m2[ e - x /2 f-k+l 
mil 




v/2 

I / —icq_l x 



ni! V \/2 



e- xc2 l 2 LH- n i{xc 2 ) 



Superscripts on p™ 2TOl and Pn 2ni which will be apparent from the subscripts, will usually be suppressed. 



APPENDIX B: Hall response in the extended picture 

To compute the DC Hall conductance we just need the zero momentum component of the current operator 

j j 

Note that the current at q = depends only on the cyclotron coordinate, just as it depended only on the oscillator 
coordinate in the small treatment. Upon coupling the system to an external potential 3>(q), H becomes 

H = E^i+ e ^)E^ ( B2 ) 

j 3 
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If we now keep the interaction to 0(q), and recall r e = R e + rj e , we obtain 

H = E + e<i>(q) E^ iq ' ^ + M + • ■ ■ (B3) 

If we complete squares on T/gj, we can read off its mean value in the presence of $, or the corresponding electric field 
iq$. From this we obtain a mean current J(0) corresponding to the right Hall conductance of ne/B. 

The exercise should make it clear that the oscillator coordinates are just the collective coordinates formed from 
rj e at small q. If we go to higher orders in q we will find that the current involves both rj e and R e . Our extended 
formalism allows one to explore corrections due to this mixing in a small-g expansion. 

Note that in this approach the magnetic moment of Simon et al need not be put in by hand since r\ + e is still; in 
the picture and can respond to a slowly varying magnetic field by a changing zero-point energy. 

The objections of Lee et aZ(1999) on the CF Hall conductance are moot since we do not add resistivities as in Eqn. 
(|107fl but rather conductivities (of the r) e and CF variables). 



APPENDIX C: Proof of Hartree-Fock nature of trial states 

Consider 



(f\H\i) = (p\d f Hdl\p) 



(CI) 



where |p) stands for the (ground) state with p-filled LL, and i, f label single-particle excitations on top of this ground 
state. We want to show that this matrix element vanishes if i ^ /, i.e., the hamiltonian does not mix these putative 
HF particle states. (This result was established for the small q theory by Murthy, 1999). The proof, which relies on 
just the rotational invariance of the potential, applies with trivial modifications to the hole states, i.e., to 

(pM^lp). (C2) 

This matrix element in question takes the schematic form 

(C3) 

measure: 



(f\H\i) = / (p\d f d\d 2 4d i 4\p)p 12 ( q )p 34 (-q) 

Jq 

where 1 stands for mini and so on, and / stands for an integral over a rotationally invariant 



(C4) 



Now we use Wick's theorem and perform pairwise contractions on the vacuum expectation value, bearing mind 
that 

• We cannot contract the indices 1 and 2 or 3 and 4 since this will require that q = at which point the measure 
(which contains the potential) vanishes. 

• If we contract i and / we already have the desired result. 
Here is a representative of the contractions we can get: 

Pi2P34<5/i(l - nf ) 623(1 - rcf)&M(l - nf ). 
where nf is the Fermi function for the LL labeled by m 



nf = 8(p — 1 — m) 

and so on. Since / = 1 and i = 4, the factor (1 — nf^l — nf ) = 1. The integrand assumes the form 



(C5) 
(C6) 



OO OO 



E E P/2(<j)P2i(-q) = 

m 2 =0 n 2 =p 



^2 Pm/m2(q)Pm 2 m I (-q) 
n 2 =0 



X! / , ™/«2(q)pn 2 ™,(-q) 



where we have also used the fact that e _lq R • e lq R = I (the identity operator) in doing the sum over m-2, and 

is some rotationally invariant function. It follows that every term in the sum over n 2 vanishes unless nf — rii due to 

the angular integral in q. 
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APPENDIX D: On the conserving nature of TDHF 

We will now verify that the constraint is a left-eigenvector of Ti, as required of the conserving approximation. 
To this end we will need the matrix elements 

p ni „ 2 (q) = <n 1 \e~ lc ^\n 2 > (Dl) 

W 2 (q) = <nx\e-^\n 2 > (D2) 

Note that only the cyclotron part of the electron and pseudovortex coordinates appear in these exponentials. An 
important result we will need is 

X]Xn 1 n(qi)Pnn 2 (q 2 ) = < m le^^e^ * ^\n 2 > (D3) 
n 

= e -^QixQ 2 < ni | e -*(T L +^2)-i7| n2 > ( D4 ) 

Q = ql* (D5) 

where we have used the completeness of the states \n > and the commutation rules obeyed by rj. Finally we separate 
the exponentials in the reverse order to get the second useful identity 

X]Xmn(qi)p„„ 2 (q2) 

n 

= e- 4QlxQ2 ^ / 5„ 1 „(q 2 )Xn„ 2 (qi) (D6) 

n 

Note that p and \ do not commute, even though p and x do. 
Let us now right-multiply the putative left-eigenvector x by TL 

^2 Xn' 1 n' 2 {q)'H{n' 1 n' 2 - 1 n l n 2 -,q) = 

™i™2 

(e(rci) - e(n 2 ))xrun 2 (q) 
+ ^(^(^)-^K))^^e^ 2i2 / 2 x„;^(q)^ ni (-q)p ni „ 2 (q) 

"i™2 

- 2^( n f M -N F K))y ^u(s)e 2 P„ 1 „' 1 (s)xn' 1 r l ^(q)/5^„ 2 (-s)e 



sxq 



(D7) 



Let us consider the direct and exchange terms separately. In the direct term, one n' index can always be summed 
freely, while the other is constrained by the Fermi occupation factor Np. The sum over the free n' gives, according 
to Eqn. (|T)3|) 

direct term = ^ N F (n 2 )p niJl2 (q) < n' 2 \e~ i{ = ~ c) ^\n' 2 > 

n' 2 

- ^^K)p nin2 (q) < nile-^^-^^Ini > (D8) 

The two terms are immediately seen to cancel. Now let us turn to the exchange terms, and consider the one that 
has the factor Np(n' 2 ), and a free sum over n'\. In this term, one can use Eqn. (|D6(1 to exchange the p and x matrix 
elements to obtain 

exchange term = - ^ N F {n' 2 ) ^ x„ in / (q) x J A^. u ( s ) e - i 5- / 5 n , n , (s) / 5 n / nj (-s) (D9) 



Notice that the phase factor e 4 ^ ' sxq has been cancelled by an opposite phase factor from Eqn. I|D6(I . Now the 
angular s integral forces rt^ = n 2 for a rotationally invariant potential, and the result contains the Fock energy of the 
state n 2 

e F (n 2 )x nin2 (q) (D10) 

Similarly, the other exchange term proportional to NF(n'i) ends up giving — t F {ni)x nin2 (q). Due to the peculiar 
nature of the Hamiltonian the Hartree energy is a constant independent of the CF-LL index, and the difference of 
the Fock energies is the same as the difference of the full HF energies. Thus, the exchange contributions cancel the 
diagonal term (e(ni) — e(ri2))xnin 2 (q), and X«in 2 (q) is indeed an left eigenvector with zero eigenvalue for 7i. Also, 
this property is independent of the form of v{q) as long as it is rotationally invariant. 
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APPENDIX E: Activation gaps 

Now we need to find the energy cost of producing a widely separated particle-hole (PH) pair. This will be done by 
evaluating 



A a = (p + P\H\p + P) + (p + H\H\p + H) 
-2<p|tf|p) 



= E(P)+E{H) 



I, 



1 f d 2 q _ 2,2 /2 

2 / ^ <q)e ■ 



(El) 
(E2) 

(E3) 



where P denotes a particle added to the state labeled p, = (n = p, m = 0) and H denotes a state in which a hole has 
been made in the state p = (n = p — 1, m = 0). Let us consider 



E{P) = (p|d„d t 1 d 2 4d 4 dt|p) 

In performing the contractions we 

• Do not make any contractions within H. This gets rid of E = (p\H\p), the ground state energy 

• Do not contract 1 with 2 or 3 with 4 since v(0) = 0. 
We end up with 



(E4) 



/ [^i£ 2 3<Ml-nf)(l-nf)(l-n; 
J 1 



- ^3^14^(1 - n 3 )(n 4 )(1 - n 2 )] p 12 psA 

Since 4 = /i = 1 in the first term, we can drop (1 — nf )(1 — nf) and for similar reasons (1 — n^ )(l — ) in the 
second giving us 



00 00 



E ^ P ) = E E Mq)/M-q) 

m2=0 n 2 =p 
00 p— 1 

-EE ^(qK 2 (-q)- 

rri2— U2—0 



(E5) 
(E6) 



Since the sum over m 2 is unrestricted, we can use completeness and e 4q R • e tq R = / to get rid of the m-index 
altogether. Thus we end up with 



e{p) = Em 2 -Em 2 

\ra=p n=0 / 

p-1 

(P|p(9)p(-9)|p)-2^|p p „| 5 



n=0 



(E7) 
(E8) 



A similar calculation for the hole state gives (upon dropping the ground state energy as usual) 



E(H) 



p-i 



-<p - l\ P (q)p(-q)\p -1)+2J2 \P P -h 



n=0 



where 



(n\p(q)p(-q)\n) = ^ Mq)t. 

n' — o 

Putting all the pieces together, we obtain the gap. 



(E9) 



(E10) 
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APPENDIX F: Critical fields for magnetic transitions 

We need to calculate 

E(p — r, r) = (p — r, r\H\ p — r, r) 



(Fl) 



the energy in a state with p — r spin-up LL's and r spin-down LL's. Since the HF calculation for the spinless case is 
very similar, this treatment will be brief. We write 



H = ^2 / d\d 2 dld4,pi2p34 

1234"^ 



(F2) 



with the understanding that a label like 1 stands for the triplet (m,mi, s\), s being the spin. The matrix elements 
pij are defined by 



P\2 



le 



— iq-R, 



and as a result 



E(p — r,r) — 



f nf( S )(l-nf( S ))|p nin J 2 £<m|I|: 



=n/p 

where we acknowledge the fact that the occupation factors nf and can depend on the spin. We have also used 
the fact that the sum over all values of m is the degeneracy of each CF-LL, n/p. Carrying out the sums over n\ and 
«2, we obtain 



n 

Pjg 
r-1 



-r-1 



E{p — r,r) = 

"p—r—1 

{ni\p(q)p(-q)\ni) - ^ \Pnm 2 

ni—O ni,ri2—0 
r-1 

+ (ni\p(q)p(-q)\n 1 ) - ^ \Pnm, 

ni ,rt2— 



It is now straightforward to compute the critical field for the transition |p — r, r) — * |p — r — 1, r + 1) by invoking 



E(p — r, r) — E(p — r — 1, r + 1) = g 



2m e 



B c 



(F3) 



APPENDIX G: Acronyms and Symbols 

This paper invokes many symbols and acronyms, not all of which have been standardized. For the convenience of 
the reader, a list is supplied in Table IjVIfl . 
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